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In this review we summarize theoretical progress in the field of active matter, placing it 
in the context of recent experiments. Our approach offers a unified framework for the 
mechanical and statistical properties of living matter: biofilaments and molecular motors 
in vitro or in vivo, collections of motile microorganisms, animal fiocks, and chemical or 
mechanical imitations. A major goal of the review is to integrate the several approaches 
proposed in the literature, from semi-microscopic to phenomenological. In particular, 
we first consider dry systems, defined as those where momentum is not conserved due 
to friction with a substrate or an embedding porous medium, and clarify the differences 
and similarities between two types of orientationally ordered states, the nematic and the 
polar. We then consider the active hydrodynamics of a suspension, and relate as well 
as contrast it with the dry case. We further highlight various large-scale instabilities of 
these nonequilibrium states of matter. We discuss and connect various semi-microscopic 
derivations of the continuum theory, highlighting the unifying and generic nature of 
the continuum model. Throughout the review, we discuss the experimental relevance 
of these theories for describing bacterial swarms and suspensions, the cytoskeleton of 
living cells, and vibrated granular materials. We suggest promising extensions towards 
greater realism in specific contexts from cell biology to ethology, and remark on some 
exotic active-matter analogues. Lastly, we summarize the outlook for a quantitative 
understanding of active matter, through the interplay of detailed theory with controlled 
experiments on simplified systems, with living or artificial constituents. 
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I. INTRODUCTION 

The goal of this article is to introduce the reader to a 
general framework and viewpoint for the study of the me- 
chanical and statistical properties of living matter and of 
some remarkable non-living imitations, on length scales 
from sub-cellular to oceanic. The ubiquitous nonequilib- 
rium condensed systems that this review is concerned 
with (Joanny and Prost, 2009 Jiilicher et at. 2007 



Toner et al. 2005 Vicsek and Zafeiris 2012) have come 



to be known as active matter (Ramaswamy, 2010). Their 
unifying characteristic is that they are composed of self- 
driven units - active particles - each capable of convert- 
ing stored or ambient free energy into systematic move- 



ment (Schweitzer, 2003). The interaction of active par- 
ticles with each other and with the medium they live 
in gives rise to highly correlated collective motion and 
mechanical stress. Active particles are generally elon- 
gated and their direction of self-propulsion is set by their 
own anisotropy, rather than fixed by an external field. 
Orientational order is thus a theme that runs through 
much of the active- matter narrative, as can be seen for 
instance in the image of a swarm of myxobacteria, shown 
in Fig. [l] The biological systems of our interest in- 
clude in vitro mixtures of cell extracts with bio-filaments 
and associated motor proteins (Fig. [2|, the whole cy- 
toskeleton of living cells, bacterial suspensions (Fig. [3|, 
cell layers (Fig. [I]), and terrestrial, aquatic (Fig. [5| and 
aerial flocks. Non-living active matter arises in layers of 
vibrated granular rods, colloidal or nanoscale particles 
propelled through a fluid by catalytic activity at their 
surface (Fig. [6|, and collections of robots. A distinctive 
- indeed, defining - feature of active systems compared 
to more familiar nonequilibrium systems is the fact that 
the energy input that drives the system out of equilib- 
rium is local, at the level of each particle, rather than at 
the system's boundaries as in a shear flow, for example. 
Each active particle consumes and dissipates energy go- 
ing through a cycle that fuels internal changes, generally 
leading to motion. Active systems exhibit a wealth of 
intriguing nonequilibrium properties, including emergent 
structures with collective behavior qualitatively different 
from that of the individual constituents, bizarre fluctua- 
tion statistics, nonequilibrium order-disorder transitions, 
pattern formation on mesoscopic scales, unusual mechan- 
ical and rheological properties, and wave propagation and 
sustained oscillations even in the absence of inertia in the 




FIG. 1 (color online) Liquid-crystalline order in a myxobacte- 
rial flock. Image courtesy of Gregory Velicer (Indiana Univer- 
sity Bloomington) and Juergen Bergen (Max-Planck Institute 
for Developmental Biology). 
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FIG. 2 Patterns organized in vitro by the action of mul- 
timeric kinesin complexes on microtubules, imaged by dark 
field microscopy. The concentration of motor proteins in- 
creases from left to right. Image (a) shows a disordered array 
of microtubules. The other two images display motor-induced 
organization in spiral (b) and aster (c) patterns. The bright 
spots in the images correspond to the minus end of micro- 
tubules. These remarkable experiments from [Surrey et al.\ 



(2001 ) led the way to the study of pattern formation in active 



systems. Adapted with permission from [Surrey et al. (2001). 



strict sense. 

Living systems of course provide the preeminent exam- 
ple of active matter, exhibiting extraordinary properties 
such as reproduction, adaptation, spontaneous motion, 
and dynamical organization including the ability to gen- 
erate and to respond in a calibrated manner to forces. 

A theoretical description of the general properties of 
living matter is not currently achievable because of its 
overall complexity, with the detailed state of a cell deter- 
mined by a hopelessly large number of variables. How- 
ever, in a given living organism there are at most 300 
different cell types, which an optimist could view as a 
very small number given the immensity of the accessible 
parameter space. Perhaps, therefore, global principles 
such as conservation laws and symmetries constrain the 
possible dynamical behaviors of cells or, indeed, of organ- 
isms and populations, such as collections of bacteria (Fig. 
|3|, fish schools (Fig. [s]) and bird flocks. Quantifying the 
spontaneous dynamical organization and motion of living 



FIG. 3 Bacterial "turbulence" in a sessile drop of Bacillus 
subtilis viewed from below through the bottom of a petri dish. 
Gravity is perpendicular to the plane of the picture, and the 
horizontal white line near the top is the air- water-plastic con- 
tact line. The central fuzziness is due to collective motion, 
not quite captured at the frame rate of l/30s. The scale bar 
is 35/i m. Adapted with permission from |Dombrowski et aL\ 
(2004| 




FIG. 4 A disclination defect of strength m — —1/2 formed 
by human melanocytes on a plastic surface. The bars are 
lOO/xm. The three images show three different situations: (i) 
the core of the defect is an area free of cells (left); (ii) the core 
of the disclination is an area with isotropically distributed 
cells (center); and (iii) the core of the defect is occupied by 
a star-shaped cell (right). The cells that form the nematoid 
order are in an elongated bipolar state . Reproduced with 
permission from Kemkemer et al. (2000). 



systems is a first step toward understanding in a generic 
way some of these principles, by focusing on specific ques- 
tions that are accessible to theory. This has proved 
to be the case for the long wavelength behavior of ac- 



tive membranes ( Manneville et 



msma 



et al. 



2001 



Prost and Bru- 



|1996[ [Ramaswamy and Rao[ |2001[ Ramaswamy 



2000[ ?); the general theory of flocking ( [Toner and 



Tu 1995 1998), (Toner et a/., 2005); or the macroscopic 



mechanical properties of the cytoskeleton as an active 
gel dJiilicher et~aL\ |2007| [Kruse al] |2005|. Agent- 



based models offer a minimal approach to the study of 
active systems, with an emphasis on order and fluctua- 
tions rather than forces and mechanics. In such a setting, 
seminal studies of flocking as a phase transition were first 



carried out by Vicsek et al. ( 1995 ) [see also the remark- 
able computer- animation work of Reynolds (|1987)], with 



important modifications and extensions by 



Chate ( 


2004) and 


Chate et al. ( 


2008a 2006) 



Gregoire and 



FIG. 5 (color online) A remarkable demonstration of 
polar order in a sardine school (Courtesy of Jon 
Bertsch, from underwater images form the Sea of Cortez: 
htt p : / / www . t halassagr aphics . com /blog / ?p= 1 6 7) . 




FIG. 6 (color online) Colloids that swim: Hydrogen peroxide 
is oxidized to generate protons in solution and electrons in 
the wire on the platinum end of micron-sized gold-platinum 
wires. The protons and electrons are then consumed with the 
reduction of H2O2 on the gold end. The resulting ion flux 
induces an electric field and motion of the particle relative 
to the fluid, propelling the particle towards the platinum end 
at speeds of tens of microns a second. On longer timescales 
the rods tumble due to thermal or other noise; the resulting 
motion is thus a persistent random walk. Reprinted with 
permission from |Paxton et aL\^00^ . 



els describe point particles with fixed speed moving on 
an inert background. The direction of motion changes 
according to a noisy local rule that requires particles to 
align with their neighbors at each time step. This family 
of models displays a well-defined transition from a disor- 
dered to an ordered phase with decreasing noise strength 
or increasing density. In the context of the cytoskele- 
ton activated by motor proteins, detailed simulations of 
ensembles of semi-fiexible filaments on which motor bun- 
dles can exert force dipoles have also been carried out by 



several authors (Head et a/., 2011 


Mogilner and Oster 


1996 


Nedelec et al. , 


1997 Pinot et al. , 2009 


). 



In this article we will not review the agent-based mod- 
els and the wealth of results obtained by numerical sim- 
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FIG. 7 (color online) Measurement of contractile stresses 
in cell (Xenopus oocyte) extracts and a reconstituted acto- 
myosin gel made using skeletal muscle myosin II thick fila- 
ments: (a) Schematic of active gel in a narrow glass cap- 
illary (inner diameter, 400/im) sandwiched between two oil 
droplets, (b) Confocal fluorescence images of the network at 
five different time points shows that the active gel first pulls 
away from the capillary walls and deforms the oil droplets 
as it contracts before reaching a breaking point when the gel 
collapses. ( c) Dark- field images of contracting cell extract 
shows the same behaviour. Analysis of the deformation of 
the oil droplet gives an estimate of the contractile force which 
is around 1/jN or lOOpN per act in bundle. Adapted with 
permission from Bendix et al. ( 2008 ) . 



ulations, but rather focus on identifying generic aspects 
of the large scale behavior of active systems and charac- 
terizing their material properties. Many of the macro- 
scopic properties of active systems are universal in the 
sense that systems operating at widely differing length 
scales, with significant differences in their detailed dy- 
namics at the microscopic level, display broadly similar 
properties. Visually similar flocking phenomena are seen 



in fish shoals (Parrish and Hamner 1997) and collections 



of keratocytes (Szabo et a/., 2006). Contractile stresses 



are evident on a sub-cellular scale in the cytoskeleton 
(Fig. [7|) dBendix et al] |2008| [joanny and PrQst[ |2009| ), 
as well as on a scale of many cells in swimming algae 
(Rafai* et al. 2010). The hope is to be able to classify 



active matter in a small number of universality classes, 
based on considerations of symmetry and conservation 
laws, each with a well-defined macroscopic behavior. We 
consider here four classes of active matter according to 
the nature of the broken symmetry of the ordered phase 
and the type of momentum damping. First the broken 
symmetry: elongated self-propelled objects are in gen- 
eral polar entities with distinct heads and tails, which 
can cooperatively order either in a polar (ferromagnetic) 
phase or in a nematic phase. In a polar phase, all the 
microscopic objects are on average aligned in the same 
direction; this is the case for bacteria or fish. The polar 
order is described by a vector order parameter p, known 
as the polarization. Nematic ordering can be obtained in 
two ways, either in systems where polar self-propelled ob- 
jects are parallel but with random head-tail orientations 
or in systems where the self-propelled particles are them- 





nematic 


polar 


dry 


melanocytes 
vibrated granular rods 


migrating animal herds 
migrating cell layers 
vibrated asymmetric 
granular particles 


wet 


bipolar catalytic rods in 
suspension 


cell cytoskeleton 

cytoskeletal extract in 
bulk suspension 

swimming bacteria in 
bulk 

Au-Pt catalytic colloids 



TABLE I Examples of active systems classified according to 
symmetry and presence/absence of momentum conservation 
(wet/dry). 



selves head-tail symmetric such as the melanocytes that 
distribute pigment in the skin. The polar /apolar distinc- 
tion in the context of living matter is brought out clearly 
by Gruler et al. (1999) and Kemkemer et al. (2000), as 
shown in Fig. |4] Nematic order is described by a ten- 
sor order parameter, Q, known as the alignment tensor. 
Next, the nature of damping: in a bulk fluid, active or 
otherwise, viscosity damps the relative motion of neigh- 
boring regions, which means that the total momentum 
of the system is conserved. In systems in contact with a 
substrate, e.g., restricted to lie on a surface or between 
two closely spaced walls, or moving through a porous 
medium, the drag due to the substrate dominates, ei- 
ther simply through no-slip or through binding and un- 
binding of constituents of the systems to the confining 
medium, and the momentum of the fluid is not conserved. 
Table [l| summarizes various examples of active systems 
classified according to symmetry and presence/absence 
of momentum conservation (wet /dry). Of course even 
in "dry" systems, where momentum is not conserved, 
hydrodynamic or medium-mediated interactions may in 
some cases be important, depending on the length scale 
of interest (Schaller et a/. , 2011a). For instance, for active 



particles in a viscous fluid of viscosity r] and experiencing 
a frictional drag 7, hydrodynamic flows can be neglected 
only on length scales larger than y^r]/j. 

A useful theoretical framework to describe the macro- 
scopic properties of active matter is provided by the 
methods of nonequilibrium statistical mechanics. In a 
generalized hydrodynamic approach, a coarse-grained de- 
scription of the large-scale, long-time behavior of the sys- 
tem is given in terms of a small number of continuum 
fields. The evolution of these fields is written in terms of a 
set of continuum or hydrodynamic equations that modify 



the well-known liquid crystal hydrodynamics ( ,de Gennes 



and Prost , 1993 Martin et al. , 1972 ) to include new 



nonequilibrium terms that arise from the activity. Gener- 
alized hydrodynamic theories have been very successful 
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in the description of many condensed matter systems, 



such as superfluids ( Dzyaloshinskii and Volovick 



hquid crystals (Martin et al, 1972), polymers 
1993), as well as of course simple fluids. 



1980), 



Milner 



One approach to obtain a hydrodynamic theory of ac- 
tive systems is to start from a microscopic model and use 
the tools of statistical physics to coarse-grain the model 
and obtain the long-wavelength, long-time scale equa- 



tions ( 


Ahmadi et at. 


2005 


2006| 


Aranson and Tsimring , 


2005 


Kruse and Jiilicher 


20031 


Liverpool 2003). This 



task is diflicult if the microscopic description is realistic 
and it can only be carried out at the cost of approxima- 
tions, such as low density or weak interactions. It allows 
one to relate the parameters in the macroscopic equa- 
tions to speciflc physical mechanisms (albeit in a model- 
dependent manner) and to estimate them in terms of ex- 
perimentally accessible quantities. The low-density limit 
has been worked out for several models, initially without 
and later including the effect of the embedding solvent. 
Some speciflc examples are discussed in section [V} 

An alternative, more pragmatic, approach is to directly 
write hydrodynamic equations for the macroscopic flelds 
including all terms allowed by symmetry, as was pio- 
neered for dry flocks by Toner and Tu (1995 1998), and 



extended to the case of self-propelled particles suspended 



in a fluid by Simha and Ramaswamy (2002a) and more 



generally to active-fllament solutions (Hatwalne et al 
[2004| ) . As one might expect in hindsight, novel terms, 
of a form ruled out for thermal equilibrium systems, ap- 
pear in these "pure-thought" versions of the equations of 
active hydrodynamics. 

A more systematic implementation of the phenomeno- 
logical hydrodynamic approach is to treat the nonequilib- 
rium steady state of an active system as arising through 
the imposition of a non- vanishing but small driving force 
on a well-deflned parent thermal equilibrium state whose 



existence is not in question pe Groot and Mazur 1984). 
For example, in the biological case of systems composed 
of cytoskeletal fllaments and motor proteins, such as the 
cell cytoskeleton, the driving force is the difference A/i 
between ATP and its hydrolysis products. If A/i is as- 
sumed to be small, the macroscopic hydrodynamic equa- 
tions can be derived in a systematic way following the 
Onsager procedure. One identifles thermodynamic fluxes 
and forces and writes the most general linear relation be- 
tween them that respects the symmetries of the problem. 
The novel terms mentioned above that arose by directly 
writing down the hydrodynamic equations of active mat- 
ter are then seen to be a consequence of off-diagonal On- 
sager coefficients and an imposed constant nonzero A/i. 
Generalized hydrodynamic theories have been very suc- 
cessful in the description of complex and simple fluids 



(Martin et al. , 1972 ). The advantage is that the equations 
are expanded around a well deflned state. The drawback 
is that many active systems, biological ones being partic- 
ularly pertinent examples, are far from equilibrium and 



one might for example miss some important physics by 
the restriction to the linear nonequilibrium regime. This 
phenomenological approach can, however, be very use- 
ful, especially when coupled to microscopic derivations 
for speciflc model systems. 

Regardless of the choice of hydrodynamic framework, 
the reasons for the choice of variables remain the same. 
When an extended system is disturbed from equilibrium 
by an external perturbation, its relaxation is controlled 
by the microscopic interactions among constituents. It 
is useful to divide the relaxation processes into fast and 
slow, and to build a theory of the slow dynamics in which 
the fast processes enter as noise and damping. Such a 
division is unambiguous when there are collective ex- 
citations with relaxation rates uj{q) that vanish as the 
wavevector q goes to zero; these are the hydrodynamic 



modes of the systems (Forster 1975 Martin et al. 1972). 



Familiar examples are diffusion and sound waves in fluids. 
A formulation in terms of hydrodynamic flelds provides 
a generic description of the nonequilibrium large scale 
physics that relies only on general properties and local 
thermodynamics, and does not depend strongly on mi- 
croscopic details of the interactions. In order to build a 
hydrodynamic theory, whether phenomenologically or by 
deriving it from a microscopic model, the flrst task is to 
identify the slow variables, which are the local densities 
of conserved quantities, the "broken-symmetry" variables 
which have no restoring force at zero wave number, and, 
in the vicinity of a continuous phase transition, the ampli- 



tude of the order parameter (Forster, 1975 Martin et al. 



1972). From there on, the procedure is well deflned and 
systematic. Conserved quantities are fairly easy to iden- 
tify. Momentum conservation for "wet" systems means 
that the momentum density is a slow, conserved variable, 
while for "dry" systems it is a fast variable. Identifying 
the type of order - polar or nematic, in the cases that 
we will examine- yields both the amplitude of the order 
parameter and the broken-symmetry modes. 

Several useful reviews on various aspects of the be- 
havior of flocks and active systems have appeared in 



the recent literature (Joanny and Frost 



2009 



Jiilicher 



et al.\ |2007t |Ramaswamyl 


2010 


Toner et al.\ 


2005 


Vic- 


sek and Zafeiris| 2012), as well as on the properties of 



swimmers (Ishikawa, 2009 Koch and Subramanian 2011 



Lauga and Powers, 2009) and bacterial suspensions and 



colonies (Ben- Jacob et a/., 2000 Gates 2012 Murray 



2003 ) . The present review aims at highlighting the unity 
of various approaches and at providing a classiflcation of 
active systems. It demonstrates the link between micro- 
scopic derivations and continuum models, showing that 
the hydrodynamic equations derived from different mi- 
croscopic models with the same general symmetry have 
the same structure and only differ in the details of param- 
eter values, and that the continuum models proposed and 
used in the literature can all be formulated in a unifled 
manner. It is hoped that this review will provide a use- 
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ful and self-contained starting point for new researchers 
entering the field. 

This review is organized as follows. We first consider 
dry active systems focusing on the ordering transition, 
the properties of the ordered phase and the differences 
between polar and nematic active matter. We also dis- 
cuss in section [ll| a system first analyzed theoretically 



by Baskaran and Marchetti (Baskaran and Marchetti 
2008a|b ) where the self-propelled particles are polar in 



their movement but nematic in their interaction and 
hence in the nature of their macroscopic order. In sec- 
tion [lllj we discuss orientable active particles suspended 
in a fluid, or active gels, for short. We give a system- 
atic derivation of the constitutive equations of a nematic 
active gel, for a system weakly out of equilibrium, and 
dwell briefly on the effects of polar order, viscoelastic- 
ity and the presence of multiple components. In section 
JVj we present some applications of the hydrodynamic 
theory of active matter to specific geometries that could 
arise in experiments. We focus on instabilities of thin 
films and on the rheological properties of active matter. 
Section |V| gives a brief account of microscopic theories 
of active matter. We offer examples to show how mi- 
croscopic theories allow one to determine in principle 
the phenomenological transport parameters introduced 
in the hydrodynamic theories and to estimate their or- 



der of magnitude. The concluding section |Vl| presents 
open questions related to hydrodynamic theories of ac- 
tive matter and gives some perspectives. 



II. "DRY" ACTIVE MATTER 

In this section we consider active systems with no mo- 
mentum conservation. This class includes bacteria glid- 



ing on a rigid surface (Wolgemuth et al., 2002), animal 
herds on land (Toner and Tu, [1^998) or, in the artifi- 



cial realm, vibrated granular particles on a plate (Aran- 



son and Tsimring, 2006 Deseigne et al\ 2010[ Kudrolli 



et al. 



2003 



2008] Narayan et al.\ 2007| Ramaswamy et al. 



Yamada et al. , 2003 ) in all of which momentum is 



damped by friction with the substrate. It is also plausi- 
ble that models without an explicit ambient fluid contain 
the main physics of concentrated collections of swimming 



bacteria (Drescher et al., 2011) and motor- filament sus- 



pensions (Liverpool 2003), where steric and stochastic 



effects could for most purposes overwhelm hydrodynamic 
interactions. Minimal flocking models have also been 
tested against observations on aerial displays by large 
groups of birds where, despite the much lower concentra- 
tion in comparison to bacterial suspensions, (presumably 
inertial) hydrodynamic effects seem negligible ( Bailer ini| 
eTaTI [20081 iGmelh and Chat^l [Mo] ). 

We refer to these systems as "dry" active systems. In 
this case the only conserved quantity is the number of 
particles (neglecting of course cell division and death) 



and the associated hydrodynamic field is the local density 
of active units. Over-damped active particles that can 
order in states with polar and nematic symmetries are 
described below. The first class consists of polar or self- 
propelled units with interactions that tend to promote 
polar order, i.e., explicitly align the particles head to head 



and tail to tail, as in the classic Vicsek model (Vicsek 



et al. , 1995). The second class consists of active particles 



that may themselves be apolar, such as melanocytes, the 
cells that distribute pigments in the skin, where activity 



induces non-directed motion on each cell ( Gruler et al. 
1999[) or polar, such as self-propelled hard rods (Baskaran 



and MarchettH |2008b| [Peruani et d] [20061 W^^Z et al. 
2010), but with interactions that tend to align particles 
regardless of their polarity, so that the ordered state, if 
present, has nematic symmetry. A cartoon of the differ- 
ent cases is shown in Fig. [S] 



A. Polar Active Systems: Toner and Tu continuum model 
of flocking 



A continuum effective theory for the flocking model 



of Vicsek et al. (1995) was proposed in 1995 by Toner 
and Tu| ( |1995[ |1998[ ) [see also (f Toner et a/^|2005[ )]. The 
Vicsek model describes a collection of self-propelled par- 
ticles with fixed speed and noisy polar aligning inter- 
actions, and displays a nonequilibrium phase transition 
from a disordered state at low density or high noise to 
an ordered, coherently moving state at high density or 
low noise strength. Toner and Tu formulated the con- 
tinuum model phenomenologically solely on the basis of 
symmetry considerations. Recently, the Toner- Tu model 
was derived by [Bertin et al.\ (|2006l |2009|) and [ihlel (|2011|) 



by coarse-graining the microscopic Vicsek model. These 
derivations provide a microscopic basis for the hydrody- 
namic theory, yielding explicit values for essentially all 
the parameters in the Toner-Tu model, except for the 
noise strength in the latter. The Boltzmann-equation 
approach of Bertin et al. ( |2006| ) and Jhle, ( |2011| ) leads 
to a deterministic coarse-grained description, with the 
stochasticity of the Vicsek model reflected in an average 
sense through the diffusion and relaxation terms. In this 
section we introduce the continuum equations in their 
simplest form and analyze their consequences. 

Since the particles are moving on a frictional substrate, 
the only conserved field is the number density p(r, t) of 
active particles. In addition, to describe the possibility 
of states with polar orientational order, one must con- 
sider the dynamics of a polarization vector field, p(r, t). 
These continuum fields can be defined in terms of the po- 
sition r^(t) of each active particle and a unit vector On{t) 
denoting the instantaneous orientation of the velocity of 
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polar active 
particle 



apolar active 
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self-propelled 
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polar order of 
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polar particles 



nematic order of 
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FIG. 8 (color online) Schematic of the various types of ac- 
tive particles and orient at ionally ordered states. Polar active 
particles (top left image), such as bacteria or birds, have a 
head and a tail and are generally self-propelled along their 
long axis. They can order in polar states (bottom left) or ne- 
matic states (bottom center) . The polar state is also a moving 
state with a nonzero mean velocity. Apolar active particles 
(top center image) are head-tail symmetric and can order in 
nematic states (bottom right). Self-propelled rods (top right 
image) are head-tail symmetric, but each rod is self-propelled 
in a specific direction along its long axis. The self-propulsion 
renders the particles polar, but for exclusively apolar inter- 
actions (such as steric effects), self-propelled rods can only 
order in nematic states (bottom center image). 



each particle as 



Kr,t) = ^<5(r-r„(t)), 



(la) 



Although we are dealing with a nonequilibrium system, 
it is convenient and instructive to write the dynamical 
equations in a form in which terms that can be viewed 
as arising from a free energy functional Fp are separated 
out: 



dtp = -V • {vopp) , 



(2a) 



dtp + Ai(p . V)p = — + f , (2b) 
7 op 



where vq is the self-propulsion speed of the active parti- 
cles, 7 a rotational viscosity, and Ai a parameter control- 
ling the strength of the advective term on the left hand 
side of Eq. (2b). The last term on the right hand side 



of Eq. ( 2b ) captures the fluctuations and is taken to be 



white, Gaussian noise, with zero mean and correlations 
{Ur,t)Mr\t'))=2A5^pS{r-v')S{t-t') . (3) 



For simplicity, and consisten cy with [Toner and Tu ( [1998 ) 



and with the derivation of Bertin et al p006fc of the 



Toner- Tu equations from the model of [Vicsek et aL\ 
(1995), we have neglected the diffusive current and the 
associated noise in the density equation in ([2|. In ([2| 
and everywhere in the following we have taken the noise 
terms to be purely additive, ignoring dependence on the 
local values of p (see, e.g., ( Dean[ 1996[ )) or of the or- 
der parameter (here p). This approximation is adequate 
for the present purposes of calculating two-point correla- 
tors in a linearized theory, but is actively under discus- 



sion ( 


Gowrishankar and Rao, 2012 Mishra, 


2009 


Mishra 


et al 




2012) for situations of strong inhomogeneity such 



as the coarsening of active nematics. 

Note that in the continuum flocking model p plays a 
dual role: on the one hand p is the orientational order 
parameter of the system, on the other vop represents the 
particle velocity field. This duality is crucial in deter- 
mining the large scale behavior of this nonequilibrium 
system. The free energy functional used in ([2| is given 

by 



P I |4 

^IpI 



K 



viV • p 



^ + -1 



2 



(4) 



where po is the average density and Sp = p — po. The 
first two terms on the right hand side of Eq. (|4| con- 
trol the mean-field continuous order-disorder transition 
that takes place as the parameter a goes through zero. 
In the derivation of Bertin et al. ([2006), a depends on 
local density p and the noise strength in the under- 
lying microscopic model, and turns negative at suffi- 
ciently large p. A reasonable phenomenological approach 
to describe the physics near the transition is to take 
a{p) = Dr (1 — p/pc)^ changing sign at a characteristic 
density pc- For stability reasons the coefficient /3 is posi- 
tive. The ratio Dr/^ has dimensions of frequency. In the 
following we will use it to set our unit of time by letting 
Dr/j = 1. We also divide implicitly all the coefficients 
in the free energy (|4| by Dr. The third term in the free 
energy describes the energy cost for a spatially inhomo- 
geneous deformation of the order parameter. The Frank 
constant K is positive. For simplicity we employ a one- 
elastic constant approximation and use the same value 
for the stiffness associated with splay and bend deforma- 
tions of the order parameter field, as also for gradients 
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in its magnitude. Splay /bend anisotropy can however, 
play a very important role in active systems, as demon- 
strated for instance by Voituriez et al (2006). The last 
two terms are allowed only in systems with polar symme- 
try [see, e.g., (Kung et al, 2006)]; they give the density 
and |pp contributions to the spontaneous splay. Upon 
integration by parts, they can be seen to enable gradients 
in p or in the magnitude |p| of the order parameter to 
provide a local aligning field for p. Alternatively, the A 
term can be viewed as a splay-dependent correction to 
a{p): splay of one sign enhances, and of the other sign 
reduces the local tendency to order. 

By using Eq. Q for the free energy, the equation for 
p takes the form 



1. Homogeneous steady states 



5tp+Ai (p • V)p = - [a{p) + I3\pf] p + KV^p 

-«iV^ + ^V|p|2-Ap(V-p) + f , 

Po 2 



(5) 



It is instructive to compare Eq. ([5| with the Navier- 
Stokes equation for a simple fluid. The term proportional 
to Al is the familiar advective nonlinearity. Unlike the 
fluid, the flocking model does not possess Galilean invari- 
ance as the particles are moving relative to a substrate. 
As a result, Xi ^ vq [which would have read Ai 7^ 1 had 
we not extracted a velocity scale in (2a)]. If we think of p 
as proportional to a flow velocity, the second term on the 
right hand side represents the effects of viscous forces. 
The third and fourth terms of Eq. ^ can be rewritten 
in an approximate pressure-gradient form as — (l/po)VP 
with a pressure P(p, |p|) vip — (Apo/2)|pp. This 
highlights the parallel and the contrast with the Navier- 
Stokes equation: in the latter, an equation of state de- 
termines the thermodynamic pressure in terms of density 
and temperature and not the velocity fleld. The flrst term 
on the right hand side of Eq. ([5| arises from the role of p 
as order parameter for a polarized state and does not have 
an analogue in equilibrium fluid flow. There is, however, 
an equilibrium analogue in the context of electrostriction, 
a property of all dielectric materials that causes them to 
change their shape under the application of an electric 
fleld, although in this case polarization does not corre- 
spond to velocity, but rather to reactive polarity. General 
symmetry arguments given by |Mishra et a/.| (2010) and 



Toner et al. 



( [2QQ5| allow for more general polar terms in 
Eq.([5|): ^V|p|^ + A2p(V • p) The derivation based on 
the free energy (|4|, using ([2b| yields A3 = — A2 = A. For 
the out-of-equilibrium system considered here the coeffi- 
cients of these two terms are not generally related, but 
the microscopic derivation of Bertin et al] ( |2009[ ) also 
yields A3 



The dynamical model described by Eqs. (2a) and (2b) 



exhibits by construction, in a mean-field treatment of ho- 
mogeneous configurations, a continuous transition from 
a disordered to an ordered state. For a > 0, correspond- 
ing to an equilibrium density po < Pc, the homogeneous 
steady state of the system is disordered or isotropic, with 
p = and a corresponding zero mean velocity. For a < 0, 
corresponding to po > pc^ the system orders in a state 
with uniform orientational order, with |po| = \/— c^oT^, 
where ao = Q^(po)- In the ordered state, which is also a 
moving state, with v = vqPo^ continuous rotational sym- 
metry is spontaneously broken. This mean-field analysis 
survives fiuctuation corrections even in two dimensions 



(Toner and Tu, 1995 1998), evading the Mermin- Wagner 
theorem forbidding - at equilibrium - the spontaneous 
breaking of a continuous symmetry in two spatial dimen- 



sions (Chaikin and Lubensky, 2000 Hohenberg, 1967 
Mermin and Wagner, 19661). In the present nonequilib- 
rium system, however, the theorem does not apply. 

The escape from Mermin- Wagner is primarily due to 
the advective nonlinearities in ([2| which effectively gener- 
ate long-range interactions in the system, even if density 



variations are ignored as shown recently by Toner (2012). 
The ordered phase itself, especially near the transition, is 
exceedingly complex as a result of coupling to the density. 
Although the mean density is an important control pa- 
rameter for the transition to a ffock, the local propensity 
to order depends, via a{p) in (|4|, on the local density 
p(r,t). This is analogous to the density-dependent ex- 
change coupling in a compressible magnet in Milosevic] 



let 



(1978). The analogy ends there: unlike in the mag- 



net, the order parameter feeds back strongly into the den- 
sity dynamics via the current on the right hand side of 
Eq. (|2]d). See section |II.A.3| of [Toner and Tu| ( |1998D , and 
the review by Toner et al. (2005) for a more complete 



discussion of this point. 



2. Properties of the isotropic state 

To study the linear stability of the homogeneous steady 
states we examine the dynamics of spatially inhomoge- 
neous fluctuations Sp = p — po and (^p = p — Po from the 
values po and po in each of the homogeneous states. The 
isotropic state is characterized by po < pc and po = 0. 
The equations for the dynamics of linear fluctuations 
about these values are 



dtp 



dtSp = -voPo V • p 



-aop 



Po 



(6a) 

VSp + KV'p + f , (6b) 



-A2 as well as A^ ^ Vq and vi = vo/2. 



where ao > 0. Fourier-transforming ([6| in space and 
time to look at modes of the form exp(iq • r — iujt) leads 
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to dispersion relations between frequency uj and wavevec- 
tor q for small fluctuations about the uniform isotropic 
phase. Polarization fluctuations transverse to q decou- 
ple and decay at a rate ao + Kq^ . The relaxation of 
coupled fluctuations of density and longitudinal polariza- 
tion is controlled by coupled hydrodynamic modes with 
frequencies 



(7) 

Linear stability requires that fluctuations decay at long 
times so that lm\uj((i)\ < 0. The isotropic state is lin- 
early stable for all parameter values, provided vqVi > 0. 
The parameter vi enters ([6| like an effective compres- 
sional modulus. Microscopic derivations of the con- 
tinuum model have yielded vi = vo/2 at low density 



(Baskaran and Marchetti 2008b Bertin et al , 2009 Ihle 



2011 ). In this case vi > and the isotropic state is always 



stable. At high density, however, caging effects can result 
in a density dependence of Vi , which can in turn lead to a 
density instability of the isotropic state as shown by |Fily| 
et al] (|2012|) andlTailleur and Catesl (|2008|). Here we re- 



strict our discussion to the situation where the isotropic 
state is linearly stable. Even in this case, it has unusual 
properties: approaching the mean- field ordering transi- 
tion, i.e., decreasing ao, is like decreasing friction. As ar- 
gued by|Ramaswamy and Mazenko ( 1982 ) for equilibrium 



fluids on a substrate, the dispersion relations can acquire 
a real part for ao < vqVi/K and a range of wavevectors, 
Qc- ^ Q ^ Qc-\-^ with (Baskaran and Marchetti 2008b| ) 



2voVi 



1 



Kap 
2voVi 



± A 1 



Kao 

VqVi 



(8) 



Density fluctuations then propagate as waves rather than 
diffuse. Near the transition, where 0+, we find 

Qc- and qc-\- — 2^v^V\lK and the propagating 
waves resemble sound waves, with cJj^ ^q^v^V\. These 



propagating sound-like density waves are ubiquitous in 
collections of self-propelled particles. We stress, however, 
that if the limit ^ is taken first, density fluctuations 
always decay diffusively. 



3. Properties of the ordered state 

We now examine the linear stability of the ordered 
state by considering the linear dynamics of fluctuations 
about po > 9c and po = pox. The detailed analysis be- 
low is for (i = 2, but the general conclusions including the 
instability of the uniform ordered state just past thresh- 
old hold in all dimensions. It is useful to write the order 
parameter in terms of its magnitude and direction, by 
letting p = pfi, where n is a unit vector pointing in the 
direction of local orientational order. Fluctuations in the 
polarization can then be written as ^p = rio^p + po^n. 
The condition |np = 1 requires that to linear order 
no • (5n = 0. For the chosen coordinate system no = x. 
Then in two dimensions ^n = (5ny and (5p = X(5p+ypo^^- 
The linearized equations are 

[dt + vopodx) Sp = -vopoV • Sp , (9a) 
{dt + XiPodx) Sp = -2\ao\ii5p + apoSp + A2Po(V • 5p) 
-{vi/po)Vdp + Aspo V(5p + KV^Sp + f , (9b) 

where a = — (a' + P'Po) — OiP' / P — a' > {) and primes 
denote derivatives with respect to density, e.g., a' = 



(Sp<^)p=PQ- Foi" the analysis below it is useful to dis- 
play explicitly the three coupled equations for the fluc- 
tuations. For compactness, we introduce a vector with 
elements = {^p, 5p^ 5n} and write them in a matrix 
form as 



a*(/.q(t)=M(q)-0q(t)+Fq(t) 

where the matrix M is given by 



(10) 



M(q) 



-nxVQPQ 

apo - iqxVi/po 
-iqyVi/{poPo) 



-^QxVoPo 
2|ao| - iqx>^Po 
iqy>^3Po 



Kq^ 



-iqyVoPoPo 

iQy>^2Po 

-iqx>^iPo - Kq^ 



(11) 



The sign of the combination A = Ai — A2 — A3 is im- 
portant in controlling the nature of the instabilities in 
the ordered phase. Notice that all microscopic deriva- 
tions of the continuous equations for various systems of 
self-propelled parti cles ( Baskaran and Marchetti[ [2008b 



Bertin et al, 2009) yield A > 0. Finally, the noise vec- 
tor Fq(t) has components Fq(t) = {0, /q (t), fq (t)}, with 
f^{t) = q-fq(t), f^(t) = fq(t)-q/^(t) andq = q/g. The 



dispersion relations of the hydrodynamic modes are the 
eigenvalues of the matrix M. The solution of the full cu- 
bic eigenvalue problem is not, however, very instructive. 



It is more useful to discuss a few simplified cases ( Bertin 



et al 2009 Mishra a/. 2010). 



a. Linear instability near the mean-field order-disorder tran- 
sition. We first consider the behavior of fluctuations as 
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the mean-field transition is approached from the ordered 
phase, i.e., for ao 0~ and po ^ . Fluctuations 6p 
in the magnitude of polarization decay at rate ~ |ao| 
and become long-lived near the transition. In fact in this 
limit one can neglect the coupling to director fluctua- 
tions and only consider the coupled dynamics of Sp and 
Sp. This decoupling becomes exact for wavectors q along 
the direction of broken symmetry, q = qit, The decay of 
director fluctuations is described by a stable, propagating 
mode, uJn = —Q^iPo — iKq^ . For small q the imaginary 
part of the frequencies oo^{c\) of the modes describing 
the dynamics of density and order parameter magnitude 
takes the form 



Im[ujl] = - [S2q^ + s^q^ + 0{q^ 



(12) 



with 52 



2\ao\ 



vod pQ 



4.vi\ao\^ 



, implying an instability at 



small q when S2 < 0, corresponding to |ao| < . Mi- 



croscopic calculations yield ao ~ and vi ^ vq (Bertin 



et al. 2009 ), indicating that the instability exists for arbi- 



trarily small Vq and in a narrow region of densities above 
the mean fleld order-disorder transition. The 0{q'^) term 



in Eq. \12\ stabilizes fluctuations at large wavevectors. 
In detail, near the mean-fleld transition (ao 0~) 
52 — — a^/8/3ao and 54 5a^/128/3^ao and the mode 
is unstable for q < q^ /3|aoP/ (5a^) ^ {pc — po)^^^- 

Numerical solution of the nonlinear hydrodynamic equa- 
tions with and without noise have shown that in this 
region the uniform polarized state is replaced by com- 
plex spatio-temporal structures. In the absence of noise, 
one flnds propagating solitary waves in the form of or- 
dered bands aligned transverse to the direction of broken 
symmetry and propagating along x amidst a disordered 



background ^Bertin et al 2009 Mishra et al 2010). As 



discussed below, these bands have been observed in nu- 



merical simulations of the Vicsek model (Chate et al. 



2008b) and are shown in Fig. 10 



like waves, with dispersion relation 

UJ± = qc± {9) - iq'JC± {9) + Oiq'), (14) 

with speeds and dampings c±{6) and /C±(6>), whose de- 
tailed forms are given in the review by Toner et al. (2005) 



and depend on the angle between q and the direction 
X of broken symmetry. For = 7r/2, i.e., q along we 
flnd c± = ±^yvoVl and JC± = {K-\-D)/2, while for 6 = 0, 
corresponding to q along x, the two modes are decou- 
pled and = qvo — iDq^ and = —q\i — iKq^. The 
difference in the propagation speeds is due to the lack of 
Galilean invariance. In a normal fluid, sound waves prop- 
agate because of mass and momentum conservation and 
inertia. In contrast, the presence of propagating long- 
wavelength density disturbances here is a signature of 
the spontaneously broken orientational symmetry which 



renders 5n "massless" (Ramaswamy, 2010). 



Going back to the full form of the matrix M in (11), 



we notice that to leading order in q the relaxation of 5p 
is governed by the equation 



q,u; 



-2|ao|(5pq,c^ + apoSp^^^ + 0{q) . (15) 



This reminds us to be precise about what it means for 
a mode to be "fast": fluctuations in the magnitude of 
the order parameter cannot be neglected, but rather 
are slaved to density fluctuations on long timescales: 
Spci^uj — (<^Po/2|ao|)^Pq,u;. Inserting this in (10) yields 
an effective dynamics for 5n and p with an instability 



of splay fluctuations if A3 > 



> (Mishra et al. 



2010). An analysis of the eigenvalues of the full matrix 



M conflrms this result. 



c. Giant density fluctuations. The linearized equations 
with noise can also be used to calculate correlation func- 
tions. Of particular interest is the static structure factor 



b. "Sound waves". Deep in the ordered phase, i.e., for 
large negative ao, fluctuations dp in the magnitude of the 
order parameter decay on non- hydrodynamic time scales. 
We assume Sp has relaxed to zero on the time scales 
of interest and let for simplicity po — 1- In addition, 
we include a diffusive current and the associated noise 



that were neglected in the density Eq. (2a). The coupled 



equations for density and polar director fluctuations, Sp 
and Sn, take the form 

{dt + vod,) Sp = -voPodySn + DV^Sp + V • f ^ , (13a) 



{dt ^ Xidx)Sn = -vidy 



KV^Sn + fy , (13b) 



where I) is a diffusion constant and fp describes white, 
Gaussian noise. The hydrodynamic modes, obtained 



from Eqs. (13a) and (13b), are stable propagating sound- 



where 



2^ 



5(q,c^), (16) 



1 

^(q,^) = ^ / dteM^ujt){Sp^{0)Sp_^{t)) (17) 

PQ^ Jo 

is the dynamic structure factor. We are interested in cal- 
culating density fluctuations in the region away from the 
mean-fleld transition, where the ordered state is linearly 



stable. Eqs. (13a) and (13b) yield 6'(q, a;). For a reader 
wishing to work through the calculation, note that to 
leading order in the wavevector the noise fp in the den- 
sity equations can be neglected, except for q along the 
direction of broken symmetry {0 = 0). Integrating over 



frequency as in Eq. (16), again to leading order in q, one 
obtains 
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^(q) 



(I'o + Ai)^ cos^ 6> + sin 6 



K+(0) K-{0)\ q 



(18) 



For 6> = the leading singularity displayed in (18) 
vanishes; to determine the finite value of 5'(q) in that 
limit requires going to higher order in q which we will not 
do here. As we will see below, the divergence of the static 
structure function at large wavelengths (as in the 
present linearized treatment) is a remarkable and robust 
property of uniaxially ordered active systems, whether 
polar or nematic, and follows generally from the orienta- 
tional order that develops in a sufficiently dense collection 



of self-driven particles with anisotropic body shape (Ra- 



Toner and Tu 



1995 



maswamy et a/.[ [2QQ3| |Simha and Ramaswan^ |2QQ2a|b| 



1998 Toner et al 2005). The l/q^ is 



also a dimension- independent property of our linearized 
theory; the numerical prefactors do, however, depend on 
dimensionality d. The divergence of the static structure 
function for g ^ implies an important violation of the 
familiar scaling of number ffuctuations in equilibrium. In 
the limit of vanishing wavevector the structure factor is 
simply a measure of number ffuctuations, with 



(19) 



where N and (A^) are the instantaneous and average 
number of particles in a region of size V, respectively, 
and AA/"^ = ((A" — (A"))^) is the variance of number 
ffuctuations. In an equilibrium system AA" ~ y^^C^^ 
so that AN/{N) - 1/vW ^ as {N) oo. The 
~ 1/q^ divergence of S{q) for g ^ predicted by 
Eq. ( p!8| ) implies that in an active system AA" grows faster 
than ^^/{N). To understand this we rewrite Eq. (18) as 
AA^ ~ Y^^C^^O^^T^~^~0)- Assuming that the smallest ac- 
cessible wavevector is of the order of the inverse of the 
system size V~^/^, with d the space dimension, we ffnd 
S{q ^ 0) - V^/^ - (A')2/^. This gives 



AN - (A/-)^ , 



1 

d 



(20) 



The linear theory reviewed above predicts a strong en- 
hancement of density ffuctuations in the ordered state, 
with an exponent a = 1 in two dimensions (Toner and 
Tu, 1995) [Toner et aL\ [2005] ). We wih see in section 
II. B| that a similar linear theory also predicts the same 
scaling in the ordered state of active apolar nematic. In 
general, we expect the exponent to be substantially mod- 
iffed as compared to the theoretical value by ffuctuation 
effects arising from the coupling of modes at different 
wavenumbers as a result of nonlinear terms in the equa- 
tions of motion. "Giant density ffuctuations" have now 



been seen experimentally in both polar ( Deseigne et al. 
2010) and apolar (Narayan et al. , 2007| vibrated granular 




FIG. 9 (color online) Swimming bacillus subtilis bacteria ex- 
hibit strong polar order and fluctuations that that are smaller 
at low density (A) than at high density (B). Nearby bacte- 
ria with arrows of the same color belong to the same dynamic 
cluster; the arrows indicate a b acterium's speed and direction. 
Adapted with permission from Zhang et al. (2010). 



of the apolar case. The scaling exponents, although mea- 
sured over a limited dynamical range, are consistent, for 
the apolar case, with the linearized treatment presented 
here and in (Ramaswamy et al. [2003 ), and with the 



renormalizat ion- group treatment (Toner and Tu 1995 



1998) [Toner et al] |2005|) in = 2 for polar ffocks. Giant 



number ffuctuations have also been seen in simulations 



of Vicsek-type models by Chate and collaborators ([Chate 
2QQ8a[^ 



et a/.[|2008a|b[ [2006| ). Much remams to be understood 
about number fluctuations in flocks, as revealed for in- 



stance by the experiments of (Zhang et a/., 2010) (see 
Fig.|9|. Recent work has also revealed other mechanisms 
that can yield large number fluctuations and even true 
phase separation in active systems in the absence of an 



ordered state (^Cates) [2012[ [Fily et al.\ [2012| jTailleur and 



Gates, 2008). We refer the reader to section II. G for fur- 



matter. See section Hi. B. II for a more detailed discussion 



ther discussion of this point. 

Recall that in thermal equilibrium the thermodynamic 
sum rule relates S{q = 0) to the isothermal compressibil- 
ity XT = - y (l^) according to S{q ^ 0) = p^ksTxT, 

with po = {N) /y the mean number density. One might 
then be tempted to say that orient at ionally ordered ac- 
tive systems are characterized by a diverging effective 
compressibility. This would in general be wrong. Aug- 
menting the free-energy functional (|4| with a linear cou- 
pling J U{r,t)5p{r,t)d^r and calculating the response 
function {S (p) / dU) ^^^j can readily be shown to yield a 
finite result for q, ^ 0. The giant density ffuctuations 
seen in the ordered state of active systems are an ex- 
cess noise from the invasion of the density dynamics by 
the soft Goldstone mode of orient ational order, not an 
enhanced response to perturbations. 
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d. Long-range order of dry polar active matter in 2d. The 
ordered polar state is remarkable in that it exhibits long- 
range-order of a continuous order parameter in two di- 
mensions. This has been referred to as a violation of 
the Mermin- Wagner theorem, although it is important 
to keep in mind that the latter only holds for systems in 



thermal equilibrium. From Eqs. (13a) and (13b) one can 
immediately obtain the correlation function of fluctua- 
tions in the direction Sn of polar order, with the result 
(|(5nq(t)p) ~ 1/^^. Taken literally, this result implies 
quasi-long-range order in two dimension in analogy with 
the equilibrium XY model. It has been shown, however. 



that nonlinearities in Eqs. (13a) and (13b) are strongly 
relevant in d = 2 and lead to a singular renormalization 
of the effective stiffness of the polar director at small 
wavenumbers so that (|Jnq(t)p) diverges more slowly 
than for most directions of q, thus preserving long- 
range order ( |Toner and Tu[ |1995| |1998[ [Toner al 
2QQ5|). A qualitative explanation of how this happens 



can be given as in Ramaswamy (2010), but a remarkable. 



pictorial yet quantitative account can be found in lecture 



notes by Toner (2009). Exact arguments for scaling ex- 



ponents in two dimensions are given by Toner (2012) in 
the limit where particle number is fast, i.e., not conserved 
locally. 



e. Numerical simulation of Vicsek-type models. Although 
a review of numerical simulations of agent based models 
is beyond the scope of this paper, for completeness we 
summarize the main findings. Quantitative agreement 
has been found between numerical experiments on micro- 
scopic models and the predictions of the coarse-grained 
theory, including long-range order in d = 2, the form 
of the propagating modes, anomalous density fluctua- 
tions, and super-diffusion of tagged particles (|Gregoire| 
|2004[ [Toner et al] |2005[ ). Early Immeri- 



and Chate 



et 



1G24 



512 




1024 



FIG. 10 (color online) Typical snapshots in the ordered 
phase. Points represent the position of individual particles 
and the red arrow points along the global direction of motion. 
Note that in the simulations sharp bands can be observed only 
if the system size is larger than the typical w i dth of t he bands. 
Adapted with permission from Chate et al. (2008a). 



continuum model describe d by Eq. ([2|) has also yielded 



traveling wave structures (Bertin et al. 2009 Gopinath 



et al.^ 2012 Mishra et al. 



2010). Finally, traveling den- 



sity waves have been observed in actin motility assays 
at very high actin density ( Butt et al.\ 2010[ Schaller 



et a/., 2010), as shown in Fig. 11 Although it is tempt- 
ing to identify these actin density waves with the trav- 
eling bands seen in simulations of the Vicsek model, the 
connection is at best qualitative at the moment. In par- 
ticular, it remains to be established whether these actin 
suspensions can indeed be modeled as dry systems, or 
hydrodynamic interactions are important, as suggested 



in Schaller et al. (2011a). 



cal studies (Vicsek et a/., 1995) and some later variants B. Active Nematic on a substrate 



(Aldana et al.^ 2007 Baglietto and Albano, 2009 Gonci 



2008 ) found a behavior consistent with a continu- 
ous onset of flocking, as a mean-field solution would pre- 
dict. However, the currently agreed picture, in systems 
where particle number is conserved, is of a discontinu- 



ous transition ( 


Chate et al. 


2007, 2008b 2006 


Gregoire 


and Chate, 2004[). It was i 


bund in fact that in a large 



domain of parameter space, including the transition re- 
gion, the dynamics is dominated by propagating solitary 
structures consisting of well-defined regions of high den- 
sity and high polar order in a low-density, disordered 
background. The ordered regions consist of stripes or 
bands aligned transverse to the direction of mean order 
and traveling along the direction of mean motion with 
speed of order vq^ as shown in Fig. [To] Away from the 
transitions, for weak noise strengths, a homogeneous or- 
dered phase is found, although with anomalously large 
density fluctuations. Numerical solution of the nonlinear 



Nematic order is the simplest kind of orientational or- 
der, with a spontaneously chosen axis which we shall des- 
ignate the X direction, and no distinction between x and 
—X. Lacking a permanent vectorial asymmetry, the sys- 
tem, although driven, has no opportunity to translate its 
lack of time-reversal symmetry into an average nonzero 
drift velocity. The reader would be forgiven for think- 
ing that such a non-moving state cannot possibly dis- 
play any of the characteristics of a flock. Remarkably, 
however, some of the most extreme fluctuation proper- 
ties of ordered phases of active particles are predicted 



(Ramaswamy et al. 2003) and observed in experiments 



( Narayan et al. , 2007 ) and earlier in simulations ( Chate 



et al. , 2006 Mishra and Ramaswamy , 2006) on this flock 
that goes nowhere on average. As with polar active sys- 
tems, several levels of description are possible for an 
active nematic - minimal agent-based stochastic mod- 
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FIG. 11 Density waves in a dense actin motility assay. The 
left image is a snapshot of the traveling wavefronts of actin 
observed in motility assays for actin density abo ve about 
20 filament s//im^, adapted with permission from Schaller 
et al. (2010). The right image (courtesy of V. Schaller) is ob- 
tained by overlaying maximal intensity projections of 20 con- 
secutive images, corresponding to a total time of 2.5 seconds. 
The time overlay allows one to trace the trajectory of the 
filaments, showing that filaments in the high-density regions 
move collectively with high orient at ional persistence, while 
filaments lying outside the bands perform uncorrelated persis- 
tent random walks. The elongated density modulation travels 
in the direction of the filaments' long direction at a speed that 
remains approximately constant over the time scale of the ex- 
periment. 



els (Chate et al. 2006) and Fokker-Planck treatments 



thereof (Shi and Ma, |2Q1Q), more detailed dynamical 



models that generalize rod-like-polymer dynamics to in- 
clude the physics of motors walking on cytoskeletal fil- 
aments (Ahmadi et a/., 2005, [2006 j , or direct contin- 
uum approaches based on partial differential equations 



(PDEs) for the slow variables (Simha and Ramaswamy 



2002a). The microscopic models can in turn be coarse- 



grained ([Ahmadi et al] |2006[ |Mishra[ |2009[ [Mishra and 
Ramaswamy I 2011 Shi and Ma, 2010) to yield the con- 



tinuum equations in averaged (Ahmadi et al. 2005; Shi 



and Ma, 2010) or stochastic (Mishra, 2009 Mishra and 



Ramaswamy, 2011) form. 



1. Active Nematic 

In this section we construct the equations of motion 
for active apolar nematics and discover their remarkable 
properties. We continue to work in a simplified descrip- 
tion in which the medium through or over which the ac- 
tive particles move is merely a momentum sink, without 
its own dynamics. Effects associated with fluid flow are 
deferred to a later section. We consider again a collection 
of elongated active particles labeled by n, with positions 
r^(t) and orientation described by unit vectors i>n{t). As 
in the case of polar systems on a substrate, the only con- 
servation law is that of number. The slow variables of 
such a formulation are thus the number density p(r, t) at 
location r and time t, as defined in Eq. (la) in section 



eter, a traceless symmetric tensor Q with components 



?a/3(r,t) 



(21) 

measuring the local degree of mutual alignment of the 
axes of the constituent particles, without distinguishing 
head from tail. In constructing the equations of motion 
we keep track of which terms are permitted in a sys- 
tem at thermal equilibrium and which arise strictly from 
nonequilibrium activity. As we are ignoring inertia and 
fluid flow, the dynamics of the orientation is governed by 
the equation 



IqdtQ 



SQ 



(22) 



describing the balance between frictional torques gov- 
erned by a rotational viscosity 7q and thermodynamic 
torques from a free-energy functional 



Fc 



Q Q 



+CqQ : VV 



Po 



4 
A 
2 



(Q : Q)' 

Sp ' 
Po 



Kq 
2 



(VQ)^ 



(23) 



that includes a tendency to nematic order for aq < 
and /3q > 0, Frank elasticity as well as variations in the 
magnitude of Q via Kq, bilinear couplings of Q and p, 
and a compression modulus A penalizing fluctuations Sp 
in the density about its mean value po- Note that all pa- 



rameters in (23) can depend on p, as indicated explicitly 
for aq. We also stress that this form of the free energy 
only applies in two dimensions. For active nematic in 
three dimensions terms of order also need to be in- 
cluded. We have included in (22) a spatio-temporally 



white statistically isotropic tensor noise fg to take into 
account fluctuations of thermal or active origin. It is im- 



plicit that all terms in (22) are traceless and symmetric. 
We have ignored several possible terms in the equation of 
motion for Q - additional couplings to p via off-diagonal 
kinetic coefficients, nonlinear terms not expressible as the 
variational derivative of a scalar functional, and terms of 
sub-leading order in gradients - none of which plays an 
important role in our analysis. The unique physics of the 



self-driven nematic state (Ramaswamy et al.\ [2003[ ) en- 
ters through active contributions to the current J in the 
conservation law 



dtp ■ 



V J. 



(24) 



There are several ways of generating these active terms: 
one can write them down on the grounds that nothing 



rules them out in a nonequilibrium state 


( Ramaswamy 


et al. 


2003 


Simha and Ramaswamy 


2002a 


, derive them 



from specific microscopic models involving motors and 



II. A and the apolar nematic orientational order param- filaments (Ahmadi et al., 2005,[2006j) or apolar flocking 
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FIG. 12 (color online) Curvature generates temporary polar- 
ity in a nematic. Activity gives rise to a current in a direction 
determined by this polarity. 




FIG. 13 (color online) Curvature gives rise to self-propulsion 
in an apolar active medium, illustrating Eq. (25). A discli- 



nation of strength +1/2 (circled in red) in a nematic phase 
of a vibrated layer of rods. The locally polar director struc- 
ture, together with the absence of time-reversal symmetry in a 
driven state leads to directed motion. The curvatures around 
a nearby defect of strength -1/2 (circled in green) lack a pre- 
ferred direction; that defect does not move. Adapted with 



permission from Narayan et al. (2007) 



agents (Chate et al. 2006 


Mishra 2009 


Mishra and Ra- 


maswamy , 


2011), or frame them in the language of lin- 



ear irreversible thermodynamics (de Groot and Mazur 



1984| iJiilicher et aL\\2007l . This last we will do later, in 
section III.B.3[ For now. Figure [12] provides a pictorial 
argument. In a system with only an apolar order param- 
eter, curvature V • Q implies a local polarity and hence 
in a nonequilibrium situation maintained by a nonzero 
activity must in general give a current 



' active 



CqV • Q , 



(25) 



where (q is a phenomenological active parameter. 
Graphic experimental evidence of curvature- induced cur- 



rent can be found in Fig. [131 for details see [Narayan et ( 



(2007), Supporting Material. As shown by Ramaswamy 



et al. (2003) and reiterated below, this (q term leads to 
giant number fluctuations in active nematics. This is de- 



spite the fact that (25) has one more gradient than the 
current in the polar case in ([2|. We will comment below 
on this subtlety after we summarize the calculation and 
the relevant experiments and simulations. 



Using (23) in (22), and including in the current J in 



( 24 ) both the active contribution ( 25 ) and a purely pas- 
sive diffusive part — ^ , with a mobility 7"^, leads 
to 



dtp ■ 



DV^p^BV^W : Q + CqVV : Q + V-F , (26) 



with D = A/(pq7), B = Cq/(po7)- We have treated (q 
as constant, and we have introduced number-conserving 
fluctuations through the random current f ^ which, again, 
we take to be statistically isotropic and delta-correlated 
in space and time. 

Despite its apparent structural simplicity the apolar 
active nematic state has been explored less than its po- 
lar counterpart. The ordering transition appears numeri- 



cally ( Chate et al. 2006 ) to be of the Kosterlitz-Thouless 
type, but analytical approaches must confront the diffi- 
culties of the coupling of order parameter and density 
and have not been attempted so far, beyond mean- field 



(Ahmadi et al. 2005, 2006) which predicts a continuous 
transition in two dimensions. Working with equations 



of a form similar to ([22| and ([26|), |Shi and Ma| ( |2010[ ) 
obtained in a Boltzmann-equation approach from the 
model of Chate et al.\ ( |2006| ,with a finite- wavenumber 
instability of the uniform nematic state just past on- 
set, and a "perpetually evolving state". Note that simi- 
lar behavior has been obtained by Giomi and collabora- 
tors for "wet" nematic fluids (IGiomi et a/.l 120111 120121) 



and has been observed in very recent experiments on in 



vitro suspensions of microtubule-kinesin bundles (Dogic 



2012). Finally, there are preliminary numerical indica- 



tions ( Mishral 2009 Mishra and Ramaswamy, 2011) of 



anomalous growth kinetics of active nematic order fol- 
lowing a quench from the isotropic phase, with a strong 
clumping of the density. We will not pursue these is- 
sues further in this review. In the remainder of this sec- 
tion we restrict ourselves to a simple understanding of 
the statistics of linearized density fluctuations deep in a 
well-ordered active nematic. Tests of the resulting pre- 
dictions, through experiments and numerical studies, will 
be discussed thereafter. 

For the purposes of this section it is enough to work 
in two dimensions denoted (x,?/), so that Q has com- 
ponents Qxx = -Qyy = Q = SCOS20, Q^y = Qy^ = 

> in 



cyy 

P = Ssm20. We take aq < 0, /3q 
yielding a phase with uniform nematic order and den- 
sity po in mean-field theory, and choose our basis so 
that this reference state has ^ = 0. Expressing (22) and 
( 26 ) to linear order in perturbations about the reference 
state, p = po -\- Sp{r^t)^ Q = diag(6', — 5*) + (5Q, where 
^Qxy = ^Qyx — 2S0 and 6Qxx = —^Qyy = 0(0'^)^ and 
working in terms of the space-time Fourier transforms 



(23), 
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it is straightforward to calculate the correla- confirmed first in a computer simulation of an apolar 



tion functions 5'^,., = 
and 

details of this calculation (iRamaswamy et a/.l 120031), we 



h Pciuj^aoj) ' Rather than dwelling on the 



flocking model by Chate et al. (2006) and then in an 
experiment on a vertically vibrated monolayer of head- 
tail symmetric millimeter-sized bits of copper wire by 



suggest that readers verify for themselves a few essen- 
tial facts, (a) Parameter ranges of non-zero measure can 
be found such that the real parts of the eigenvalues of 
the dynamical matrix for this coupled problem are neg- 
ative, i.e., active nematics are not generically unstable, 
(b) Terms involving Cq from (23) are irrelevant - sub- 



Narayan et al. (2007). Of particular interest is the fact 



that the experiment observed a logarithmic time decay 
of the autocorrelation of the local density, a stronger test 
of the theory than the mere occurrence of giant fluctu- 
ations. Issues regarding the origin of the fluctuations 



dominant in wavenumber in (26) and leading only to a 
finite renormalization of parameters in ^f2\. (c) Ignor- 



in the experiment were nonetheless raised by |Aranson 



et 



(2008), with a response by Narayan et al. (2008). 



Faced with such large density fluctuations, it is natural to 
ask whether the active nematic state is phase-separated 



in the sense of Das and Barma (2000). This idea was 



tested by Mishra and Ramaswamy (2006) to a Simula- 



ing Cq, the dynamics of Q through (22) becomes au- 
tonomous, independent of p. In the nematic phase it is 
straightforward to show that has equal-time correla- 
tions 5'q = /((ia;/27r)5'q^ ex l/q^, and a lifetime ex l/q^. 

The term VV : Q through which 9^^ appears in the matic curvature ( [Lebwohl and Lasher[|1973D , a nd indeed 



tion model in which particles were advected by the ne- 



density equation (26) can be re-expressed as dxdyO. This 
autonomy of means that this term can be viewed as 
colored noise with variance cx q'^qy/q'^, and relaxation 
time ~ 1/^^. The zero- frequency weight of this noise 
is thus ^ qxqy/q^^ i.e., of order q^ but anisotropic. By 
contrast, the conserving noise V • has zero-frequency 
weight that vanishes as q^ and is thus irrelevant in the 
face of the active contribution VV : Q, whose long-time, 
long-distance effects are now revealed to be analogous 
to those of a number-non-conserving noise. The deter- 
ministic part of the dynamics of p, on the other hand, 
is diffusive and hence number-conserving. This combi- 
nation implies immediately that the equal-time corre- 
lator ^ noise strength x relaxation time (^x^l^ ^ 
i.e., direction dependent but of order l/g^, as argued in 



greater detail in section II.A.3.C for a linearized descrip- 
tion of density fluctuations in a ^olar flock. 

There is a puzzle here at first sight. We remarked 
above that the apolar contribution to the active current 



( 25 ) is one order higher in gradients than the polar term. 



Why then do density fluctuations in the linearized the- 
ory for polar and apolar flocks show equally singular be- 
haviour? The reason is that in polar flocks density gradi- 
ents act back substantially on the polarization p through 
the pressure-like term ^'iVp/po in ([5|, reducing the ori- 
entational fluctuations that engendered them. The effect 
of density gradients on d in the apolar case is far weaker, 
a quadrupolar aligning torque ~ d^dzP which does not 
substantially iron out the orientational curvature that 
engendered the density fluctuations. 

It should be possible to check, on a living system, the 
rather startling prediction that flocks that go nowhere 
have enormous density fluctuations. Melanocyte suspen- 
sions at high concentration show well-ordered nematic 
phases when spread on a substrate as seen by [Gruler 
et al. (1999). However, we know of no attempts to mea- 



sure density fluctuations in these systems. The predic- 
tion of giant number fluctuations has been tested and 



the steady state showed the characteristics of Das and 



Barma (2000)'s fluctuation-dominated phase separation. 



2. Self-propelled hard rods: a system of "mixed" symmetry? 

The Vicsek model contains an explicit alignment rule 
that yields an ordered polar, moving state. The align- 
ment interaction is in this case explicitly polar, in the 
sense that it aligns particles head-to-head and tail-to- 
tail. A different model of self-propelled particles where 
the alignment interaction has nematic, rather than polar, 
symmetry has been studied in detail by [Baskaran and 
Marchetti| ( |2008a|b| [2QTq| . These authors considered the 
dynamics of a collection of self-propelled hard rods mov- 
ing on a frictional substrate and interacting via hard core 
collisions. Numerical simulations of self-propelled hard 



rods have also been carried out by Ginelli et al. (2010); 



Kraikivski et al. (2006); Peruani et al. (2006 2008fc; and 



Yang et al. (2010) and have revealed a rich behavior. The 
model consists of rods interacting via excluded volume in- 
teractions and self-propelled at speed along their long 
axis. The detailed analysis by [Baskaran and Marchetti 



(2010) of the collisions of two self-propelled hard rods 



(see Fig. 14) shows that although self-propulsion mod- 
ifies the collision by enhancing longitudinal momentum 
transfer and effectively tends to align two colliding rods, 
the alignment is in this case "apolar", in the sense that 
it tends to align particles without distinguishing head 
from tail. .Baskaran and Marchetti, (,2010) used tools of 



nonequilibrium statistical mechanics to derive a modified 
Onsager-type theory of self-propelled hard rods, obtain- 
ing a Smoluchowski equation that is modified in several 
ways by self-propulsion. They then obtained hydrody- 
namic equations by explicit coarse-graining of the kinetic 
theory. One result of this work is that self-propelled hard 
rods do not order in a polar moving state in bulk. In spite 
of the polarity of the individual particles, provided by 
the self-propulsion, the symmetry of the system remains 
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FIG. 14 (color online) The cartoon shows that the binary col- 
lision of self-propelled hard rods provides an apolar aligning 
interaction, in the sense that it aligns rods regardless of their 
polarity. 



nematic at large scales. Self-propulsion does, however, 
enhance nematic order, as shown earlier by |Kraikivs"E| 



et al. (2006) in a model of actin motility assays, where 



actin filaments are modeled as rigid rods moving friction- 
ally on a in two-dimensional substrate and the propulsion 
is provided by myosin motors tethered to the plane. It 
is well known (in a mean field approximation neglecting 
fluctuations) that thermal hard rods of length ^ undergo a 
continuous isotropic-nematic t ransition in two dimensio n 
at a density pons = 3/(27r-^^) (Doi and Edwards, 1986), 
as first shown by Onsager in 1949 (Onsager, 1949). Self- 



propelled rods perform a persistent random walk (ten 



Hagen et aL\ |2QQ9l [Kareiva and Shigesada[ |1983| |Ku- 



drolli[ 2010) consisting of ballistic flights randomized by 



rotational diffusion. At long time the dynamics is diffu- 
sive and isotropic, with effective diffusion constant = 
D + VQ/2Dr^ where D is the longitudinal thermal diffu- 
sion coefficient and the rotational diffusion rate. Self- 
propelled rods of length i then behave as rods with an ef- 
fective length ieff ^ V^e/Dr = l^l^vl/{2{DDrY/'^ 
and undergo an isotropic-nematic transition at a density 

2 

s/(l + n2 )? where we have used D ~ i'^Dr. 



Pin = POn 



The microscopic derivation carried out in|Baskaran and 



Marchetti" f2008a) supports this estimate. This result 
has been verified recently via large scale simulations of 



2010). 



self-propelled hard rods (Ginelli et a/., 2010 Yang et al. 



The statistical mechanics of self-propelled hard rods 
and the derivation of the modified Smoluchowski equa- 
tions is described in detail in IBaskaran and Marchettil 



(2010) and will not be repeated here. Instead the hydro- 



dynamic equations for this system are introduced phe- 
nomenologically below and their consequences are dis- 
cussed. 



Once again, the only conserved field is the density 
p(r,t) of self-propelled rods. Theory and simulations 
have indicated that although self-propelled particles with 
apolar interactions do not order in polar states, their 
dynamics is characterized by an important interplay of 
polar order, described by the vector order parameter p 



defined in Eq. (lb), and nematic order, characterized by 
the alignment tensor Qij of Eq. (21). In other words. 



although the polarization p does not describe a broken 
symmetry, it seems necessary to incorporate its dynamics 
to capture the rich physics of the system into a contin- 
uum description. We then consider coupled equations for 
the conserved particle density, p, and the two coupled ori- 
entational order parameters, p and Q. Once again, the 
equations can be written in terms of a "free energy" as 

dtp = -V • {vopp) + DV^p + - VV:-^ + V • F , 

(27a) 

dtp + Al (p • V)p - ^ip • Q + ^2Q : Qp = - 1 ^ + f , 

(27b) 

+ f ^ , (27c) 



5,q + a;(p- v)Q 



1 

7Q 



6R 



pQ 



SQ 



ST 



where the subscript ST denotes the symmetric, trace- 
less part of any second order tensor T, i.e., [TJ^t has 
components Tfj^ = \ {Tij + Tji) - \6ijTkk, and 



K 



^P{p,P,S) (V • p) - V2Qa^{daP^)] , (28) 

We have grouped several terms (including some that were 
made explicit for instance in Eq.(|4|) by writing a local 
spontaneous splay in terms of a "pressure" P{p^p,S). 
The reader should note that it is essential to retain the 
general dependence of this pressure on p, p and S to 
generate all terms obtained from the microscopic theory. 
Here p = |p| and the magnitude S of the alignment tensor 
has been defined by noting that in uniaxial systems in two 
dimensions we can write Q^/s = 25 (nQ,n/3 — ^(^a/?), with 
n a unit vector. Then Qa-fQ-f(3 — S'^^ajS- The kinetic 
theory of Baskaran and Marchetti yields a^g > at all 
densities, indicating that no isotropic-polar transition oc- 
curs in the system. In contrast, Q^q(p) is found to change 
sign at a characteristic density p/at, signaling the onset 
of nematic order. The closure of the moments equations 
used in Baskaran and Marchetti ( 2008a|b ) to derive hy- 
drodynamics only gives terms up to quadratic in the fields 
in the continuum equations, but is sufficient to establish 
the absence of a polar state and to evaluate the renormal- 
ization of the isotropic-nematic transition density due to 
self-propulsion. Cubic terms, such as the ~ /3qQ^ needed 
to evaluate the value of the nematic order parameter in 
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the isotropic state is stable. It can also support finite- 



PlN 
Pon 




FIG. 15 (color online) The mean- field isotropic- nematic tran- 
sition line for self-propelled hard rods. Ada pted with permis- 
sion from Baskaran and Marchetti (2008a). 



wavevector propagating sound-like waves (Baskaran and 
Marchetti, 2008b). The properties of the ordered state, 



on the other hand, are more subtle, and not yet fully ex- 
plored. If we consider a homogeneous perturbation of the 
nematic state, we find that fluctuations Spx of the polar- 
ization along the direction of nematic order are unstable 
if X[So > o^pQ- This global instability occurs outside the 
hydrodynamic regime as it involves fluctuations in the 
polarization and its significance remains a bit of a puz- 
zle. Deep in the nematic state, neglecting fluctuations 
of both polarization and magnitude of the nematic order 
parameter, the homogeneous nematic state is unstable 
above a critical self propulsion speed. The instability 
arises from a subtle interplay of splay and bend defor- 
mations and diffusion longitudinal and transverse to the 
direction of nematic order and has been discussed in de- 



tail in Baskaran and Marchetti (2008b). 



the ordered state and the active ~ * Qp can be ob- 
tained from a higher order closure and are included here 
for completeness. A higher order closure also confirms 
the absence of a term proportional to in the polar- 
ization equation, which could, if present, yield a bulk 
polar state. We note that a p • Q • p term in the free 
energy is permitted and would yield both the p • Q term 
on the left hand side of the p equation, as well as a pp 
term in the Q equation, with related coefficients. The 
joint presence of these terms can lead to the existence 
of a polar ordered state. The hard-rod kinetic theory 
of Baskaran and Marchetti produces only the p • Q term. 
The fact that the theory does not generate the p • Q term, 
which would be obligatory had the dynamics been gov- 
erned entirely by a free-energy functional, underlines the 
nonequilibrium nature of the dynamics we are construct- 
ing. It would be of great interest to find the minimal 
extension of self-propelled hard-rod kinetic theory that 
could produce a phase with polar order. 



a. Homogeneous Steady States and their Stability. The 
only homogeneous steady states of the system are an 
isotropic state with p = pQ and p = Q = and a ne- 
matic state where p = 0, but Q is finite. The transition 
occurs at pjN where apQ ~ {piN — p) changes sign. As 
in the case of the active nematic, we work in two di- 
mensions denoted x and so that Q has components 
Qxx -Qyy Scos20, Qxy Qyx Ssin20. In the 



nematic state aq < 0, /3q > in (27c), yielding a phase 
with uniform nematic order and S = >/— q^q/(2/3q). We 
choose our coordinates so that this reference state has 
^ = 0. The phase transition line is shown in Fig. [15] as 
a function of density po and self- propulsion speed vq for 
the microscopic hard rod model discussed in [Baskaran 
and Marche"ttil ( |2008a|b] . 

As in the case of polar and nematic active systems. 




FIG. 16 (color online) Snapshots from simulations of self- 
propelled hard rods. The red dots mark the front end of the 
rods, i.e., the direction of self-propulsion. The rods are oth- 
erwise head-tail symmetric. The left image is for a density 
prodLlods — 0.7744, where Lrod is the length of the rod, and 
1/Pe = 0.005, where Pe = LrodVo/D\\ is the Peclet number, 
with vo the self-propulsion speed of each rod and the diffu- 
sion coefficient for motion along the long direction of the rod. 
The right image shows a close-up of clusters of sizes n=3, 10, 
and 22, highlighting the partially blocked structure and the 
"smectic-like" ordering of rods within a cluster. The clusters 
are chosen from a simulation with parameters 1/Pe = 0.00095 
and the same density as in the le ft image. Adapted with per- 
mission from Yang et al. (2010). 



Numerical simulations of collections of self-propelled 
rods with steric repulsion have recently revealed a rich be- 
havior, quite distinct from that of polar Vicsek-type mod- 
els. As predicted by theory, self-propelled rods with only 
excluded volume interactions do not order in a macro- 
scopically polarized state, but only exhibit nematic or- 
der, which appears to be long-ranged in two dimensions. 
Again, simulations have also confirmed that self- propul- 
sion enhances the tendency for nematic ordering as well 



as aggregation and clustering (Ginelli et al, 2010 Peru- 



ani et al., 2006 Yang et al., 2010). After an initial tran- 



sient, the rods form polar clusters that travel in a directed 
fashion. Large clusters can form by collisions of smaller 
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FIG. 17 Typical snapshots of nematic bands above the order- 
disorder transition. in the ordered phase. Arrows indicate the 
polar orientation of particles; only a fraction of the particles 
are s hown for clar i ty rea sons. Adapted with permission from 
Ref dGinelh et a/.| 120101). 



ones and break up due to collisions with other clusters 
or due to noise. Eventually the system reaches a station- 
ary state, in which the formation rate of any cluster size 
equals its breakup rate, and the rods aggregate in large 
stationary clusters. Above the order-disorder transition, 
phase separation manifests itself with the formation of 
nematic bands, consisting of high density regions where 
the particles are on average aligned with the long direc- 
tion of the band, but move in both directions, exhibiting 



no polar order ( ,Ginelli et al. 2010), as shown in Fig. 17 
Self-propulsion has also been shown to increase the seg- 
regation tendency in mixture of self-propelled rods with 



distinctly different motilities ( McCandlish et al., 2012). 



In closing this section we remark that mixed- symmetry 
models of the sort just presented may well be the best 
description of the collective crawling of rod-like bacteria 



such as those studied by Wu et al. (2011 ). Cluster forma- 
tion qualitatively similar to that observed in simulations 
of SP rods has been observed in recent experiments in 



myxobacteria Peruani et al. (2012) 



C. Current Status of Dry Active Matter 

Although much progress has been made in the under- 
standing and classification of dry active matter, a number 
of open questions remain. 

Clustering and phase separation, with associated giant 
number fluctuations, are ubiquitous in active systems. It 
was first suggested that giant number fluctuations may 
be a distinct property of the ordered state (nematic and 
polar) of active systems, intimately related with the ex- 
istence of a spontaneously broken orientational symme- 
try ( [Ramaswamy et "al 2003 Simha and Ramaswamy 
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FIG. 18 Cluster formation in myxobacteria mutants 
(SA2407) that cannot reverse their gliding direction. At long 
times the dynamical clustering process reaches a steady state 
that depends on cell density. The images show the clusters 
obtained at two different packing fractions 77 = pa, with p 
the two-dimensional cell density and a — AAjim^ the average 
area covered by a bacterium, with 77 = .16 (left) and 77 = 0.24 



(right). Adapted with permission from jPeruani et a/.| (2012). 



2002a Toner and Tu 1998). In ordered states such large 



fluctuations can indeed be understood as arising from 
curvature-driven active currents unique to the ordered 



state of active systems ( [Narayan et al. , 2007). More re- 
cently giant number fluctuations consistent with a stan- 
dard deviation growing linearly with the mean number 
of particles, rather than with the 1/2 exponent expected 
in systems where the central limit theorem applies, have 
been reported in active systems of symmetric disks with 
no alignment rule that do not exhibit any orientational 



broken symmetry (Fily et al., 2012). In this case large 



density fluctuations associated with strong clustering and 
phase separation seem to arise from the general mecha- 



nism proposed in Gates et al. (2010) and Tailleur and 



Gates] ( |200 8) and reviewed recently in Gates (2012) asso- 
ciated with the breaking down of detailed balance in sys- 
tems that are driven out of equilibrium by a local input of 
energy on each constituent, as in self-propelled systems 
or bacterial suspensions. Strong density inhomogeneities 
leading to persistent clustering have also been seen in 
layers of vibrated granular spheres (Aranson et al.^2m^ 



Prevost et al. 2004). Again, although in this case other 



mechanisms such as the inelasticity of collisions may play 
a role, more work is needed to elucidate the generic as- 
pect of this ubiquitous phenomenon. 



Although some controversy still remains ( Aldana et al. 



2007 Gonci et al. 2008 Vicsek et al. , 1995), there is now 



strong numerical evidence ([Ghate et al. , 2007) that the 



order-disorder transition in dry active systems is discon- 
tinuous, with an associated wealth of coexistence phe- 
nomena. One unusual aspect of the transition is that 
density fluctuations destabilize the ordered state right 
at the mean-field order transition. This behavior seems 
to be associated with a phenomenon that has been re- 



ferred to as "dynamical self-regulation" (Baskaran and 



Marchetti, 2012 Gopinath et al., 2012) associated with 



the fact that in this class of active systems the parameter 
that controls the transition, namely the density of active 
particles, is not tuned from the outside, as in familiar 
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equilibrium phase transitions, but rather it is dynami- 
cahy convected or diffused by non equihbrium active cur- 
rents controhed by the order parameter itself. It would of 
course be very interesting to attempt to use formal field 
theory and renormalization group methods to shed some 
light on this unusual non equilibrium transition. 

The behavior of collections of self-propelled or active 
particles is strongly affected by the presence of bound- 



tive stress distribution observed experimentally (Trepat 



aries and obstacles ( Elgeti and Gompper 2009 Wensink 



and Lowen 2008). Ratchet effects have been demon- 



strated experimentally and theoretically for active par- 
ticles interacting with asymmetric obstacles or on asym- 



metric substrates, without any external forcing ( Angelani 



et al] |2009[ [DiLeonardo et al] |2010| ) . For bacteria under- 
going run and tumble dynamics or self-propelled particles 
performing persistent random walks, the guiding effect 
of asymmetric boundaries can yield a variety of recti- 



fication effects (Galadja et a/., 2007 Wan et al 2008), 



that can even power submillimiter gears, as demonstrated 



in recent experiments (DiLeonardo et al 2010 Sokolov 



et al. 2010| ). Understanding the interaction of bacteria 
with obstacles and with passive particles is crucial for 
harnessing the collective power of these living systems 
and developing micron-scale mechanical machines pow- 
ered by microorganisms. Whether these phenomena can 
be described in the dry limit discussed in this section or 
require a proper treatment of flow as introduced in the 
next section and even hydrodynamic interactions is still 
an open question. 

Finally, there has been a surge of recent interest in 
the properties of dense active matter and the living crys- 
talline or glassy states that may be formed in these 
systems. The nonequilibrium freezing of active parti- 
cles may be directly relevant to the behavior of suspen- 
sions of self-propelled Janus colloids or other artificial 



microswimmers (Palacci et al, 2010). Recent work has 



demonstrated that active particles do form crystalline 
states, but the freezing and melting is in this case a 
true non equilibrium phenomenon that cannot be de- 
scribed simply in terms of an effective temperature for 
the system (Bialke et al, 20121). In vitro experiment on 



confluent monolayers of epithelial cells suggest that the 
displacement field and stress distribution in these living 
systems strongly resemble both the dynamical hetero- 
geneities of glasses and the soft modes of jammed pack- 
ings (|Angehni eTfl/.j |2Q1Q| |201l] [Petitjean et al] |201q1 



et al. , 2009p in living cellular material. 



III. ACTIVE GELS: SELF-DRIVEN POLAR AND APOLAR 
FILAMENTS IN A FLUID 

In this section we describe an alternative method for 
constructing hydrodynamic theories for a class of active 
materials. This approach involves a systematic deriva- 
tion of the hydrodynamic equations based on a gener- 
alized hydrodynamic approach close to equilibrium fol- 
lowing closely the work of Martin Parodi, and Pershan 



for nemato-hydrodynamics ( de Gennes and Prost 1993 
Martin et al. 1972[). As in the previous sections, the the- 



ory is mostly based on symmetries and does not involve 
significant microscopic considerations; it is thus applica- 
ble to a whole range of systems, which share the appro- 
priate polar or nematic symmetries and which are liquid 
at long times. We focus here on an 'active gel', defined as 
a fluid or suspension of orientable objects endowed with 
active stresses, with momentum damping coming from 
the viscosity of the bulk fluid medium, rather than from 



m 



friction with a substrate or a porous medium (Jiilicher 
et a/.[|2007[ |Kruse et al.\ |2004| |Simha and Ramaswamy[ 
2002a|b ) . The equations that emerge are those proposed 



Simha and Ramaswamy ( 2002a ) for self-propelling or- 



ganisms, but the development in Jiilicher et al. (2007) 



and Kruse et al. ( 2004 ) was carried out in the context of 



the cytoskeleton of living cells, a network of polar actin 
filaments, made active by molecular motors that consume 
ATP. 



In section III.B we consider an active system with a 
polarization field p but we suppose that the physics of 
this system is invariant under a change of p to —p. We 
therefore consider an active nematic gel. The polariza- 
tion vector can in this case be considered as the direc- 
tor field for the nematic order. An alternative approach 
is to describe the ordering not by a director field but 
by a nematic alignment tensor, Q. This approach, also 



part of the phenomenological treatment in Simha and Ra- 



maswamy (2002a), is presented in Salbreux et al. (2009). 



A brief discussion of polar active gels and of the effects 



of polarity on the dynamics is given in section III.C 



Poujade et al, 2007; Trepat et al, 2009). This observa- A. Hydrodynamic equations of active gels 



tion has led to new interest in the study of active jammed 
and glassy states (Henkes et al. , 2011) obtained by pack- 



ing self-propelled particles at high density in confined 
regions or by adding attractive interactions. In these 
models the interaction with the substrate has so far been 
described simply as frictional damping. Although a rich 
and novel dynamical behavior has emerged, it has become 
clear that a more realistic description of stress transfer 
with the substrate will be needed to reproduce the ac- 



For the sake of simplicity in this section we only dis- 
cuss a one component active polar gel considering there- 
fore that the complex composition of active materials 
such as the cytoskeleton can be described by an ef- 
fective single component. The original formulation of 



Simha and Ramaswamy (2002a) for the hydrodynamics 



of self-propelled orientable suspensions was already an 
explicitly multicomponent formulation of active liquid- 



20 



crystal hydrodynamics, but without a formal link to the 
nonequilibrium thermodynamic approach to active sys- 
tems of Jiilicher et al. (2007). Generalizations of the 



law reads 



latter approach to multicomponent systems are possible 
and have been recently proposed by [Callan- Jones and| 



Jiilicher (2011) and Joanny et al. (2007). The multicom- 



ponent theory takes properly into account the relative 
permeation between the various components which are 
neglected in the simpler one-component description. As 
in the previous sections, we retain as slow variables, the 
number density p, the polarization p and the momentum 
density g = pmv where v is the local velocity of the gel 
and m the (effective) mass of the molecules. 



1. Entropy production 

The derivation of generalized hydrodynamic equations 
is based on the identification of fl uxes and forces from 
the entropy production rate S a la de Groot and Mazur 
(1984). We do not wish to consider heat exchange here; 



we will assume that the active gel has a constant tem- 
perature, meaning that it is in contact with a reservoir 
at a finite temperature T. In this case, the entropy pro- 
duction rate is related to the rate of change in the free 
energy of the active gel TS = — ^ • 

For a passive gel at rest, the free energy density / is 
a function of the two intensive variables, the density p 
and the polarization p and its differential is df = jidp — 
hadpa- The field conjugate to the density is the chemical 
potential p and the field conjugate to the polarization is 
the orientational field h. For a passive system moving at 
a velocity v, the density of kinetic energy ^pmv^ must 
be added to the free energy. 

For an active gel one must also take into account the 
fact that energy is constantly injected into the gel locally. 
A simple intuitive way to introduce the energy injec- 
tion is to assume that as in the case of the cytoskeleton, 
this is due to a nonequilibrium chemical reaction such 
as the consumption of ATP. If the energy gain per ATP 
molecule is denoted by A/i, and the rate of advancement 
of the reaction (the number of ATP molecules consumed 
per unit time and unit volume) is denoted by r, then 
the associated rate of change of the free energy per unit 
volume is —rAja. Taking into account all contributions, 
we find the entropy production rate of an active gel at a 



constant temperature T ( Kruse et al. , 2004 , 2005 ) 



TS 



d 1 2a dp 
--(-pmv)-p^ 



-haPa^rAp} . (29) 



2. Conservation laws 

The two conserved quantities in an active gel are the 
density and the momentum. The density conservation 



dp 
di 



V • (pv) = . 



(30) 



The momentum conservation law can be written as 



dt 



(31) 



where the momentum flux in the system is 11^^^ = 
pmVaVis — crj^^, where the first term is associated with 
the so called Reynolds stress and aj^,^ is the total stress 
in the system. For most active gels, in particular for bi- 
ological systems, the Reynolds number is very small and 
we ignore in the following the Reynolds stress contribu- 
tion. 



3. Thernnodynannics of polar systems 

For generality in this section we consider active gels 
in three dimensions. The polarization free energy of an 
active polar material is a functional of the three compo- 
nents of the polarization vector p. However, if the system 
is not in the vicinity of a critical point there are only two 
soft modes associated with rotations of the polarization. 
The modulus of the polarization is not a hydrodynamic 
variable and is taken to be constant. Without any loss 
of generality, we assume that the polarization is a unit 
vector, and that the effect of its modulus is integrated in 
the phenomenological transport coefficients. The polar- 
ization free energy of the active gel is then the classical 
Frank free energy of a nematic liquid crystal (de Gennes 



and Prost 1993) 



(V-p)2 



(p-(Vxp))2 



^(px(Vxp))2 + ^;V-p-i/.i;p2 



(32) 



The first three terms correspond to the free energies of 
splay, twist and bend deformations. The three Frank con- 
stants Ki are positive. We have also added in this free 
energy a Lagrange multiplier /ly to insure that the polar- 
ization is a unit vector. This free energy is very similar 
to the free energy of Eq. (|4| , although in Eq. (|4| we had 
made the additional approximation that the Frank con- 
stants are equal. In the case where the polarization is a 



critical variable, one would also need to add to Eq. (32) a 



Landau expansion in powers of the polarization modulus 
as done in Eq. (|4|. 

The orientational field is obtained by differentiation of 
the free energy ( [32| ). It is useful to decompose it into a 
component parallel to the polarization h\\ and a compo- 
nent perpendicular to the polarization h± 

In the simple case where the system is two-dimensional, 
the polarization can be characterized by its polar angle 
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6. In the approximation where the Frank constants are 
equal the perpendicular molecular field is = KV^O. 

In a non-isotropic medium the stress is not symmet- 
ric. There is an antisymmetric component of the stress 
associated to torques in the medium. As for nematic 
liquid crystals, this anti-symmetric component can be 
calculated from the conservation of momentum in the 



fiuid (de Gennes and Prost, 1993) and is equal to cr^. 



4. Fluxes, forces and time reversal 

Using the conservation laws and performing integra- 
tions by parts, the entropy production can be written 
as 



(33) 



TS = dv{aal3Val3 + Paha + ^A/i} , 



flux 


force 








ha 


r 


A/i 



TABLE II Fluxes and associated forces controlling the en- 
tropy production in a one-component nematic active gel. 



care must be taken in considering the time reversal sym- 
metry. The fiuxes must all be separated into a reactive 
component with a signature opposite to that of the con- 
jugate force and a dissipative component with the same 
signature as the conjugate force. As an example, the re- 
active component of the stress is the elastic stress and 
the dissipative component is the viscous stress. Only the 
dissipative component of each flux contributes to the en- 
tropy production. 



where Gajs is the symmetric deviatoric stress tensor de- 
fined by 



(34) 



with cr^^ the antisymmetric part of the stress tensor and 
P the pressure. We have also introduced the strain rate 
tensor Vcxj3 and the antisymmetric part of the velocity 
gradients associated to the vorticity cJq,^, defined as 



Finally, Pa = with 
Dpa 



Dt 



dpa 
dt 



V(3df3Pa + ^a^P^ 



(35a) 
(35b) 

(36) 



is the convected derivative of the polarization. Note that 
if the system is chiral, kinetic momentum conservation 
must be taken into account properly and the local rota- 



2012) 



tion is no longer given by 00^/3 (Furthauer et al. 

This form of the entropy production allows for the 
identification of three forces: v^^ which has a signature 
— 1 under time reversal, which has a signature +1 un- 
der time reversal and A/i which also has a signature +1 
under time reversal. The conjugated fiuxes are respec- 
tively (JociB^ Pa and r. The fiuxes and associated driving 
forces that control the hydrodynamics of active gels are 
summarized in Table JIl 

In a linear generalized hydrodynamic theory, the con- 
stitutive equations of the active gel are obtained by writ- 
ing the most general linear relation between fiuxes and 
forces respecting the symmetries of the problem, such as 
translational and rotational symmetries with one vector 
Poi and one tensor ^c^/? = PcxPi3 — \^cxi3 in 3(i. Particular 



B. Active nematic gels 

1. Constitutive equations 

We first consider an active polar liquid for which the 
relationship between fiuxes and forces is local in time. 

It is convenient to split all tensors into diagonal and 
traceless parts part: (Ja^ = crSais + ctq,^, with a = 
(l/3)crQ,Q,, a^a = and d the dimensionality. Similarly, 
we let Va(3 = f + whcrc u = d^Vj is the diver- 
gence of the velocity field. Finally, all fiuxes are written 
as the sums of reactive and dissipative parts. 



cr, 



P _ 

r = r'' + r"^ . 



a(3 

p:- 



pd 



(37a) 
(37b) 
(37c) 



a. Dissipative fluxes. Only fiuxes and forces with the 
same time signature are coupled and (1^13 is only coupled 
to Va(3- This leads to the constitutive equations 

d 



(7 

;^d 



TjU 



(38a) 
(38b) 



We ignore here the tensorial character of the viscosity 
and assume only two viscosities, as for an isotropic fiuid: 
the shear viscosity r] and the longitudinal viscosity fj. 
The two other fiuxes are coupled and the corresponding 



constitutive equations read (Kruse et 

pd ^ 
" 71 
r*^ = AAfi + epaha 



2004 2005) 



epaAij , 



(39) 
(40) 



7i is the rotational viscosity and we use here the Onsager 
symmetry relation which imposes that the "dissipative" 
Onsager matrix is symmetric. 
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b. Reactive Buxes The reactive Onsager matrix is anti- 
symmetric and couples fluxes and forces of opposite time 
reversal signatures. 



-(A/j.^TyiPaha , 



(41a) 



-CAfiqaf3 + ^P(3ha- ^^7^7^^^) , 



— U 



(41b) 
(41c) 

(41d) 



In most of the following we consider incompressible flu- 
ids so that 1^ = V • V = 0. In this case the diagonal 
component of the stress can be included in the pressure 
which is a Lagrange multiplier ensuring incompressibility 
and one can set ( = ui = f] = 0. 

To summarize, the hydrodynamic equations for an in- 
compressible one-component active fluid of nematic sym- 
metry are given by 

mp (S^ + V • V) V = - VP + V • tr , (42a) 

(St + V • V) Pa + UJa(3P(3 = -J^lVa(3P(3 + — + ^A/ipa , 

7i 

(42b) 

to be supplemented with the incompressibility condition 
V • V = 0. Assuming the Frank constants are all equal 
to the molecular field is given by 



ha = KV'^Pa + h^\Pa , 



(43) 



with /ly a Lagrange multiplier to be determined by the 
condition |p| = 1. Finally, it is convenient for the fol- 
lowing to write the deviatoric stress tensor given by the 
sum of trace and deviatoric parts of the dissipative and 
reactive components by separating out passive and active 
parts as 

= CtI^ + CTa^ , (44) 

with passive and active contributions given by 

= 2r]Va^ + y ^P^ha - ^Pjh^Sa^) , (45a) 

<^ = -CA^qa^ . (45b) 

In many biological applications inertial effects are negligi- 
ble and the Navier-Stokes equation (42a) can be replaced 



by the Stokes equation obtained by simply neglecting all 



inertial terms on the left hand side of Eq. (42a) and cor- 



responding to a force balance equation, given by 

- VP -}- V • o- = . (46) 



2. Microscopic interpretation of the transport coefficients 

The Onsager approach that we have outlined intro- 
duces several transport coefficients. Some of these co- 
efficients exist for passive nematic liquid crystals such 
as the viscosities, the rotational viscosity 71 or the flow 
-coupling coefficient ui. The two important new coeffi- 
cients are the transport coefficients associated with the 
activity of the system, e and (. The active stress in the 
system is aa = — CA/i. In the case of the cytoskeleton 
this can be viewed as the stress due to the molecular 
motors, which tend to contract the gel. The sign of the 
activity coefficient ( is not imposed by theory. A negative 
value corresponds to a contractile stress as in the actin 
cytoskeleton. A positive value of ( corresponds to an ex- 
tensile stress as observed in certain bacterial suspensions. 
Conceptually, active stresses in living matter were first 
discussed by Finlayson and Scriven (1969). These au- 



thors, however, specifically steer clear of uniaxial stresses 
arising from motor-protein contractility. The first incor- 
poration of self-propelling stresses into the generalized 
hydrodynamics of orient able fluids was by |Simha and Ra- 
although it has long been understood 



maswamy 



(2002a 



that the minimal description of a single force-free swim- 



and Kessler 1992). 



mer is a force dipole (Brennen and Winet, 1977 Pedley 



The other active coefficient, e, is an active orientational 
field that tends to align the polarization when it is posi- 
tive. In the limit where the modulus of the polarization 
is p = 1 one can always consider that e = and introduce 
an effective activity coefficient ( + eji vi . In the following 
we therefore choose e = and use this effective value of 
the activity coefficient C,. 



3. Active currents in nematic and polar systems from forces 
and fluxes 

Let us take the opportunity to show that the active 



currents in nematic (section II. B) and polar (section 
II. A) systems follow naturally from the forces-and- fluxes 



framework. In this short treatment we ignore fluid flow. 
Consider a mesoscopic region in our active medium, in 
a chemical potential gradient V<^> corresponding to the 
concentration p, and subjected to a nonzero chemical po- 
tential difference Aji between a fuel (ATP) and its reac- 
tion products ( ADP and inorganic phosphate) . For small 
departures from equilibrium, the fluxes r (the rate of con- 
sumption of ATP molecules) and J must be linearly re- 
lated to A/i and V$. The presence of local orientational 
order in the form of p and Q allows one to construct 
scalars p • V<l> and V • Q • V<l>; thus r in general gets a 
contribution (CpP + CqV ' Q) ' V*^? where and C,q are 
kinetic coefficients depending in general on p and other 
scalar quantities. The symmetry of dissipative Onsager 
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coefficients then implies a contribution 

J active = (CpP + CqV • Q)A/X 



(47) 



to the current. In the presence of a maintained constant 
value of A/i the current ([47|) rationahzes, through the ("p 
and terms, the form ofpal (with (pAjj^vop) and (25) 
(with Cq^I^ ~^ Cq)' particular, it underlines the fact 
that active currents do not require an explicit polar or- 
der parameter. Even without p, the term in (q in (47) 
tells us that curvature in the spatial arrangement of ac- 
tive filaments gives rise to particle motion. In retrospect 
this is not shocking: a splayed or bent configuration of 
a nematic phase has a vectorial asymmetry, as argued 
through Fig. 



12 



In a system out of equilibrium this 
asymmetry should reflect itself in a current. As shown by 
Ramaswamy et al. (2003) and discussed in section 



II.B 



leads to giant number fluctuations in active nematics. 



4. Viscoelastic active gel 

An active gel is not in general a simple liquid but 
rather a viscoelastic medium with a finite viscoelastic re- 
laxation time, which is only liquid at long time scales. In 
a passive visco-elastic medium the constitutive relation 
between stress and strain is non local in time. The sim- 
plest description of a visco-elastic medium is the so-called 
Maxwell model where the system only has one relaxation 



time (Larson, 1988). Within this model the constitutive 



equation is 



Dt 



(48) 



The Maxwell model involves two material constants, the 
viscoelastic relaxation time r and the shear modulus E. 
The long time viscosity of the medium is then 77 = Et. 
In order to respect rotational and translational invari- 
ance, we use here a convected Maxwell model with a 
convected time derivative of the stress tensor = 

^^1^ + v^d^aa(3 + ^Q!7^7/3 + Notc that there 

are several ways of defining the convective derivative of 
tensors. For simplicity we use the same notation to 
denote convected derivatives of vectors and tensors. 

The generalization of the Onsager hydrodynamic ap- 
proach of the previous section to a viscous elastic polar 
passive medium lead to the constitutive equations for an 



active gel (Jiilicher et al. 2007). 



(49a) 

= — (1 + ^^)ha - J^lVa(3P(3 , (49b) 

JJt 7i JJt 

where for simplicity we have only considered an incom- 
pressible active gel where the modulus of the polarization 



is unity. Note that the memory of the system not only 
plays a role for the stress but also for the dynamics of 
the orientation and that we have supposed that the two 
corresponding relaxation times are equal. 

It is important to note that the dynamical equation 
for the polarization is very similar to Eq. (|2b|) used in 



section [n] which has been obtained using the same sym- 
metry arguments; Eq. (2b) ignores memory effects and 



therefore the visco-elasticity of the polarization response; 



all the extra terms in Eq. ( 2b ) do not appear here because 



the Onsager approach that we use only derives the linear 
hydrodynamic theory. 



C. Active polar gels 

For the sake of simplicity, we have only presented here 
the derivation of the hydrodynamic theory of active gels 
in the simplest case where the system is a single compo- 
nent fluid and has nematic symmetry, p being the direc- 
tor, and ignore any type of noise. Several extensions of 
this theory have been proposed. 

For a polar system, there is an extra polar term in the 
free energy (32), Fp = J dr v{p) {\/ - p). This sponta- 
neous splay term is a surface term if the coefficient v is 
a constant. If v depends on the local density, this term, 
that was included in Eq. (|4|, yields a "pressure-gradient" 
proportional to Vp in the equation for the polarization. 
Other non linear polar terms that can be added on the 
right hand side of the dynamical equation for the polar- 
ization (49) are proportional to p • Vp, Vp^ and pV • p. 



These terms have already been considered in Eq.(2b). 



The first of these terms cannot be derived from a free en- 
ergy and is therefore an active term proportional to A/i. 
Its effect in a dynamics linearized about an ordered state 
was considered for active liquid-crystalline suspensions 
by Simha and Ramaswamy (2002a). The two other terms 



can be derived from a free energy and have both active 
and passive contributions. Within the Onsager linear 
hydrodynamics scheme, other polar terms in the equa- 
tions show up only at nonlinear order or at sub-leading 
order in a gradient expansion. Active stresses propor- 
tional to Ajii{diPj + 9jPi) unique to polar fluids are ob- 
tained from the microscopic theory (Marchetti and Liv- 



erpool[ |2007[ ) , but are considered nonlinear in the driving 
forces in the context of the Onsager approach. The ef- 
fect of these polar terms has been studied in detail by 
Giomi and collaborators (Giomi et al, 2008). In gen- 



eral these polar terms are important when describing an 
active suspension as opposed to the one-component sys- 
tem considered here. In this case these terms yield spa- 
tial inhomogeneities in the concentration of active par- 
ticles dGiomi and Marchetti} |2012[ [Tjhung et~al] |2QllD 



that are not obtained in active suspensions with nematic 



symmetry (Giomi et a/., 2011, 2012). 



The effect of noise in active polar gels can be studied by 



24 



introducing random Langevin forces in the constitutive 
equations. For treatments including thermal and active 
non-thermal noise sources in a systematic way, se e |Basu 
et a/.|(|2008|);|Chen et a/.| (|2007|); iHatwalne et a/.| (|2004|); 



of simple elements. Several qualitative features of these 
experiments including the defect patterns could be sim- 
ply understood using continuum models describing the 
polar orientation of the rigid filaments and the density of 



Lau et al. (2003); Lau and Lubensky (2009); and Sarkar processive motors (Lee and Kardar, 2001 Sankararaman 



and Basu (2011 ). The precise description of the statistics 



et 



2004) 



of active noise requires a microscopic description of the 
gel which is not generic and goes beyond the scope of this 



A detailed study of the nature of defects and their 
phase transitions within the framework of the active 



review ( |Basu et al\ |2008| ). [Hatwalne et al\ ( |2004[ ) intro- gel theory was done by |Kruse et al\ (|2004[ |2006a|) who 



duce the active noise in the isotropic phase of an active 
gel through the fluctuating active stress in the Navier- 
Stokes equation, and use scaling arguments to estimate 
the apparent temperature it would generate in a tracer 
diffusion measurement. 

Most active systems described in the earlier sections of 
this review are multicomponent systems containing a sol- 
vent and active objects. In many instances these systems 
can be described by an effective one-component theory as 
described here. However, in particular when considering 
viscoelastic effects the one-component theory ignores the 
permeation of the solvent through the active gel. A de- 
tailed two component theory of active gels that takes into 
account properly permeation effects is given by ICallaiT 
Jones and Jiilicher (2011). 



D. Active Defects 

Ordered phases of active matter, like their counter- 
parts at thermal equilibrium, should exhibit topological 
defect configurations, generated either through specific 
boundary conditions or spontaneously in the bulk. As in 
equilibrium systems, the nature of these defects should 
depend on whether the system has polar or apolar sym- 



metry (de Gennes and Prost 1993 Kung et 



2006 



Palffy-Muhoray et al. 1988). The selection criterion for 



the defect strength in active systems is not obvious, as 
one cannot a priori invoke free-energy minimization as 
in equilibrium systems. However, it does appear that 
experiments see strength +1 and strength +1/2 defects, 
respectively, in polar ( ,Nedelec et al\ 1997) and apolar 
(Narayan et al, 2007) active systems. Activity confers a 



particularly interesting property on defects, namely rota- 



tional (Kruse et al. 2006a Nedelec et al., 1997) or trans- 



lational (Dogic, 2012 Narayan et al. 2007) movement. 



with sense or direction determined by the chirality or 
polarity associated with the defect. 

The first quantitative experiments which explicitly 
demonstrated how active mixtures of long rods (mi- 
crotubules) and motors (kinesin) could spontaneously 
form defects such as asters and spirals were reported in 
Nedelec et al\ ( |1997D and [Surrey aT\ ( |200lD ; see Fig. 
2] These patterns showed a remarkable resemblance to 
the microtubule-based spindle patterns in the cell, thus 
suggesting that the gross features of spindle patterning 
could be understood as arising from a self-organization 



showed that flows arising from active stresses leads to 
systematic rotation of chiral defects. We summarize here 



the calculation, working as in Kruse et al. (2004) with 
the ordered state described by a vectorial order param- 
eter p in two dimensions. Polarity enters nowhere in 
the analysis, except as justification for working with a 
strength +1 defect. We parametrize a two-dimensional 
defect configuration of the polarization field p (of unit 
magnitude) with topological charge ±1 using polar co- 
ordinates (r, 6>). Thus charge +1 defects such as asters^ 
vortices and spirals^ may be represented by an angle ?/^, 
with pr = cosT/^ and pe = sinip^ such that t/j = (or tt) is 
an aster, i/j = ±7r/2 a vortex and i/j = i/jq (any other con- 
stant) a spiral. At equilibrium, a situation corresponding 
to defect configurations in a ferroelectric nematic liquid 
crystal, the optimal value of tjj is obtained by minimiz- 
ing the Frank free-energy functional F ( 32 ) with K2 = 



(no twist as we are working in two dimensions). With 
appropriate boundary conditions, it easy to see that the 
stable defect configurations are (i) asters when Ki < Ks, 
so that splay is favored, (ii) vortices when Ki > i^3, i.e., 
bend is favored, and (iii) spirals when Ki = K^. In an 
active system, however, the stability of defect configura- 
tions is obtained by solving the dynamical equations for 
the polarization p together with the conditions of force 



balance and overall incompressibility, as in section III.B 



Suppose Ki < i^3, so that the aster is stable in the 
absence of activity, A/i = 0. Now introduce activity. 
Linear stability analysis shows that at sufficiently large 
contractile active stresses (Afi < 0, the aster gets desta- 
bilized giving rise to a spiral with an angle i/jq set by 
the flow alignment parameter ui (assuming stable flow 
alignment). An entirely similar analysis, starting from a 
stable vortex for Ki > i^3, with A/i = 0, shows again an 
instability for large enough A/i. The reason this happens 
here and does not happen for systems without activity 
is that the active stresses associated with the perturbed 
director configuration give rise to flows whose effect on 



the director is to reinforce the perturbation. Figs. [19] and 
2Q| illustrate the flow associated with the spiral instabil- 
ity and the stability domains of the defects, respectively. 
The active spiral has a sense of direction and will there- 
fore rotate. The angular speed can be obtained by solving 
the steady state equations for i/j and vq, leading to 



vo{r) = ujQr log 



(50) 
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FIG. 19 (color online) Sketch of a rotating spiral defect in an 
active nematic fluid with vanishing elastic anisotropy 6K — 0. 
The gray arrows mark the director, and the solid lines show 
its spiral structure. The hydrodynamic velocity field is in the 
azimuthal direction, and its profile is indicated by the dark 



arrows. Adapted with permission from ^Kruse et al. (2004). 




FIG. 20 (color online) Stability diagram of topological defects 
in an active nematic fluid. For small values of the activity C^M 
one gets asters and vortices, respectively for 5K — Ki — K3 < 
and > 0. When of sufficiently large magnitude both 
asters and spirals are unstable to the formation of a rota ting 
spiral. Adapted with permission from Kruse et al (2004). 



tra term vi J d^rpV • p in the free energy (32), with the 
form of a spontaneous splay that depends on the local 
concentration, and restore the concentration equation, 
dtp = — V- J, where the filament current J = vopp—DVp, 
has an active advective and a diffusive contribution. Now 
for large enough A, the defect configurations are gen- 
erated by the internal dynamics (and insensitive to the 
boundary for large systems) and so have finite size. When 
advection is negligible, the defect size is set by the ratio 
of the spontaneous splay coupling strength ^4 to a Frank 
modulus, as it would be in an equilibrium polar liquid 
crystal. On the other hand, when advection is apprecia- 
ble so as to have density clumping in regions where p 
points inwards towards a common center, the defect size 
is set by the ratio of diffusion to active advection, D/vq. 
We note that since the advection current J oc pp, a small 
perturbation of the vortex configuration renders it unsta- 
ble and the only stable defects are inward pointing asters 



and spirals ( Gowrishankar and Rao, 2012). 



Since defects now have a finite size, it is possible to 
have an array of defects which interact with each other. 
Such studies have shown that under certain conditions, 
one obtains a stable lattice of asters accompanied by a va- 



riety of phase transitions (Gowrishankar and Rao 



Voituriez et a/., 2006 Ziebert and Zimmermann 



2012 



2005). 



A detailed study of defect-defect interactions and the dy- 
namics of defects and their merger in this active context 
are open problems for the future. 

So far we have discussed charge 1 defects generated 
in polar active media. Apolar active media, described 
by a local orientational tensor Q, exhibit ±1/2 strength 
disclinations, topologically identical to those obtained in 
equilibrium nematic liquid crystals (Fig. 7). The orien- 
tation field around a defect of strength +1/2 has a po- 
larity, whereas that around a —1/2 has a 3-fold sym- 
metric appearance. On general grounds, the +1/2 de- 
fect should move spontaneously whereas the —1/2 defect 
should show no such tendency. Precisely this behavior 
seems to be observed in the active nematic phase in a 



vibrated granular-rod monolayer (Narayan et al. 2007). 



where, as is inevitable on dimensional grounds, cjq scales 
as the ratio of the active stress to a viscosity, with a 
detailed form that includes a dependence on the director 
kinetic coefficient and the flow-alignment parameter. In 
a flnite system of size imposing a vanishing velocity 
at the outer boundary due to the presence of a wall, we 
can set the length scale ro = R. 

In the above analysis we ignored the dynamics of the 
concentration fleld p. This is probably acceptable when 
the active units are long and rigid or when their concen- 
tration is so high that excluded volume considerations do 
not permit significant density inhomogeneities. We can 



include p (Gowrishankar and Rao 2012| ) through an ex- 



E. Current status on active gels 

The nonequilibrium thermodynamic description of ac- 
tive systems is a systematic approach based on symme- 
tries and in particular on invariance against time inver- 
sion. It is however based on a linear expansion of fluxes 
in terms of forces and can in principle only describe sys- 
tems close to equilibrium where Ap tends to zero. One of 
the main application though is to biological systems that 
are mostly far from equilibrium systems. There is no sys- 
tematic extension of the theory to systems far from equi- 
librium. One must rely either on a microscopic descrip- 
tion that generates non-linear contributions by coarse- 
graining to large length scales and long time scales or on 
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experimental results that emphasize specific non-linear 
aspects, which can then be introduced in the theory. Mi- 
croscopic or mesoscopic descriptions of molecular motors 
are a good example of the first case and lead to motor 



Active gel models have also been used to describe cross- 



forces or velocities that are not linear in A/i ( Jiilicher and 



Prost, 1997 Liverpool et al. 



linked motor-filament systems that behave 


as soft solids 


at large scales (jBanerjee and Marchetti| 


2011| 


Levine 


and MacKintoshf 20091 |MacKintosh and 


^jcvine 


|2008| 


Yoshinaga et al.\ 2010') and also to understand the ori- 



Several experiments gin of sarcomeric oscillations, both from microscopic 



suggest that the treadmilling associated to the polymer- 
ization and depolymerization of actin in a cell depends 
on the force applied on the filaments or the local stress in 



a non-linear way (Mogilner and Oster, 1999 Prost et al. 



2007). The force- dependent treadmilling is essential for 
many cellular processes such as cell migration or cell ad- 
hesion and a full description including these effects in the 



active gel theory has not been proposed yet ( Keren et al. 



2008) 



Another intrinsic limitation of the current active gel 
theory is the assumption of linear rheology and the use of 
the Maxwell model with a single relaxation time. A large 
body of experimental work shows that for many types of 
cells as well as for actomyosin gels there is a broad distri- 
bution of relaxation times leading to a complex modulus 
decreasing as a power law of frequency with a small ex- 
ponent between 0.1 and 0.25 (Fabry et al., 2001). An 
ad hoc power law distribution of relaxation times could 
be introduced in the active gel hydrodynamic theory but 
despite some efforts, this type of law is not understood 
on general grounds (Balland et al. , 2006). The non-linear 



rheology of actin networks has also been studied in detail 
and actin networks in the absence of molecular motors 



have been found to strain thicken (Gardel et al. 2004). 
This behavior can at least in part be explained by the 
inextensibility of the actin filaments. In the presence of 
myosin motors, the active stress induced by the molecu- 



lar motors itself stiffens the active gel ( Koenderink et al. 



2009). 



The hydrodynamic description of active polar or ne- 
matic gels is very close to that of nematic liquid crystals. 
Nevertheless the existence of an active stress leads to sev- 
eral non-intuitive and spectacular phenomena. The most 
spectacular is the flow instability described below that 
leads to spontaneously flowing states. Other spectacular 
effects are associated with active noise in these systems. 
In all active systems the noise has a thermal component 
and a non thermal active component. The properties of 
the active noise cannot be inferred from the macroscopic 
hydrodynamic theory and must be derived in each case 
from a specific microscopic theory. The study of tracer 



diffusion in a thin active film for example (Basu et al. 
[2011) leads to an anomalous form of the diffusion con- 
stant that does not depend on the size of the particle 
but only on the thickness of the film. We believe that 
there are still many unusual properties of active gels to 
be discovered and that in many cases, this will require 
detailed numerical studies of the active gel hydrodynamic 



equations such as the one performed in Marenduzzo et 



(2007). 



as well as from a continuum viewpoint ( 


Banerjee and 


Marchetti, 2011 


Giinther and Kruse, 2007 


). Further, ac- 



tive gel models have been shown to successfully account 
for the experimentally observed traction stresses exerted 



by cells and cell sheets on compliant substrates (Banerjee 
and Marchetti||2011| [Edwards and Schwarzl|2011||Mertz 



et 



2012). 

Finally, a large part of the theoretical activity on ac- 
tive gels aims at a quantitative description of biologi- 
cal phenomena at the scale of the cell and of cellular 
processes involving the cytoskeleton. Some success has 
already been obtained in discussing lamellipodium mo- 



tion ( Kruse et al. , 2006a ) or the formation of contractile 



rings during cell division (Salbreux et al. 2009). An im- 



portant step is the connection of the parameters of the 
hydrodynamic theory with the more microscopic parame- 
ters that can be monitored experimentally which requires 
an explicit coarse graining of the microscopic theories as 
discussed in section [Vl At larger scale one can build a hy- 
drodynamic theory of tissues that shares many features 



with the active gel theory described here (Ranft et al. 



2010). 



IV. HYDRODYNAMIC CONSEQUENCES OF ACTIVITY 

In this section we describe a number of remarkable hy- 
drodynamic phenomena induced by activity. Most of the 
section is devoted to the description of materials prop- 
erties of active gels, such as thin film instabilities and 



rheology. In section IV. D| we also highlight some of the 
remarkable successes of the hydrodynamic theory of ac- 
tive gels in describing phenomena observed in living cells. 
A more complete review of these latter class of phenom- 



ena can be found in Joanny and Prost (2010). 



A. Instabilities of thin liquid active films 

1. Spontaneous flow of active liquid films 

One generic property of active orient able liquids, 
whether polar or apolar, is the instability of any homoge- 
neous non- flowing steady state toward a inhomogeneous 
spontaneously flowing state, as shown by |Simha and Ra^ 
maswamy (2002a). We illustrate this instability in the 



simple geometry of a thin active nematic liquid film ( Voi- 



turiez et al., 2005). 

We study a thin film of thickness h supported by a 
solid substrate. For simplicity we only consider the two- 
dimensional geometry sketched in Fig. [21] where the sub- 
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Lc<L 



\\\\\\\\\\\\\\\\\\\\\\ 



FIG. 21 (color online) Spontaneously flowing film of active 
liquid. 



strate is along x and the normal of the film is along 
y. We choose anchoring conditions parallel to the film 
surface so that for y = O^h the polarization is along x, 
Px = Py = 0. An obvious solution for the equations of 
motion of an active liquid is a non flowing state v = 
with a constant polarization parallel to x. We now dis- 
cuss the stability of this steady state. We look for a state 
of the system which is invariant by translation along x 
so that all derivatives with respect to x vanish and with 
a polarization that is not parallel to the film surfaces 
Px = COS 0{y)^py = sin 0{y). The velocity is along the 
X direction and the shear rate tensor only has one non 
vanishing component Vxy = w = \dyVx' 

The total stress can be written as = —PSaf3 + 
^ o'a/s where the last term is the deviatoric stress 



given by the constitutive equations (|49[) and the previ- 



ous term is the antisymmetric component of the stress. 
Force balance in the film is written as dydy^ = 0. Taking 
into account the fact that on the free surface the shear 
stress vanishes, we obtain = 0. Using the constitu- 
tive equation, this leads to 

-h^ = 4r]w-CAiJ.sm20^iyi{h\\sm20^h±cos20) , (51) 

where we have introduced the parallel and perpendicu- 
lar components of the orientational field h\\ = hx cosO + 
hy sin 6 and h± = hy cos 6 — hx sin 6 

The constitutive equations for the polarization of 



Eq. (49) give 



hx 

uiw smt 

71 



— — uiw cos c 

71 



(52) 



Combining Eqs. (51) and (52), we obtain the perpen- 



dicular field and the velocity gradient 

CA/isin2i9(l + cos 20) 
' + 2ui cos 26 



hi = 



71 



1 



W = 



CA/i sin 2(9 



4?7 
71 



21^1 cos 20 



(53) 



In the approximation where the Frank constants are 
equal, the perpendicular molecular field is h± = 



5F_ 
56 



K 



dy 



If the angle 6 is small, an expansion of Eq.(53) 



for the perpendicular field to lowest order in 6 gives 
j2 = where the characteristic length L is de- 
fined by 



^ ^ -2CA/i(l + ^i) 

+ (1 + ^^)2)] 



(54) 



Assuming 1 + z^i > 0, if the activity coefficient ( is neg- 
ative corresponding to a contractile stress, > and 
the length L is indeed real. The polarization angle varies 
then as ^ = ^0 sin (|;). This satisfies the anchoring con- 
dition on the solid surface ^ = but the anchoring con- 
dition on the free surface y = h can only be satisfied 
if h = ttL. If the film is thin h < irL = Lc there is 
no solution with a finite and the non-flowing steady 
state is stable. For a film thicker than the critical value 
h = ttL a solution with a finite value of the polarization 
angle exists and the non-fiowing steady state is unsta- 
ble. If h > Lc there is spontaneous symmetry breaking 
and two possible solutions with amplitudes ±^o- The am- 
plitude ^0 can be obtained by expansion at higher order of 
53|) and vanishes ifh = ttL. In this case the second of 
gives a finite velocity gradient and the film spon- 



Eq, 



Eq. (53 



taneously flows with a finite flux. Note that in general 
the onset of spontaneous flow is controlled by the sign of 
the combination ({1 -\- ui). The flow coupling coefficient 
(also known as flow alignment parameter) ui can in gen- 



eral have both positive and negative values (de Gennes 



[and Prost[|1993, ). It is controlled by the shape of the ac- 
tive units and the degree of nematic order. Deep in the 
nematic state, vi < —1 corresponds to elongated rod-like 
particles, while vi> 1 corresponds to disk-like particles. 
The onset of spontaneous flow is therefore controlled by 
the interplay of the contractile/tensile nature of the ac- 
tive forces and the shape of the active particles. A de- 
tailed description of this can be found in [Edwards and 
Yeomansl (|2009D andlGiomi et a/. I (120081). 



This transition is very similar to the Frederiks transi- 
tion of nematic liquid crystals in an external electric or 



magnetic field (de Gennes and Prost 1993). The active 
stress CA/i plays here the role of the external field. If 
the thickness is larger that the critical value Lc, a dis- 
tortion of the polarization appears. Any distortion of 
the polarization creates a gradient of active stress that 
must be balanced by a viscous stress which implies the 
appearance of a finite flow. 

The Frederiks transition could also be considered at 
a constant film thickness varying the active stress CA/i. 
The film spontaneously flows if the active stress is large 
enough (in absolute value). 

Finally, to properly describe the spontaneous flow tran- 
sition for polar active films one needs to consider a two- 
fluid model that allows for variations in the concentra- 
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tion of active particles. In this case spontaneous flow is 
accompanied by spatial inhomogeneities in the concen- 



tration or "banding" not seen in active nematics ( Giomi 



et al. 



2008). In addition for stronger values of activity 



in polar films steady spontaneous flow is replaced by os- 
cillatory flow that become increasingly complicated and 
even chaotic for strong polarity. 




FIG. 22 (color online) A film of ordered polar active suspen- 
sion, illustrating the anchoring boundary condition on the 
filaments, parallel to the free surface and the base. 



2. Instabilities of thin films 



The Fredericks active film instability discussed in sec- 



tion IV. A. 1 occurs in situations where the geometry of 
the film is fixed and its surface cannot deform, whereas 



many beautiful phenomena in thin film flow ( Oron et al. 



1997 Sarkar and Sharma, 2010) involve distortions of the 



free surface. In the context of cell biology, a study of a 
film of active fluid with a dynamic free surface is the nat- 
ural starting point for a complete description of a mov- 
ing lamellipodium, which is a thin, flat, fluid projection, 
full of oriented actin, that forms the leading edge of a 
crawling ceh (ISmah et a/.l 120021 1 Verkhovsky et a/.l|l999). 



The spreading of bacterial suspensions ( Bees et al. 



2000 



[2002|) is another biological instance in which the processes 
discussed in this section could intervene. However we dis- 
cuss here simpler situations where mechanisms such as 
cell division in the bacterial case and actin treadmilling 
in the lamellipodium case are not included and that in 
a first step could only be compared to model biomimetic 
experiments. Despite these limitations, the problems of 
an active drop or an active film are interesting as novel 
variants of classic fluid mechanics problems and as set- 
tings for phenomena of relevance to biology. 

We review in this section the hydrodynamics of thin 
films of fluid containing orientable degrees of freedom and 
endowed with locally uniaxial active stresses, bounded on 
one side by a solid surface, and on the other side by a 
surface free to undulate in response to flows in the film. 
We present in some detail the case of an unbounded film 



( Sankararaman and Ramaswamy[|2009[ ), to highlight the 
exotic physical effects that arise from polar orientational 
order in an active system. In this example, we consider an 
active system comprising a solvent fluid and dissolved ac- 
tive particles. In the language of the previous section, we 
thus consider a multicomponent active fluid. We assume 
that the polar material velocity with respect to the back- 
ground fluid (related to the relative current between the 
two components) is strictly slaved to the polar order and 
equals vqP- This assumption is meaningful in the case 
of bacterial colonies but would have to be reconsidered 
in the discussion of a cell lamellipodium. We also sketch 
results for the case of a finite, partially wetting drop with 
small equilibrium contact angle and apolar orientational 



order (Joanny and Ramaswamy, 2012), where topologi- 



cal defects play a role. In both cases, we consider only 



local orientation field is anchored parallel to the free sur- 
face and the rigid base, and free to point in any direction 
in the plane of anchoring. Our focus is of course on the 
effects specifically due to the active stresses and currents. 

Following [Sankararaman and Ramaswamy ( 2009| ), we 
consider a fluid film (Fig. [22| containing active particles 
with number density p(r, t) and orientation field p(r, t) 
as a function of time t and three-dimensional position 
r = {rj_^z)^ where z denotes the coordinate normal to 
the horizontal coordinates = and the solid sur- 

face lies at z = 0. In the case of the drop, where we 
consider only apolar order, we identify the vector p with 
the nematic director field, with p — p symmetry. The 
free surface is located at 2: = h{rj_^t). The flow of the 
film is characterized by the 3-dimensional incompressible 
velocity field v(r,t). We focus on the case where the 
film or drop has macroscopic order, i.e., the mean (p) 
is nonzero, which means that variations in the direction, 
not the magnitude, of p plays the central role. Our aim 
is to obtain an effective equation of motion for the thick- 
ness h and the z-averaged density and polarization. This 
requires solving the Stokes equation in the presence of 
stresses generated by the active particles. At the end of 
the section, we offer a qualitative physical explanation as 
well. 

The dynamics of the height h is related to the veloc- 
ity through the kinematic condition h = Vz — v± ■ ^ ±h 



(Stone, 2005). The incompressibility of the suspension 



implies volume conservation, so that the height dynam- 



ics becomes a local conservation law ( Stone 2005 ) 



dth^V^-{hv^) = 



(55) 



in the ± plane, where is the in-plane velocity field 
averaged over the thickness of the film. With our as- 
sumptions, the flux of active particles in the laboratory 
reference frame is j = p(v + vop) so that the concentra- 
tion p obeys the continuity equation 

dtp = -V • [p(v + ^op)] , (56) 

which simply generalizes ( [2a[ ) to the case where the par- 
ticles are moving not through an inert background but 
through a suspension with velocity field v. The velocity 
field V obeys the Stokes equation 



VP - V • (0-^ + aP) , 



(57) 



planar alignment (de Gennes and Prost 1993) where the 



with viscosity 77, 
—(Apppp and cr^ 



pressure P, and stresses cr^ = 
arising from activity with strength 
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("A/i and nematic elasticity fde Gennes and Prost 1993), 



respectively, as given in Eqs. (45b) and (45a). Note that 
we have imposed here an explicit linear dependence of the 
active stress on the local concentration p of active par- 
ticles as expected in a bacterial suspension whereas for 
filament- motor systems the active stress increases faster 
than linear with the filament density. The polar order 



parameter p obeys (adapting Simha and Ramaswamy 
(2002a) to the case where p is coupled to the free-surface 



tilt) 

^ + ^l^«/3P/3 + Ai(p-VK = + ^ + (^S) 

where is the time-derivative in a frame comoving and 
corotating with the fluid, defined in Eq. (36), vi is the 
flow alignment parameter introduced in Eq. (41b) and 



familiar from liquid-crystal hydrodynamics ( |de Gennes vx(^) = 



and Prost [ [19931), Ai is the coefficient of the advective 
nonlinearity in Eq. ( 2b ) , and is the strain rate tensor 
This is the liquid limit (r 



defined in Eq. ( |35a 
of Eq. 



0) 



(49) where the polar active term proportional 



to Ai has been included. Finally, Fp is the free-energy 
functional for p and p given in Eq. (|4| , whose content will 
be discussed further below. Equation ([58| generalizes Eq. 
([5| to the case where a fluid medium is present. 

It has been argued in [Sankararaman and Ramaswamy] 
( |2009| ) that in a geometry of finite thickness in the z 
direction, symmetry could not rule out a term of the form 
{C/h)^ in the averaged equation of motion for p. 
Such a term arises through the interaction of p with the 
free surface, and the coefficient C encodes the preference 
of p to point uphill or downhill with respect to a tilt of 
the free surface. Possible microscopic mechanisms and 
estimates of magnitudes for such a term are discussed in 



Sankararaman and Ramaswamy ( 2009 ). A key additional 



remark needs to be made: such an effect, while allowed 
by symmetry, must explicitly involve properties of the 
free surface and the base; it cannot emerge simply from a 
averaging of the bulk 3d hydrodynamics. We therefore 



add such a term to ( 58 ) 



We now proceed to solve the Stokes equation (57) for 



the velocity in terms of p, p and /i, restricting ourselves 



to the lubrication approximation |Batchelor| ( 2000 ) and 
Oron et a/.|([T997|) = 0, \d^v\ > |V±v|. 



We study perturbations about a reference configura- 
tion of the film with uniform concentration po and height 
ho, spontaneously ordered into a state with nonzero mean 
polarization (p) = po^k. As we have chosen the active- 



particle current relative to the medium in (56) to be en- 
tirely along p, with no explicit diffusive contribution, 
and as the particles cannot escape the fluid film, the 
normal components of p must vanish at the bounding 
surfaces at 2; = and z = h in agreement with the pla- 
nar alignment condition that we impose. In a perturbed 
state with a non-uniform film thickness, this means that 
Pz{z = h) dxh to linear order in V/i. The z direction 



being the smallest dimension in the problem, it is con- 
sistent to assume that the variation of p with respect 
to z is at mechanical equilibrium via nematic elastic- 
ity. The instantaneous value of is then related to the 
thickness profile by Pz = {z/h)dxh. In the z-averaged 
description that we are aiming for, pz = {l/2)dxh and 
dzPz - h-^dxh. 

Let us represent the perturbed state as p^ = x + 
^ <^ 1 The divergence of the active stress has com- 
ponents dacr^x = C^KdyO ^ dxp/ po ^ h-^dxh) , dacr^y = 
(^ApdxO and dacr^z — to linear order. We ap- 



ply the thin film approximations of (Oron et al. 1997 



Stone 2005 ) to first calculate the pressure and then ob- 



tain the linearized expression of the averaged in-plane 
velocity 



hz - z^l2 
V 



-iV^V\h - -CAphdlV^h - f_i 



. (59) 

where = (Ap[{dyO + dxp/ Po + ^ ^dxh)ii + dxOy] con- 
tains the dominant contributions of activity. Inserting 
this result in the incompressibility condition ( [55| and lin- 
earizing h = hQ-\- 6h, p = Po-\- Sp, leads to the dynamical 
equations 



dtSha 



_CAphl 
(1 - lhlq')qlSh^] 



S p 
Po 



3p 



(60) 



for the in-plane spatial Fourier transforms 5hc^(t), 5pc^{t), 
6>q(t) of the perturbations in height, concentration and 
orientation. The group of terms multiplied by 
the right-hand side of (60) displays four consequences of 
activity, viz., from right to left: (i) pumping by curva- 



ture (see Fig. 23); (ii) anisotropic osmotic flow, through 
the interplay of a uniform orientation field and a con- 
centration inhomogeneity; (iii) splay-induced flow from 
tilting the free surface (the spontaneous flow instability 
discussed in the previous section, but adapted to the case 
of a deformable surface) and (iv) an active anisotropic 
contribution to the effective tension. Following these ac- 
tive terms is conventional surface tension. Terms (iii) 
and (iv) are, respectively, destabilizing and stabilizing 
for contractile stresses, and the other way around for ex- 
tensile stresses; the reader will see that this is physically 
reasonable. We now turn to the equation of motion 
for the polar orientation p. After linearizing (58), and 
averaging over z, in the hydrodynamic limit, we obtain 

iC 

-qyShq^ - {D^ql + D^qy + iXiqx) 0^ 



ho 

- {KQy - ^qxqy)Spci, (61) 

where D± = D - {Xi ± l)hlCAp/4r], D - K/rj being a 
director diffusivity with K the Frank constant (assuming 
that the Frank constants are all equal), and ^ = {ui — 
l)hlCAp/Apor]. 
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FIG. 23 (color online) Curvature of the director of an active 
liquid crystalline phase gives rise to flows. The arrows indicate 
direction of flow if the active stresses are contractile. 
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FIG. 24 (color online) Tilt of the free surface leads to active 
flow; the arrows directions are for the contractile case. 



Linearizing the active particle conservation law (56) 



about the ordered uniform state gives to leading order in 
wave number: 



dt^Pci = -ipoVoQyOd-ivoqx^Pq 



-0{ql6p^,qlSh^). (62) 



The most accessible and interesting instability, arising 
from the combination of activity and the tilt coupling 
C, can be understood by ignoring the concentration and 
motility {vq = 0), but retaining contractile or extensile 
active stresses (CA/i ^ 0). The dispersion relation has 
the complex form: 



UJ = ± 



1 + zsgn(ga;C'CA/i) 
V2 



1/2 



\CCAf,q,\'/'\qy\ 



(63) 

The relative signs of C^P C determine the direc- 
tion (±x) of propagation of the unstable mode. Fig. [24| 
attempts to explain the mechanism of this intriguing in- 
stability. The effects of concentration fluctuations, and 
the suppression of the instability as the motility vq is in- 
creased, are discussed in |Sankararaman and Ramaswamy 
(2009). 

So far we have assumed an unbounded film. Let us now 
consider briefly the spreading kinetics of a finite drop, as 



discussed by Joanny and Ramaswamy (2012), containing 



apolar active filaments in a state of nematic order. We 
choose a planar anchoring condition as in the case of the 






FIG. 25 Left: two half- strength defects (boojums) at the 
ends of a drop. Right: a radial strength- 1 aster defect. Both 
defects are viewed from above the plane in which the drop 
lies. The lines are the director orientation. 



unbounded film, i.e., no component normal to free surface 
or base. The geometry of the drop forces a topological de- 
fect in the interior that itself induces a deformation of the 
drop. The simplest case to consider is a defect consisting 
of two "boojums" with orientation pattern as in Fig. 25 
(left), with the assumption that the drops spreads uni- 
axially in the x direction. Two possible defect structures 
are two "boojums" with orientation pattern as in Fig. 



[25] (left) and an aster as in Fig. 25 (right). The active 
stress amounts to a peculiar kind of disjoining pressure 
Pact = C^P^^{h/ho), which can be rationalized on di- 
mensional grounds by noting that the activity strength 
itself has units of stress so that the dependence on thick- 
ness has to be logarithmic. The resulting static shape of 
a drop depends on the sign of the active stress: the drop 
is flat if the active stress is extensile, and elevated, if it 
is contractile. For a linear structure like Fig. [25] (left), 
if we assume spreading only along the long axis at fixed 
width w the result, for a fixed volume Q and viscosity r^, 
is a drop of linear dimension 



R{t) 



1/4 



(64) 



A drop of volume Q with an aster defect as in Fig. 25 
(right) spreads isotropically with 



V V 



1/6 



(65) 



Further details including the behavior of other defect 
configurations are discussed in [Joanny and Ramaswamy] 
([2012]). 



B. Polar active suspensions with inertia 

In section IIIII we summarized the continuum flock- 



ing model of [Toner and ^ ( [1995| [l998^ for polar 



self-propelling particles moving on a frictional sub- 
strate. Deep in the ordered phase, the theory predicted 
novel propagating sound-like waves, coupling orientation 
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and density fluctuations, with direction-dependent wave 
speeds and giant density fluctuations. The discussion 
thereafter, when the ambient fluid medium was intro- 
duced, focused on destabilizing effects arising from the in- 
terplay of activity and the hydrodynamic interaction for 
bulk suspensions in the Stokesian regime where viscosity 
dominates. However, we claimed in section |T] that our hy- 
drodynamic approach applied to flocks from sub-cellular 
to oceanic scales. Bacterial suspensions, cell aggregates. 



and three waves, generalizations of those in Toner and 



and the cytoskeleton or its extracts (see section III) are 



well approximated in the Stokesian limit where inertia 
is altogether ignored. For collections of large swimmers, 
such as flsh, where inertial effects are important and the 
role of viscosity can be ignored an alternative emphasis is 
appropriate. Simha and Ramaswamy (2002a) originally 



Tu (1998), coupling splay, concentration and drift each 
with exceedingly complicated direction-speed wave rela- 
tions (Simha and Ramaswamy, 2002a). The bend/twist 
waves result from the interplay of V x p± and V x , 
which provides a qualitative difference with respect to 
dry flocks. These propagating modes can be observed 
on length scales where it is reasonable to ignore their 
damping - due to viscosity for the total momentum and 
velocity fluctuations, and diffusion for concentration fluc- 
tuations. This gives a large range of length scales for 
large fast swimmers like flsh. Possibly experiments like 



those of Makris et al (2006), which do speak of flsh den- 



formulated their general theory of ordered states and fluc- 
tuations for active particles suspended in a fluid with a 
view to describing both viscosity- and inertia-dominated 



flows. Let us briefly review the discussion of Simha and 



sity waves, could test the existence of these modes in 
detail. Secondly, within this linearized treatment, equal- 
time density correlations of the density show the same 
features in the ordered phase as we already saw for dry 
systems in section II.A.3.c[ with the number variance, 
scaled by the mean predicted to diverge as A^^/^+^Z^ 



Ramaswamy (2002a) in the case where inertia is taken in d dimensions as in Eq. (20). Presumably nonlineari 



into account through the acceleration term in the mo- 
mentum equation. We ignore the advective nonlinearity 
V • Vv, so the treatment amounts to the unsteady Stokes 
equation, but this is precisely the level of description in 
which the propagating modes of translationally ordered 
crystalline or liquid-crystalline phases are discussed by 



ties, which power-counting readily shows to be relevant 
below 4 dimensions, but which are even more painful to 
analyze here than in dry flocks, will change the exponents 
but not the essential fact of supernormal fluctuations. It 
is not clear that the power-law static structure factor seen 



Martin et al (1972) 



by Makris et al. ( 2006 ) in shoals is evidence for these gi- 



The slow variables in a polar active suspension are 
the number density, the broken symmetry variable p for 
polar order (deflned in section-I A. ) and the momen- 
tum density of the suspension g coming from momentum 



ant fluctuations; a fair test of flocking theories requires a 
school rather than a shoal. 



C. Rheology 



conservation. We proceed initially as in section III A. 2, 
through the momentum equation (31) dtga = —df3lla(3^ 
with the stress tensor Ilai3 



The active hydrodynamic framework of section III 



not 



^P'^q;'^/3, where only allows us to predict spontaneous-flow instabilities as 



(j^^ contains the various contributions listed in Eqs. (38), 



(41a), (41b) and (44), including the crucial active stress. 



in sections [IV. A. 1| and |IV. A.2] but also the response of an 
active fluid to an imposed flow, i.e., the rheology of active 



In systems where inertia is important, it is essential to re- 
tain the acceleration dtg even while ignoring the inertial 
contributions mpVaVfs to the stress which are non-linear 
in g. The equation for the vector order parameter p is 
as in section III The number density is governed by Eq. 



soft matter (Giomi et 



tillan 



2010). 



2010 



Haines et al, 2009; Hat- 



walne et al. \ |2004[ [Liverpool and Marchetti^ ,2006^ ,Sain- 



For concreteness, we associate the purely coarse- 



(56). 



The hydrodynamic modes implied by these equations 
of motion are obtained by linearizing and Fourier trans- 



forming in space and time (Simha and Ramaswamy] 



2002a). Imposing overall incompressibilty with the con- 



dition V • g = 0, the number of modes is flve. We briefly 



state here the main results of this analysis (Simha and 
Ramaswamy, 2002a). First, when viscosity is ignored 



grained description of section III with a microscopic pic- 
ture of a suspension containing active particles of linear 
size at concentration p, each particle carrying a force 
dipole of strength where / is the propulsive thrust, 
with the activity of an individual particle correlated over 
a time tq, and collective fluctuations in the activity cor- 
related over length scales ^ and timescales r. 

We apply the approach here to the isotropic phase 
of active particles in a fluid, to extract linear rheolog- 



and acceleration is included, polar ordered suspensions 
are not in general linearly unstable, that is, a parameter 
range of nonzero measure exists for which stable behav- 
ior is found. The dynamic response displays a whole 
new range of propagating waves as a result of the in- 
terplay of hydrodynamic flow with fluctuations in ori- 
entation and concentration: a pair of bend-twist waves 



ical properties (Hatwalne et a/., 2004 Liverpool and 



Marchetti , 2006 ) and the autocorrelation of spontaneous 



stress fluctuations (Chen et al. 2007 Hatwalne et al. 



2004) when noise is included. 

We consider an apolar system described by orientation 



order in terms of the tensor Q, deflned in Eq. (21).. We 



obtain predictions for the rheology of active suspensions 



32 



through the coupled dynamics of Q and the hydrody- 
namic velocity field v = ^/ ptot^ where g and ptot are 
the total densities of momentum and mass of the sus- 
pension. Neglecting inertial contribution to the stress 
tensor, conservation of total momentum of particles + 
fiuid is expressed by dt^ — V • cr^, where the stress ten- 

In addition, we must allow 



sor IS as m section 



III.B.l 



for noise sources of thermal and non-thermal origin. The 
former, though mandated in the equilibrium limit by the 
fiuctuation-dissipation relation connecting them to the 



solvent viscosity in (38), are quantitatively negligible in 



comparison to effects arising from activity. Such active 
fiuctuating stresses can further arise in two ways: di- 
rectly, as additive white-noise contributions to the stress 



in (38)- (44), which we shall ignore as they add no new 



physics, or indirectly, via the stochastic dynamics of the 
orientation field, which is the case of interest. 

Consistent with Section [TITl we write the active stress 
in terms of the alignment tensor in the form 



-CA/iQ 



a/3- 



(66) 



We note for later reference that we will eventually con- 
sider the concentration dependence of active stresses, 
C (X p, which is physically reasonable and also follows 
from an explicit realization in terms of dipolar force den- 



sities associated with the active particles (Baskaran and 



Marchetti 



2009 



Hatwalne et al\ |2004| |Simha and Ra- 



maswamy[ |2002a[ ). We remind the reader that we are 
dealing with force-free, neutrally buoyant, self-propelling 
particles. Thus there are no external forces on the 
system, so that the simplest active particle, on long 



timescales, is a permanent force dipole (Brennen and 



Winet, 1977 Pedley and Kessler, 1992). Although we 



have already discussed the magnitude and sign of ^A/i 
in section it is useful to note here that for a system 
with concentration p the quantity C^/^pj p ^ if character- 
izes the strength of the elementary force dipoles associ- 
ated with, for example, individual swimming organisms. 
Negative and positive C^P refer respectively to contrac- 
tile swimmers, or "pullers", and extensile swimmers, or 
"pushers" , whose distinct rheological behavior we outline 
below. 

The rheology is to be obtained from the momentum 
equation together with the equations of motion for the or- 
der parameter field and particle concentration. We begin 
with a description of the linear rheology of active matter 
in the isotropic and orientationally ordered phases, and 
follow it up with a brief survey of the nonlinear rheology 
of active matter. 



1. Linear rheology of active isotropic matter 

To appreciate what is unique about active matter rhe- 
ology, it is useful to recall the linear rheology of pas- 
sive nematogens. The stresses arising from distortions 



(a) 



(b) 



(c) 



(d) 



# 



FIG. 26 (color online) Discs (a) and (c), and rods (b) and 
(d) with active force densities attached along their symmetry 
axes, under shear (horizontal arrows). The parameter C^P < 
in (a) and (b) and > in (c) and (d). 



of the orientational order parameter Q and concentra- 
tion p are derived from a free-energy functional Fq[Q, p] 
(Eq. (23)), giving rise to a passive deviatoric order- 



parameter stress (Forster 1974) 



3H - H Q - Q H 



(67) 



where H = -SFq/SQ + (1/3)1 TrJFg/JQ is the nematic 
molecular field. The mean deviatoric passive stress (67) is 



zero in the isotropic phase at equilibrium (and in the ne- 
matic phase as well). In addition, fiuctuations of the de- 
viatoric stress are small as one nears the transition to the 
nematic phase, as can be seen by making small perturba- 
tions SQ in the isotropic phase, leading to a change in the 
free energy Fq oc {aQ/2) J Tr((5Q)^, which in turn give 
rise to stress fiuctuations ~ ag SQ with a coefficient ag 
that decreases on approaching the isotropic-nematic (IN) 
transition to the ordered phase. Thus even though fiuctu- 
ations of Q are large as one approaches the IN transition, 
their contribution to rheology is small, resulting simply 
in a renormalisation of r, the order parameter relaxation 
time. As shown below, this pretransitional feature is fun- 
damentally different in active isotropic systems. In the 
following we take a purely coarse-grained approach and 



follow closely the work of Giomi et al. (2010); Hatwalne 
et al. (2004); and Liverpool and Marchetti (2006); for a 



more microscopic treatment see 
Saintillanl (|2010l). 



Haines et 



(2009) and 



For the active system, the deviatoric reactive stress 
has to be obtained from the equations of motion for the 
order parameter and concentration, rather than simply 
from the free energy functional. The linearized equations 
for Q in the isotropic phase, are of the form. 



dt 



1 



■fa^. (68) 
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where Eq. (68) can be regarded as the modification of 



Eq, (22) to include the effect of shear flow to linear order. 
Here r is the orientational relaxation time, which could 
for example be the run time in a collection of run-and- 
tumble bacteria, or the rotational diffusion time, per- 
haps modified by collective and/or active effects in an 
actomyosin extract, I) is a diffusivity (~ ^^/'Tq) related 
to the ratio of a Frank constant to a viscosity, ui is a 
"reversible" kinetic coefficient or flow coupling parame- 
ter, taken for simplicity to be of the same order in the 



one entering the equation for the polarization p ( Forster 



1974), fa(3 is a traceless, symmetric, spatiotemporally 
white tensor noise representing active fluctuations, and 
the dots include the coupling of orientation to flow. Note 
that this form of the linearized equation is valid for the 
passive system too, with the time scale r given by the 
order parameter relaxation time proportional to 
which gets larger as one approaches the IN transition. 

One can now proceed to calculate the linear viscoelas- 
tic properties of the active suspension. In addition to the 



active deviatoric stress (66), we include the contribution 
from the viscous dissipative stress given in Eqs. (38). 



For simplicity we neglect the tensorial nature of the vis- 
cosity of fluids with orientational order. Using (|66|, ( [67[ ) 
and (68), and applying them to spatially uniform {q = 0) 
oscillatory shear flow at frequency uj in the xy plane one 
obtains, to linear order in the fields. 



(69) 



The rheological response is defined by the complex 
modulus G{uj) = with exy{uj) = Vxy{uj) / {—iu) 

the strain. The corresponding storage (in-phase) and 
loss (out-of-phase) moduli G' {uj) and G"{uj)^ defined by 
G{uj) = G'{uj) ^iG"{uj)^ characterize the elastic and vis- 
cous response of the system to an oscillatory shear flow. 
These moduli can be read out from Eq. (69). 



The active isotropic system is rheologically a Maxwell 
fluid to linear order. This is best highlighted by 
the behavior of the apparent shear viscosity rjapp = 
\Ym^^QG{uj) / {—iuj)^ which shows an active excess vis- 
cosity rjact oc — CA/irp, corresponding to either an en- 
hancement or reduction, depending on the sign of (. This 
active thickening (thinning) can be understood as follows 
(Fig. |26|: in an imposed flow, in the absence of activity, 
discs (rods) tend to align their symmetry axis along the 



compression (extension) axis of the flow (Forster, 1974). 
When activity is switched on, the flow induced by the 
intrinsic force dipoles clearly opposes the imposed flow 
in Fig. 26 (a) and (b), and enhances it in (c) and (d). 



Activity thus enhances viscosity in Fig. [26] (a) and (b), 
since (A/j. < 0, and reduces it in (c) and (d) (CA/i > 0). 
For ("A/i < (69) shows that the viscosity grows substan- 



tially as the system approaches a transition to orienta- 
tional order as r is increased, and in fact should diverge 
if r could grow without bound. Even more strikingly, a 



swimming cells 



dead cells 



H8H 
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FIG. 27 The contractile activity of "puller" swimmers can 
increase the viscosity of a suspension. The data represent the 
effective viscosity of chlamydomonas suspensions relative to 
the viscosity 770 of the culture medium as a function of the 
volume fraction of bacteria. Solid symbols represent live cell 
data and crossed symbols represent dead cell data. The mech- 
anism producing flow-orientation, however, remains unclear. 
Adapted with permission from (Rafai et al. 2010). 



system of pushers, i.e., extensile swimmers, should show 
a prodigious reduction in viscosity as r grows; indeed, 
nothing rules out a negative viscosity, that is, sponta- 
neous flow in an initially quiescent isotropic system. Ex- 



periments by Sokolov and coworkers (Sokolov and Aran- 
2009) have indeed shown that the extensile activ- 



son 



ity of Bacillus subtilis, a "pusher" swimmer, can substan- 
tially lower the viscosity of a suspension. By contrast, in 
a passive, i.e., thermal equilibrium, system approaching 
a continuous or weak first-order transition to a nematic 
phase, the excess viscosity ~ aqr is roughly constant 
since r ex l/ag, as required by the constraints of ther- 
mal equilibrium. 

The active excess viscosity obtained within the linear 
theory should be compared to the well-known result of 
Einstein] ( |1906[ |1911} that the fractional excess viscosity 
due to the addition of passive particles to a fluid is pro- 
portional to the particle volume fraction (j) = npi^/G^ to 
lowest order in 0, with a coefficient 5/2 for spheres. For 
this purpose it is convenient to take the active coupling 
proportional to the concentration p and define an active 
stress per particle W as 



CA/x = fep = w<t> 

where W = f/i"^. Then 

Vact = -Wt4>. 



(70) 



(71) 



Equation ( 71 ) can be viewed as an additive correction to 



the 5/2, proportional to VKr, which of course can be of 
either sign. Behavior consistent with these predictions 
is seen in recent experiments measuring the activity- 
induced thickening in a system of Chlamydomonas algae 



(pullers, W < 0), as shown in Fig. IV.C.l ( (Rafai et al. 
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2010), and extreme thinning in a system of Bacillus sub- 
tilis bacteria (pushers, W > 0) (iSokolov and Aranson 



2009 ). In the case of chlamydomonas, there is some ques- 
tion as to the direct apphcabihty of the mechanism dis- 
cussed above, as recent experiments have shown that the 
flow field generated by these microorganisms, although 



contractile, is more complex than dipolar (Drescher et al. 



2010 Guasto et a/., 2010). In general the linear rheology 



is controlled by the interplay of the nature of the active 
stresses, determined by and the flow alignment coeffi- 
cient ui. A remarkable duality has been identified that 
shows that shows that tensile {( > 0) rod-shaped flow- 
aligning particles > 1) are rheologically equivalent 
(to linear order in the strain rate) to contractile {( < 0) 
discotic flow-tumbling particles > 1) (Giomi et al. 



2010) 



Equation (|69|) also predicts strong viscoelasticity as r 
increases. For passive systems, (A/j. = and aq oc r~^, 
and so G'{uot ^ 1) decreases as i^i/^/r, suggesting as ear- 
lier, that there is little viscoelasticity near an equilibrium 
IN transition. For active contractile (CA/i < 0) systems, 
by contrast, the active contribution is of 0(1) and inde- 
pendent of r even close to the IN transition. Thus, as r 
grows. 



G'(a;r > 1) - -(A/i 



(72) 



independent of r! In addition to this enhanced elastic- 
ity, the dynamic range over which elastic behavior is seen 
increases. This is quite a dramatic rheological manifes- 
tation of activity, since at equilibrium^ one would expect 
such strong viscoelastic behavior from a fluid or suspen- 
sion near translational freezing, as at a glass transition, 
not near orientational ordering. Put more radically, the 
orient at ionally ordered phase of contractile active par- 
ticles is a peculiar yield-stress material, with nonzero 
shear and normal stresses in the limit of zero shear rate 



Marchetti 



m section 



(Giomi et a/., 2010 Hatwalne et a/., 2004 Liverpool and 



2006 Marenduzzo et 



2007) as discussed 



Another observable manifestation of activity is an en- 
hanced noise temperature as inferred for example from 
tracer diffusion measurements in bacterial suspensions 
( |Wu and Libchaber[ j2000|, and frequency-dependent 
shear viscosity arising from fluctuations in stress and con- 



centration (Chen et al., 2007 Hatwalne et al., 2004). Ac- 



tivity is the transduction of chemical energy, say in the 
form of ATP hydrolysis equal to about 20k bT per ATP 
molecule. This suggests that fluctuations from this ather- 
mal noise would be significant. Following [Hatwalne et al.\ 
(12004), we may estimate the strength of nonequilibrium 



stress fluctuations that result from fluctuations of the 
force generation of active particles, through the variance 
of the active stress at zero wavenumber and frequency. 
The divergence of the active stress cr^ = — ^A/iQ ap- 
pears as a forcing term in the momentum equation. As- 
suming the active objects, whether biofilaments or bac- 



teria, are collectively in a spatially isotropic state, this 
forcing can be viewed as a noise on scales large com- 
pared to the correlation length and time ^ and r de- 
fined at the start of this section. Recall that in fluc- 



tuating hydrodynamics (Landau and Lifshitz 1998), for 



a thermal equilibrium fluid with shear viscosity 77 and 
temperature T, the variance of the random stresses at 
zero frequency is rjT. For an active fluid, then, the 
apparent temperature as probed by the motion of a 
tracer particle can be estimated by the zero-frequency, 
zero- wavenumber variance of the active stress rjT^ff ~ 
(CA/i)2/d^rdt(Q(0,0):Q(r,t)). The simplest dimen- 
sional argument would then give Tg// ~ (CA/i)^^^r/7^. 
An Ornstein-Zernike form (i/r) exp(— r/^) for the equal- 
time correlations of Q, which requires the introduction 
of a microscopic length which we take to be the active- 
particle size ^, gives instead Tg// ~ (CA/i)^^^^r/7^. From 
([to]), using rough estimates for the bacterial system of 
|Wu and Libchaber| ( |2000| ) / - v/r]£ with r] 10"^ Poise, 
£ ^ 1 jam, (/) ^ 0.1 ^ 20 jam corresponding to a "run" 
for a time r = Is and speed v 20/im/s yields Tg// 
about 400 times room temperature, consistent with the 
measurements in ( Wu and Libchaber{ |2000| ). Note that 
this enhancement is independent of the sign o/^A/i, i.e., 
the diffusivity increases regardless of whether the viscos- 
ity decreases or increases, a tell-tale sign of the nonequi- 
librium nature of the system. A further interesting con- 



sequence (Hatwalne et al. 2004) of this excess noise is 



a huge enhancement of the amplitude of the well-known 
f-d/2 long-time tails in the autocorrelation of tagged- 
particle velocities. We know of no experiment that has 
probed this last feature. 

More extended analyses of fluctuations in active sys- 
tems, demarcating the dynamical regimes lying within 
and beyond the conventional fluctuation-dissipation the- 
orem include work by Kikuchi et al. (2009) and Mizuno 
(I2OO7I). 



et 



We now turn to nonlinear fluctuation effects (Chen 



et al. , 2007 



Hatwalne et 



2004). From (66), the fluc- 



tuations in the deviatoric stress gets contributions from 
fluctuations bilinear in the orientation Q and the concen- 
tration p. The stress autocorrelation is therefore a con- 
volution of Q and p correlations. The former is evaluated 



from (68), while the latter can be calculated from the the 
linearized equations for the concentration, dtSp = — V- J, 
where the current Ja = —DdaSp — W p^disQ^iB + Z^? 
where D is the diffusion constant, po the mean concentra- 
tion, a random noise, and W is an activity parameter. 
The resulting key flnd ings of jChen et al.\ ( |2007| ) and Lau 
and LubenskyI \200^ \ are an excess fluctuation (buj ' 
in the stress fluctuations, and no excess response, i.e., 
viscosity, in the range studied, again a sign of nonequi- 



librium behaviour. [Hatwalne et al.\ (|2004j) consider stress 
contributions nonlinear in Q and show they should lead 
to excess viscosity as well of similar form, but the effect 



is presumably below detectable levels in the Chen et al. 
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(2007) experiment. 



2. Linear rheology of active oriented matter 

In extending the study of rheology to active oriented 
matter, we immediately encounter a problem. As re- 
marked in Section IV.A.1[ long-range uniaxial orient a- 
tional order, whether polar or apolar, in active Stokesian 
suspensions of polar particles is always hydrodynamically 
unstable to the growth of long wavelength splay or bend 
fluctuations, depending on the sign of C^M- This insta- 
bility has no threshold in a spatially unbounded system 



(Simha and Ramaswamy 2002a), and the growth rate 
is highest at wavenumber g = 0. Frank elasticity with 
stiffness K stabilizes modes with q greater than 



go (X VlCA/i|/i^, 



(73) 



so that there is a band of unstable modes from to go- 
Note that go is simply an approximate form of the length 



scale defined by Eq. (54). In any case, in the limit of infi- 



nite system size there is no stable reference state, no ideal 
active nematic or polar liquid crystal whose rheology one 
can study as a geometry-independent material property. 
We must therefore ask what suppresses this generic in- 



stability ( Ramaswamy and Rao 2007 ) . A mechanism for 
suppression of the instability is confinement by bound- 



ary walls. The existence of a crossover wavenumber (73) 



implies that the instability exists only if the sample's nar- 
rowest dimension h > ^/K/\(Ajii\] equivalently, for fixed 
h the activity must cross a threshold ^ K/h'^. This is 
the essential content of the treatment of lVoituriez et aL\ 



(2005) discussed in section IV. A. 1 Confinement along 



with the director spontaneously aligned along x and free 



to turn in an u nbounded xz plane, shows (Ramaswamy 
and Rao , 2007 ) a similar threshold but with a g^ depen- 



dence of the growth rate at small in-plane g. In either 
case, it is clear that confinement can produce a stable 
active liquid crystal whose linear rheology one can study. 
Alternatively, the instability can be suppressed by impos- 
ing an external shear flow or by the presence of partial 
translational order, either columnar or lamellar. Orien- 
tational stabilization can also be achieved in flow align- 
ing systems > 1). by imposing a uniform shear 
flow with shear rate 7. This stabilizes those wave vec- 
tors whose growth rate is smaller than the shear rate 7. 
As 7 is increased more and more modes are stabilized, 
until at a threshold 7c, the oriented phase is completely 
stabilized by the shear flow, yielding a stability diagram 
controlled by two variables, the flow alignment parame- 



ter z/i and the ratio of the shear to active stress (Giomi 



et al 2010 Muhuri et a/., 2007). Translational order 



both partial as in smectic or columnar liquid crystal, or 
full as in three-dimensional crystalline systems can also 
yield a stable active system, whose properties have been 
the subject of recent studies ( Adhyapak et a/.^|2012| ). 



Stabilizing the orientational phase of active matter 
now sets the stage for a study of its unusual rheology 



and Marchetti 



(Giomi et al. 2010 Hatwalne et al., 2004 Liverpool 



2006). We have already seen that on 



approaching the orientationally ordered state from the 
isotropic fluid, a suspension of active contractile elements 
{W < {)) exhibits solid-like behaviour without transla- 
tional arrest. In the orientationally ordered phase, the 



orientational order parameter (Q) 7^ 0, which by (66) im- 



mediately leads to a nonzero steady-state average of the 
de viator ic stress, in the absence of any external deforma- 
tion. This prestress is a truly nonequilibrium effect, and 
has no analogue in a passive equilibrium nematic fluid 
which has a purely isotropic mean stress, i.e., a pressure, 
despite its orientational order. This prestress implies that 
in a flow experiment, the shear stress will not vanish at 
zero deformation rate. The same features of being first 
order in shear rate and having a nonzero value at zero 
shear rate are exhibited by the normal stresses ayy — cFxx- 
We emphasize that this strange kind of yield stress is a 
manifestation of the rheology of an active oriented fluid 
without any form of translational arrest. 

We end this section on the rheological properties of ac- 
tive oriented matter by drawing the attention of readers 
to an interesting analogy between active stress of con- 
tractile filaments with macroscopic orientational order 
and fragile jammed granular matter, as discussed by |Ra- 
maswamy and Rao (2007). 



3. Nonlinear rheology of active nematics 

We now turn to the rheological properties of orien- 
tationally ordered active fluids beyond the linear regime. 
As detailed in the previous section, the generic instability 



of orientationally ordered active suspensions ( Simha and 



Ramaswamy, 2002a) can be suppressed by confinement 



or by imposing an external shear ( jMarenduzzo et al. 



2007 Muhuri et a/., 2007). It is therefore meaningful 



to explore the dynamics and rheology of these phases 
in confined geometries. Extensive numerical studies of 
active nematic and polar films ( Gates et a/. [ [2008 Field- 
ing et al\ |2011[ [Giomi et al.\ |2010| [Marenduzzo et al. 
2007) under shear reveal a rich variety of phenomena 
influenced by boundary conditions and geometry. This 
complex behavior results again from the interplay be- 
tween local stresses generated by activity (quantified by 
the parameter ("), the flow aligning property of these par- 
ticles (characterized by the parameter vi) and the typ- 
ical self-propulsion velocity ^'o in polar active fluids. A 
complete understanding of their response to shear neces- 
sitates exploring the space spanned by these parameters 
under various boundary conditions. 

We present only a summary of the main results ob- 
tained in the literature. While detailed theoretical and 
numerical studies of linear active rheology including the 
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effect of polarity can be found in (Giomi et al, 2010 



2008), the nonlinear rheology has so far been studied 
mainly for apolar systems. For contractile active fluids 
in an orient at ionally ordered state, numerical studies of 



the hydrodynamic equations (Gates et al. 2008) when 



only one-dimensional variation is allowed show the on- 



set of solid-like behavior of (Hatwalne et al, 2004) and 



the existence of a yield shear stress as in ( Liverpool and 



Marchetti, 2006 Ramaswamy and Rao 



2007). Stud- 



ies of extensile fluids, allowing variation in one (Gates 



et al. 2008) and two (Fielding et al., 2011) spatial di 



rections display the onset of spontaneous flowing states, 
with complex flow states including bands in Id and rolls 
and turbulence in 2d. Related experimental findings in- 
clude large scale chaotic flows in bacterial systems ( Aran- 
Dombrowski et al.\ [2004 ) , as remarked 



son et al. 



2007 



earlier. An important general feature (Gates et al. 2008 



Fielding et al. 2011) is that dimensionality matters: im- 
posed restrictions to Id spatial variation lead to signifi- 
cantly different rheology from that seen when 2d varia- 
tion is allowed, for example. Active stresses generally ap- 
pear to stabilize shear bands, but have the opposite effect 
close to the isotropic-nematic transition. In systems with 
free boundary conditions, the approach to the isotropic- 



nematic contractile transition shows ( Gates et al. , 2008 ) 



an active enhancement of viscosity. As claimed in (Hat 



walne et al. 2004), orientational ordering in active sys- 
tems indeed resembles translational arrest in equilibrium 
systems. 

The relation between the rheology of active fluids in 
external shear and the onset of spontaneous flow in ab- 



sence of shear is discussed in (Giomi et al., 2010). This 
work has also analyzed in detail the nonlinear rheology 
of active fluids, revealing strongly non-monotonic stress 
vs strain-rate curves beyond a threshold value of activity. 



as suggested in Hatwalne et al. ( 2004[ ) , with macroscopic 
yield-stress behavior, and hysteresis. Finally, at even 
higher activity, the theoretical stress-strain curve has a 
discontinuous jump at zero strain rate, corresponding to 
a finite "spontaneous stress" in the absence of applied 
shear. 

Given the complexity and richness of the nonlinear rhe- 
ology of active fluids, we refer the reader to the literature 
for further details. 



D. Applying the hydrodynamic theory to phenomena in 
living cells 

To demonstrate that the hydrodynamic theory of ac- 
tive gels described in the previous sections indeed pro- 
vides generic tools for addressing questions relevant to 
living cells, we briefly summarize here a few examples of 
concrete successes of the theory. A more complete re- 



view can be found in Grill (2011) and Joanny and Frost 




FIG. 28 Velocity field determined by speckle microscopy in 
a keratocyte lamellipo dium. Adapted with permission from 
Vallotton et al\\2mb\. 



The first example is directly relevant to cell motility, 
in particular to the migration of fish keratocytes, eukary- 
otic cells extracted from fish scales that are among the 
fastest moving cells, migrating on a substrate at a steady 
speed of about lOjam/min. This cells have a characteris- 
tic fan shape, with a large fiat region extending in front 
of the nuclei, towards the direction of forward motion, 
known as the lamellipodium. This region is filled with 
a cross-linked actin gel, where myosin motor complexes 
use the energy from ATP hydrolysis to grab on neighbor- 
ing actin filaments and exert stress. Actin polymeriza- 
tion takes place at the leading edge of the lamellipodium, 
while the filaments disassemble in the rear region, in a 
process that plays a crucial role in driving the motil- 
ity. Using active gel theory to model the lamellipodium. 



Kruse and collaborators (Kruse et al. 2006b) were able 



to evaluate quantitatively the profile of the retrograde 
actin flow that accompanies the forward motion of the 
lamellipodium. Such a retrograde flow has been seen in 
experiments, as shown in Fig. 28 (Vallotton et al, 2005). 

A second example is that of shape oscillations, ob- 
served ubiquitously in many cells. An example is shown 
in Fig. [29j This phenomenon has been described success- 
fully by the active gel theory, augmented by assuming a 
coupling to calcium channels in the cell membrane that 
in turn are gated by the deformations in the actin corti- 
cal layer. Galcium concentration couples to local myosin 
activity that in turn controls the stretching and compres- 
sion of the actin layer, in a feedback loop that results in 
sustained oscillations ( [Salbreux et al.\ |2007[ ) . 

The active gel theory has also been used recently to 
model cortical flow in the C. elegans zygote (Mayer 
as shown in Fig 



et 



2010) 



30 The theory supported 



( |2010| ). 



experiments in successfully identifying two prerequisites 
for large-scale cortical flow necessary to initiate the an- 
teroposterior cell polarization which directs the asymme- 
try of the first mitotic division: a gradient in actomyosin 
contractility to drive flow, and a sufficiently large viscos- 
ity of the cortex to allow flow to be long-ranged. 

Another recent success of active hydrodynamics has 
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FIG. 29 Shape oscillations of no n- adhering fibroblasts. The 
second frame shows the periodic oscillation of the projected 
area of the cell and the associated Fourier spectrum. The 
period of the oscillation of these cells is very well defined and 
of the order of 30s. The oscillation period is found to decrease 
when myosin activity increases. The latter can in turn be 
modulated by t he addition of various drugs. Adapted with 
permission from Salbreux et al. (2007). 





FIG. 31 (color online) The top images (A) show a a circu- 
lar actin wave in the lamellipodium in a Drosophila cell fixed 
to a substrate and treated to enhance polarization of actin 
filaments. The wave travels along the cell's periphery, with 
the arrow marking the region of maximum actin density at 
various times. (B) Phase contrast microscopy picture of a 
control Drosophila cell fixed to a substrate with a typical cir- 
cular shape. (C) In the theoretical model, the cell perimeter 
is represented by a circle of radius R. (D) The calculated con- 
centration of active filaments as a function of the polar angle cj) 
and time. Adapted with permission from ,Asano et a/.| ( |2009| . 




% A 



FIG. 30 (color online) This figure displays the imaging of ten- 
sion in the cortex of a C. elegans zygote obtained before and 
after cortical laser ablation (COLA) performed on the acto- 
myosin meshwork. The left image shows a schematic of COLA 
performed with a pulsed ultraviolet laser along a 6/im line 
(light blue in both images) along the anterior-posterior (AP) 
axis, in the posterior of the cell. The center frames show the 
pre- (top) and post-cut image (bottom) of posterior COLA. 
The image to the right shows the tension measurement in the 
zygote's actomyosin cortex upon laser ablation along the blue 
line. The black arrows are the displacements between pre-cut 
(purple) and post-cut images (green). The average speed of 
this initial recoil measured in a direction orthogonal to the cut 
line is proportional to the normal component of the 2D ten- 
sion tensor. The w hite b ar is 4/im. Adapted with permission 
from [Mayer et a/.| ( [20Tot . 



come from its application to the organization and dy- 
namics of molecules on the surface of living metazoan 
cells ( Gowrishankar and Rao, 2012). The lateral compo- 



sitional heterogeneity of the plasma membrane at sub- 
micron scales, termed "lipid rafts", has been the subject 



of intense research (Lingwood and Simons, 2010). Most 



attempts to understand this heterogeneity are based on 
equilibrium thermodynamics. Using a variety of fluo- 
rescence microscopy techniques, and confirmed by other 
methods, it has been shown that the dynamics, spatial 
distribution and statistics of clustering of a key 'raft' 



component, viz., GPI-anchored proteins, is regulated by 



the active dynamics of cortical actin filaments ( Goswami 



et al. 2008). Following this, Gowrishankar and collab- 



orators (Gowrishankar and Rao, 2012) have developed 



a model for an active composite cell membrane based 
on active hydrodynamics, which shows how the dynam- 
ics of transient aster-like regions formed by active po- 
lar filaments (actin) can drive passive advective scalars 
(molecules such as GPI-anchored proteins) to form dy- 
namic clusters. In addition to successfully explaining the 
many striking features of the dynamics and statistics of 
clustering of GPI-anchored proteins on the cell surface, 
this model makes several predictions, including the exis- 
tence of giant number fluctuations which have been ver- 



ified using fluorescence microscopy (Gowrishankar and 



Rao, 2012). These studies suggest that rafts are actively 



constructed and that the active mechanics of the corti- 
cal cytoskeleton regulates local composition at the cell 
surface via active currents and stresses. 



Finally, Fig. [3T] shows polarized actin waves propagat- 
ing around the periphery of a Drosophila cell that has 
been fixed to a substrate, preventing it from moving, un- 



der condition of enhanced actin polymerization (Asano 



et al. , 2009 ) . In this case the hydrodynamic theory of a 



polar active gel successfully reproduces the existence of 
a critical value of polarity above which the polarization 
waves will occur, as observed in experiments. 
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V. MICROSCOPIC DESCRIPTIONS OF ACTIVE 
MATTER 

A. Review of microscopic models 

While the macroscopic equations of motion of active 
matter can be obtained from general considerations of 
symmetry, determination of the magnitude and often 
even the sign of the coefficients requires additional phys- 
ical assumptions. Unlike for systems near equilibrium, 
the assumptions are not easily identified. Furthermore 
the relaxation of the constraints required by equilib- 
rium allows the possibility of a much larger number of 
terms. Physical insight can be provided by using the 
tools and methods of nonequilibrium statistical mechan- 



ics (Zwanzig 2001) to derive the continuum equations 
via systematic coarse-graining of simplified microscopic 
models of the active processes driving the system d Aran- 1 



son and Tsimring, 2005; Baskaran and Marchetti, 2008a| 


Bertin et al 


1 ,2006j ,2009j (Kraikivski et a/.j ,2006, ,Krusel 


et al] 20011 


IKruse and Juhcher| 2000 Liverpool 2003| 


Liverpool et al] 2001). One approach is to start from 



stochastic equations for the microscopic dynamics and 
then systematically project the microscopic degrees of 
freedom on to macroscopic variables, such as the density 
p(r, t) and polarization p(r,t), defined in Eq.( |la|lb ). 



This procedure yields hydrodynamic equations for the 
macroscopic variables on length scales long compared to 
the size i of the individual active elements. For simplic- 
ity, we restrict ourselves here to over damped systems, 
where the medium through which the active particles 
may move (fluid or substrate) is inert and only provides 
friction. In this case the momentum of the active parti- 
cles is not conserved and the microscopic degrees of free- 
dom are the positions r^(t) and orientations i>n{t) of the 
active elements, plus (possibly) other microscopic degrees 
of freedom describing internal active processes. Of course 
medium-mediated hydrodynamic interactions can play a 
crucial role in controlling the dynamics of collections of 
swimmers. To incorporate such effects one needs to con- 
sider a two-component system and explicitly describe the 
exchange of momentum between particles and solvent. 
Both the translational and rotational velocities of each 
particle must be incorporated in the microscopic model, 
as well as the solvent degrees of freedom. This more gen- 
eral case is outlined briefly in section [V.B.4[ In general 
approximations are required both to identify a tractable 
microscopic model and to carry out the coarse-graining 
procedure to derive equations for the macroscopic vari- 
ables. The form of the resulting equations is general, as it 
is dictated by symmetry, and the microscopic description 
allows a calculation of the phenomenological parameters 
and transport coefficients involved in the hydrodynamic 
equations. This is important for example for biological 
systems because it allows to predict the variation of the 
transport coefficients with the biological parameters that 



can be changed in the experiments. 

The purpose of this section is to review a number of 
such derivations and summarize their similarities and dif- 
ferences. In equilibrium statistical mechanics there is 
a long tradition of simplified microscopic models which 
have eventually led to a number of paradigmatic minimal 
models (Ising, XY, Heisenberg) that capture the essen- 
tials of the behavior of variety of equilibrium systems. 
Similarly here we describe some minimal microscopic de- 
scriptions of active matter that have provided insight into 
the complex physics of these systems. An interesting as- 
pect of these microscopic realizations is that they allow 
one to identify the common behavior of different exper- 
imental systems. Two examples which we will focus on 
are (1) mixtures of cross-linking motor constructs and 
protein filaments (|Nedelec[ |1998| |Nedelec eTa/.l |T997| 



Surrey et a/., 2001) and (2) collections of self-propelled 



particles. Experimental realizations of the latter may be 



actin filaments in motility assays (Butt et al. 2010 



and Spudich, 1986 Schaller et al. 2010) or suspensions 



of swimming micro-organisms ( Dombrowski et al. , 2004 



Zhang et al. 2010). 



Kron 



1. Self-propelled particles 

A microscopic realization of active systems is provided 
by interacting self-propelled particles (SPP). While these 
can be thought of as simplified agent based models of 



flocks of birds and shoals of fish ( Vicsek et al. , 1995 ), they 
can also act as realistic models for less complex systems, 
such as polar protein filaments (e.g. F-actin) in glid- 
ing motility assays on surfaces decorated with molecular 
motors (e.g. myosin) ( |Kron and Spudich , 1986 Schaller 



et al. 



2010), suspensions of swimming micro-organisms 



^ Dombrowski et al] [2004 Zhang et al. 2010), or even 



layers of vibrated granular rods ( Narayan et al] 12007). 
It should be stressed, however, that there have been sug- 
gestions that hydrodynamic interactions may be impor- 



tant in motility assays at high filament density (Schaller 



et 



2011b) 



Each self propelled particle has a position and an ori- 
entation. It also has an individual self -propulsion ve- 
locity of magnitude vq and direction specified by the 
particle's orientation. SPP interact with neighboring 
particles either via a local rule ( [Vicsek et al] |1995[) or 
via physical steric or other interactions (Baskaran and 
2008a|b| |2009| |Leoni and Liverpool| |2010|). 



Marchetti 



Mean field models for the stochastic dynamics of the par- 
ticles can be expressed in terms of the one particle den- 
sity, c(r, z>, t). This measures the probability of finding a 
self-propelled particle with position r and orientation 
at time t. The equation of motion for c(r, z>, t) takes into 
account the interplay of fluctuations (diffusion), interac- 
tions and self-propulsion. These mean-field models are 
valid in the regime of weak interactions, low density and 
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weak density correlations. The kinetic equation for the 
one-particle density can then be solved directly either 
analytically or numerically for specific geometries and 
initial condition (Saintillan, 2012) Alternatively, to de- 



these two approaches have the same structure. The work 
by (Aranson and Tsimring, 2005 2006| ) incorporates, 
however, terms of higher order in the gradients of the 



continuum fields neglected in ( Ahmadi et al. , 2006 



scribe macroscopic behavior, the kinetic theory can also 
be further coarse-grained by projecting on to macroscopic 
fields such as p(r, t) and p(r, t) to obtain the equivalent 
of Eq. ([5| , where the parameters a^p^K^vi^ -^1,2,3 are ex- 
pressed in terms of the microscopic parameters of the sys- 



land Jiilicherl 20001 2003 Liverpool 2003). These mod 



Kruse 



els give rise to contracted states and propagating den- 



sity/polarization waves in both one (Kruse et al, 2001 



2009 Ihle, 2011). A concrete example of this procedure 2006). Finally, a number of other mean- field implemen- 



is given in section V.B.3 below. 



Kruse and Julicher[ 2000) and higher dimensions ( Liv- 
erpool! |2003[ [Liverpool and Marchetti 2005), as well as 



tem ( [Baskaran and Marchetti 2008a 2010 Bertin et al., aster and spiral patterns (Aranson and Tsimring, 2005 



2. Motors and filaments 

Another microscopic realization of active matter is pro- 
vided by suspensions of polar protein filaments cross- 
linked by active cross-links consisting of clusters of molec- 



ular motors (Backouche et al 2006 Mizuno et al 



Nedelec 1998 Nedelec et al 1997 Surrey et al 



2007 



2001). 



The motivation for this in vitro work is to provide an 
understanding of the mechanics of the eukaryotic cell cy- 
toskeleton from the bottom up. Molecular motors are 
proteins that are able to convert stored chemical energy 
into mechanical work by hydrolyzing ATP molecules. 
The mechanical work is done by the motors moving in 
a uni-directional manner along the polar filaments. Par- 
ticular motor proteins are associated with specific po- 
lar filaments, e.g. kinesins 'walk' on microtubules while 
myosins 'walk' on filamentous actin. Since the filaments 
are 'polar, they can be characterized by a position and an 



orientation. Some authors (Ahmadi et al 



2006 



Kruse 



and Jiilicheri 2000, 200^31 12006| [Liverpool! |2003| ILiverpooI 



and Marchettiy ,2005, ) have modeled the motor clusters as 
active cross linkers (see schematic in Fig. 32) capable 
of walking along the filaments and exchange forces and 
torques among filament pairs, hence yielding additional 
(active) contributions to the translational and rotational 
velocities of filaments. Another model has been proposed 



by Aranson and coworkers ( Aranson and Tsimring , 2005 



2006) who have described the filament dynamics via a 



stochastic master equations for polar rigid rods interact- 
ing via instantaneous inelastic "collisions" . Both models 
express the stochastic dynamics of the active system in 
terms of the one particle density of filaments: c(r, z>,t). 
The kinetic equation for c(r, i>, t) takes account of the in- 
terplay of diffusion and active contributions induced by 
the cross linking motor clusters or the inelastic collisions. 
The kinetic equations can again be further coarse-grained 
by projecting on the continuum fields, such as p(r, t) and 
p(r, t) to obtain the equivalent of Eq. [s] where the pa- 
rameters a, /3, ^1, Ai^2,3 can again be related to micro- 
scopic parameters characterizing the filaments and motor 
clusters. Although the parameter values depend on the 
microscopic model, the continuum equations obtained by 



tations of the dynamics of motor filament suspensions, 
some including explicitly motor dynamics or additional 
passive cross linkers, have been used by other authors to 
understand pattern formation in these systems (Lee and 
Kardar! |2001! |Sankararaman et aT] ^2004; Zi ebert et al 



2007 



Ziebert and Zimmermannl 



2004, .2005^ 



In the next sections we briefly outline one procedure 
for going from the fluctuating microscopic dynamics to 
the macroscopic equations of motion. We restrict our 
discussion to dry systems. For wet systems we refer the 



reader to published work ( Baskaran and Marchetti] 



2009 



Leoni and Liverpool! 2010! |Liverpool and Marchetti 



2006 Marchetti and Liverpool 2007). 



B. From stochastic dynamics to macroscopic equations 

We consider a collection of N identical (a simplifica- 
tion) anisotropic particles, each described by a position 
Ynit) and orientation Vnif)^ for particle n. The micro- 
scopic stochastic equations of motion for the positions 
and orientations are given by 

dt^n = Vn(rn, 1>n) + Cn(^) , (74a) 
dtOn = Cc?i(rn, On) X I>n + Sn{t) . (74b) 

The first terms on the rhs are the deterministic contribu- 
tions to the motion arising from both passive interactions 
(e.g., steric repulsion, attractive interactions, etc.) and 
active velocities. Their specific form depends on the par- 
ticular model system considered. The second terms on 



the rhs of Eqs. ( 74 ) are stochastic forces arising from a 



variety of noise sources, including but not limited to ther- 
mal noise. They are assumed to be Gaussian and white, 
with zero mean and correlations 

{^na{t)Uf3{t')) = 2A^f3{l>n) ^nm S{t - t') , (75a) 
{Qna(t)Qm^(t')) = 2Ar Snm S{t - t') . (75b) 

The translational noise correlation tensor is of the form 
A(l>) = A||Z>z> + (^ — z>z>). In a thermal system 
A II = I^ii and Ax = with D\\ and D±^ describing 
Brownian diffusion along the long direction of the par- 
ticle and normal to it, respectively. Also in this case 
Ar = Dr^ with Dr is the rotational diffusion rate. In 
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FIG. 32 (color online) The figure shows a schematic of a pair 
of filaments of length ^ with mid-point positions and orienta- 
tions ri, ij\ and r2, V2 driven by an active cross- link. 



an active system the noise strengths are in general an 
independent quantities. 

For particles with fixed self-propulsion speed v^^ along 
their long axis the deterministic part of the velocities in 



Eqs. (74) has the form 



m 

(76a) 
) , (76b) 

m 

with C(^) = C||^^ + C±(^ ~ ^ friction tensor and 
("^ a rotational friction coefficient. Again, in a Brownian 
system in thermal equilibrium at temperature T friction 
and diffusion (which in this case is also the strength of the 
noise) are related by the Stokes-Einstein relation, D^,^ = 
[C~^]q;/3 ~ ^bT /C,R. In active systems these 

relations do not in general hold and noise and friction 
should be teated as independent. 



The pairwise forces and torques in Eqs. (76) can be 



expressed in term of the total passive and active interac- 
tions in the system. The case of particles propelled by 
internal torques has also been considered in the litera- 



ture, but will not be discussed here (Fily et al. , 2012). 



1. Smoluchowski dynamics 

For simplicity we again restrict ourselves to the case 
of particles with overdamped dynamics. Starting from 



the Langevin equations (74) for the individual filaments. 



standard techniques can be used to derive an equa- 
tion for the one-particle probability distribution func- 
tion (Zwanzig, [2001 ). First one obtains an infinite hi- 
erarchy of equations for the M (< N) particle distribu- 
tion functions, cm (^i , i>i , r2 , i>2 , • • • , tm , i^M , t) , the prob- 
ability of finding M particles with positions and orienta- 



tions {ri, i>i, r2, i>2, . . . , tm, ^^m} at time t, regardless of 
the position and orientation of the other N — M parti- 
cles ( [Baskaran and Marchetti[ [2010 ). A closure approx- 
imation is then required to truncate the hierarchy. This 
is implemented by expressing a higher order distribu- 
tion function in terms of the lower ones. The traditional 
method is to write the 2-particle distribution function as 



a product of two 1-particle densities ( Aranson and Tsim- 



rmgl [20051 [Bertin a/.l [20091 [Kmse and Julicher[ |2000| 



Liverpool[[2003D , as in the familiar molecular chaos ap- 
proximation used to obtain the Boltzmann equation. 

This procedure gives a nonlinear equation for the 1- 
particle density, c(r,i>,t) = ((^(r - rn)^(i> - i>n)), 
where the bracket denotes a trace over all other degrees 
of freedom and an average over the noise, in the form of 
a conservation law, given by 







(77) 



where IZ = i> x do \^ the rotation operator. The transla- 
tional probability current, Jc(r, z>,t), and the rotational 
probability current, ^c(r,2>',t), are given by (Doi and 
Edwardsl |1986D 



V c - A • Vc, Jc = u;c- ArIZc, 



where v and are given in Eqs. (76) and A and A^^ 
are respectively the translational and rotational noise 



strengths introduced in Eqs. (75) 



It should be stressed that this closure scheme may yield 
different physical approximations for different classes 



of microscopic dynamics (Bialke et a/., 2012 Fily and 



Marchetti , 2012 ) . An illustrative example can be found in 



Baskaran and Marchetti (2010) which systematically de- 



rives a kinetic equation for self-propelled hard rods start- 
ing with a microscopic Langevin dynamics that includes 
inertia by first obtaining a Fokker-Planck equation for the 
joint probability distribution of both position/orientation 
and velocities, and finally derive the Smoluchowski limit 
where friction is large relative to inertia at the level of the 
kinetic equation by making a non-thermal assumption on 
the local distribution of velocities. This procedure yields 
a different Smoluchowski equation from that obtained by 
simply taking the overdamped limit at the level of the 
Langevin equation. While the structure of the hydro- 
dynamic equations obtained by coarse-graining the one- 
particle kinetic equation is the same in both cases (as 
it is dictated by symmetry considerations), the values of 
the parameters in the continuum equations and particu- 
larly their dependence on vq and noise strength depend 
on whether the overdamped limit is taken right at the 
outset, i.e., at the level of the microscopic dynamics, or 
at the level of the kinetic equation. Much work remains 
to be done to understand this subtle point and therefore 
we will not discuss it further in this review. A further in- 
teresting open question is how to systematically obtain a 
stochastic equation for the 1-particle density which prop- 
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FIG. 33 (color online) An example of self-propelled particles: 
cartoon of a motility assay with filaments of length I with 
positions and orientations ri,i>i and r2,i>'2 driven by motors 
tethered to a plane. 



erly includes the noisy dynamics of the density fluctua- 
tions for active systems (E)eanj [l996). 

We leave some of these very interesting open questions 
for the reader and discuss below how to use the Smolu- 



chowski dynamics as described by Eqs. (77) in the deriva- 
tion of hydrodynamic equations. 



2. From Smoluchowski to hydrodynamics 



One approach, explored extensively by Saintillan, 



Shelley and collaborators (Saintillan and Shelley , 2007 
2008a|b ), is to solve directly the Smoluchowski equa- 



tions, either analytically or numerically. This work has 
been used to investigate the stability of both aligned and 
isotropic suspensions of active particles with hydrody- 
namic interactions (see also below) and has revealed a 
rich dynamics with strong density fluctuations. This ap- 
proach will not be discussed further here. A recent review 



can be found in Saintillan (2012). 



Here we focus instead on obtaining the description of 
the dynamics of the system in terms of a few macroscopic 
fields introduced phenomenologically in the first part of 
this review. A crucial assumption in deriving this contin- 
uum or hydrodynamic theory is the choice of the contin- 
uum fields as those whose fluctuations are long lived on 
large length scales. They include fields associated with 
conserved quantities and possible broken symmetries of 
the system. In a collection of active particles with no 
momentum conservation, the only conserved quantity is 
the number of particles and hence the density, defined 



in Eq. (la) is a slow variable. In addition, to allow for 



the possibility of broken orientational order, with either 
polar or nematic symmetry, we consider the dynamics of 
a polarization field (Eq. ( |lb[ )) and an alignment tensor 
(Eq. ([21])). The fields are defined as moments of the one 



particle distribution function as. 

p(r,t) = J di> c(r,z>,t) , 
p(r,t)p(r,t) = J dO c(r,z>,t) , 
p(r,t)Q(r,t) = J Q(i>) c(r,i>,t) . 



(78a) 
(78b) 
(78c) 



with Q{u) = z>z> — ^1 and d the system's dimension- 
ality. In general one could write an exact expansion of 
c(r, z>, t) in terms of all its moments and transform the 
Smoluchowski equation into an infinite hierarchy of equa- 
tions for the moments. The hydrodynamic description is 
obtained by assuming that higher order moments relax 
quickly on the time scales of interest and by truncating 
this expansion to include only conserved quantities and 
order parameter fields. Equivalent ly, considering for sim- 
plicity the case of d = 2, one can simply write 

c(r, !>, t) = { 1 + 2p(r, t)-u + 4Q(r, : Q(i>) } , 

(79) 

insert this ansatz into the Smoluchowski equation. 



Eq. (77), and obtain the hydrodynamic equations for 



filament concentration, polarization and alignment ten- 
sor. For the details of the calculation, which also in- 
volves using a small gradient expansion for the fila- 
ment probability distribution and evaluating angular av- 
erages, we refer the reader to the literat ure (|Ahmadi] 



et a/.[[2QQ6{|Baskaran and Marchetti"l|2008a[[Bertin et~ 



2009 



Ihle| (20111 [Liverpool! |2003D . The result is a set 



of coupled equations for the hydrodynamic variables, 
p(r, t), p(r, t), Q(r, t) that has the form given in section 
[ll| with parameter values that are calculated explicitly 
in terms of microscopic properties of the system. 

In particular, one finds that the continuum equations 
governing the dynamics of motor-filament systems and 
collections of self-propelled particles have the same struc- 
ture, but important qualitative differences arise in the 
hydrodynamic parameters. In models of SPP the acti- 
vating internal processes drive each individual unit (as 
for instance in bacterial suspensions or actin filaments in 
motility assays, where the filaments are driven by the ac- 
tion of a carpet of myosins tethered to a plane) and active 
contributions to the dynamics arise even at the single- 
particle level. In cell extracts of cytoskeletal filaments 
and associated motor proteins activity arises from inter- 
actions among the filaments mediated by motor clusters 
that act as active crosslinkers, exchanging forces among 
the filaments. As a result, the various active parameters 
in the hydrodynamic equations are proportional to the 
square of the density of filaments, resulting in different 
behavior at large scales. For this reason motor- filament 
suspensions have also been referred to as systems of "mu- 
tually propelled particles" (MPP) QGiomi aH|2012D , to 
distinguish them from SPPs. 
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3. An example: derivation of continuur 
Vicsek-type particles 



equations for aligning 



To illustrate the method, we show in this section 
the derivation of the continuum equations for a simple 
model of self-propelled point particles on a substrate in 
two dimensions with a polar angular interaction that 
tends to align particles as in the Vicsek model. We 
consider N particles with polarity defined by an axis 
z>n = (cos 7 sin ^n)- The dynamics is governed by 
Eqs. (I74|) that now take the explicit form 



dV 



m ^ 



(80a) 



(80b) 



where vq is the fixed self-propulsion speed and we have 
taken = (Safd- The noise is Gaussian and white, 
with zero mean, as given in Eqs. (75) and Aa/s = ASa/s- 
Forces and torques in Eqs. ( |80a| ) are expressed here as 
derivatives of a pair potential, given by C^^V{xn,Xm) = 
— ;^^0(i? — \rn — Tm\) cos{Om " ^n), that tcuds to align 
particles of the same polarity, with Q{x) the Heavside 
step function and R the range of the interaction. Here 
we use the shorthand Xn = (r^,6>n) and 7 is the inter- 
action strength with dimensions of inverse time. As we 
are interested in a hydrodynamic model that describes 
spatial variations on length scales large compared to the 
range R of the interaction, we assume the interaction to 
be local. This polar aligning interaction was used re- 
cently in Farrell et al. (2012) in a related model that 



highlights the role of caging effects. As outlined in sec- 
tion V.B.I , one can then use standard methods (Zwanzigj 



2001 ) to transform the coupled Langevin equations into a 



Smoluchowski equation for the one-particle distribution 
c(r, ^, t), given by 



a^c+^o^^rVic = AV?c+Ai^S^^c+^Vrc(xi,t) / V iV{xi,X2)c{x2,t)^ARdQ^c^^de^c{xi,t) j de^V{xi,X2)c{x2,t) 

(81) 



To obtain hydrodynamic equations, we now proceed as 
outlined in section [V.B . 2| and transform the kinetic equa- 
tion into a hierarchy of equations for the angular mo- 
ments of the one-particle probability density, c. In two 
dimensions the moments are most simply written by in- 
troducing an angular Fourier transform ( Bertin et aL\ 
|2009 ) in terms of the Fourier components fk{T,t) = 

comparing with Eqs. ( [TS] ) it is 
easy to see that /o = fi = -\- iwy^ with w = pp the 
polarization density, and the real and imaginary parts 
of /2 are proportional to the two independent compo- 
nents of the alignment tensor, Q. Using 27rc(r, 6>,t) = 
Xlfe /fce"*^^ and retaining for simplicity only terms up to 
linear order in the gradients, we obtain a hierarchy of 



equations 

9tfk + -^dxifk+i + fk-i) + -^dy{fk^i - fk-i) 



(82) 



where all discrete sums run from —00 to -hoc and Vn 



smf 



in {Sq^i 



We consider equations 



for /o and /i, assume that /2 is a fast variable, so the 
dtf2 — 0, and neglect all higher order Fourier compo- 
nents, i. e., assume //c = for A: > 3. Eliminating /2 
in favor of /o and /i, it is then straightforward to ob- 



tain (Farreh et a/., 2012) 



3^07 
16Ap 



(W • V)W : 



jp- Ar]w 



dtp- 
8A 



^0 



R 



Vp- 



5^07 
32A;, 



5^07 
16A^ 



w(V-w) + 0(V^) 



(83a) 
(83b) 



Microscopic derivations of continuum equations of the to transform the equation for w into an equation for p, 
type presented here naturally yield an equation for the 
polarization density w = pp rather than an equation 
for the order parameter field p as given in Eq. (|5|. Of 
course it is straightforward to use the density equation 
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with the result, 

dtP+Xi{p ■ V)p = - [a{p) + /3|p|2] p - vi ■ V/j 



with 



'V|pp + A2p(V-p) + 0(V2) 



a{p) = Ar- --yp , 



/? = 
Ai = 
Ao = —A: 



vi 



2p V 



P^U-- 1 



(84) 

(85a) 
(85b) 
(85c) 
(85d) 
pp. (85e) 



Eq. ( 84 ) has the same structure as Eq. ([5| , although noise 
has been neglected here. The various coefficients are ex- 
pressed in terms of microscopic parameters and are in 
general found to depend on density and order parameter. 
In particular, the coefficient vi of the Vp term is now a 
tensor, describing anisotropic pressure gradients that can 
play a role in the ordered state. These effects were ne- 
glected for simplicity in the phenomenological model. We 
also stress that there is an important difference between 
the parameter values obtained in the present model and 



those obtained by Baskaran and Marchetti (2008a) and 



Bertin et al. ( 2009 ) as in both the latter models the effec- 



tive interaction strength (here 7) depends linearly on vq. 
As a result, both Baskaran and Marchetti (2008al) and 
Bertin et al. ( 2009 ) find that A^ ^ v^^ while here A^ ^ vq. 



This difference arises because in lBaskaran and MarchettH 



(2008a) and Bertin et al. (2009) the authors systemati- 



cally describe binary collisions, yielding an instantaneous 
change of direction with probability one in a small time 
interval, while the model of Farrell et al. (2012) presented 



here (and chosen for illustrative purpose because of its 
simplicity) considers a continuous evolution of the posi- 
tion due to the application of a force in the Langevin 
equation. These differences may of course be important 
for large values of v^). 

This example shows how the derivation of hydrody- 
namics from a microscopic model yields explicit values 
for the various parameters in the continuum equations in 
terms of microscopic parameters. Of course, as stressed 
above, these values are model dependent. In addition, 
there is a price to be paid in that in the derivation one 
has made two important assumptions. The ffist is the as- 
sumption of low density, that allows us to replace the two- 
particle probability distribution by the product of two 
one-particle distribution functions in the kinetic equa- 
tion. The second is effectively an assumption of weak 
interaction that enters in the moment closure approxi- 
mation. The two assumptions are not independent and 



essentially amount to the so-called assumption of "molec- 
ular chaos" used for instance in the derivation of the fa- 
miliar Boltzmann equation. 



4. Hydrodynamic interactions 

The examples presented above ignore the momentum 
exchange between the active particles and the solvent 
that one would have for instance in suspensions of swim- 
ming organisms. This physical effect is important for 
understanding the properties of a wide variety of experi- 
ments and should be included for a complete microscopic 
description of active matter. 

The coupling to the solvent has been treated in the 
literature in at least two ways. 

One approach is to include an additional equation 
for the solvent velocity field v(r,t) which is coupled 
to the microscopic equations of motion ([Liverpool and] 



Marchetti 2006 Marchetti and Liverpool 2007) or to the 



corresponding Smoluchowski equation (Pahlavan anci] 
Saintillan[ |2011[ [Saintillan and Shelley[ |2008a[ ). The^c- 
tive particles are convected and rotated by the local fluid 
flow while also imparting additional forces (stresses) onto 
the fluid. This gives rise to a coupled equation of motion 
for the concentration of filaments c(r, z>, t) and the ve- 



locity field v(r, t) which can be solved directly (Pahlavan 
and Saint illan[ [2011 Saintillan and Shelley[ |2008a ) or 



coarse-grained as described above to give coupled equa- 
tions for p(r, t), p(r, t), Q(r, t), v(r, t). For the case where 
the suspension is incompressible (p = constant), the mo- 
mentum (v) equation has the form p(9t + v- V)v = V-tr. 
This equation may also be often treated in the Stokes ap- 
proximation, V • (T = 0, appropriate for systems at low 
Reynolds number. The active forces due to the swim- 
mers are incorporated as an additional (active) contribu- 
tion to the stress tensor, so that (Ja(3 = cr^^ + where 
the passive contribution cr^^ has the form obtained in 



equilibrium liquid crystals (de Gennes and Prost, 1993). 



The active contribution is evaluated by noting that 
in the absence of external body forces, the force distri- 
bution exerted by the active particles on the fluid can 
be written as a multipole expansion with lowest order 
non vanishing term being a dipole. The active stress 
<j^^ can then be identified with the active reactive fluxes 



proportional to ("A/i (see Eq. 41). To leading order in 



a gradient expansion it contains two additive contribu- 
tions. The ffist, cr^^ (X CA/ig'a/3, has nematic symmetry 
and is present in both polar and nematic systems, where 
it is written as <j^^ oc (AjiiQa/s- In addition, for po- 
lar systems only, the active stress also contains terms 



a^|oc ('AiJ.daP(3 -\-C'^f^d^Pa (Liverpool and Marchetti 



2006[ [Marchetti and LiverpooHTMPfl ). When the hydro- 



dynamic equations are written on the basis of the en- 



tropy production formulation discussed in section [III.B[ 
this term is discarded as one of higher order in the driv- 
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ing forces. It is also of higher order in the gradients as 
compared to a^^, although only of first order, rather than 
quadratic, in the polarization field. Finally, it is the lead- 
ing non- vanishing contribution to the active stress that 
has polar (as opposed to nematic) symmetry and it does 
play a role in controlling the onset of oscillatory states in 



these systems ( Giomi and Marchetti[ 2012 Giomi et al. 
[2OO8). This shows that a microscopic derivation is needed 
to go beyond linear hydrodynamics. 

Alternatively, one may integrate out the solvent veloc- 
ity field to generate effective hydrodynamic interactions 
at the two-body level which give rise to additional long- 
range pair-wise contribution to the deterministic veloci- 
ties and angular velocities v^, u?^ in Eqs. ( 74 ) ( Baskaran 



and Marchetti 



2009 Leoni and Liverpool 



2010). These 



interactions arise from the forces generated on the fluid 
by the active particles. The local force distribution due 
to the active elements can be expanded in a multipole 
expansion. The lowest non vanishing term in this ex- 
pansion is a dipole, which also gives the longest range 
interactions. The simplest static models of swimming 
organisms then are force dipoles which can be character- 
ized as contractile or extensile depending on the direction 
of the forces making up the dipole. Coarse- graining the 
equations leads to a set of coupled integral equations for 
p(r, t), p(r, t), Q(r, t) which in extended (large) domains 
can be used to obtain generalized linear hydrodynamic 
modes whose behavior is identical to that obtained when 
keeping an explicit velocity field. Hydrodynamic equa- 
tions for collections of static force dipoles have been ob- 



tained (Baskaran and Marchetti 2009) which show all 



the phenomenology described before. Swimming objects 
however are dynamic, undergoing an internal cyclical mo- 
tion that leads to self-propulsion. The collective behavior 
of simplified dynamic models of swimmers with such in- 
ternal cycles has also been studied and used to generate 



effective hydrodynamic equations (Leoni and Liverpool 



2010). Averaged over an internal cycle such dynamic 



models have an average force distribution that also has a 
multipolar expansion whose lowest term is generically a 
force dipole. The hydrodynamic equations obtained are 
then similar to those obtained for static dipoles. How- 
ever by tuning internal parameters one can also study 
self-propelled swimmers whose averaged force distribu- 
tion starts at quadrupolar order. Here the symmetry 
broken between contractile and extensile objects is re- 



stored and only nematic phases are possible (Leoni and 



Liverpool 2010, 2012) 



C. Current Status of Microscopic Theories of Active 
Matter 



Perhaps the most important open question in deriving 
active theories from microscopic or mesoscopic models 
concerns understanding the nature of the noise. Ther- 



mal noise is often negligible in active systems at low fre- 
quency, but noise is nonetheless ubiquitous and plays a 
crucial role in controlling the large scale behavior, as it 
can both act to destroy large scale coherence or drive 
synchronization, depending on the specific situation. So 
far practically all work on active systems has either ne- 
glected noise or modeled it as a Gaussian, white random 
force, akin to thermal noise, but of unknown strength. In 
general one expects the noise to depend on activity. Mul- 
tiplicative or non-Markovian random forces could also 
be at play in active systems. For instance, the multi- 



plicative character of the noise introduces ( Mishra et al 
|2012) singular features into the fluctuation- driven phase- 



separation ( Mishra and Ramaswamy , 2006 ) that charac- 
terizes a stable active nematic. Microscopic models are 
needed that attempt to take into account the stochas- 
tic nature of active forces and to derive a coarse-grained 
model where the effective noise amplitude is expressed 
in terms of local nonequilibrium processes that may lead 
to temporal or spatial correlations at large scales. An 
example is the work of Lacoste and Laii| ( |2005 ) that 
couples the shot noise associated with the on/off switch- 
ing of energy-dissipating pumps in active membranes to 
the membrane curvature and demonstrates that pump 
stochasticity plays a crucial role in controlling super dif- 
fusive behavior in the membrane. 

Of great interest is also the detailed understanding of 
hydrodynamic interactions in bacterial suspensions and 
their role in controlling the large-scale behavior. Recent 
experiments have probed for the first time the flow fields 
induced by swimming unicellular organisms. Measure- 
ments on Chlamydomonas have revealed qualitative dif- 
ferences as compared to the puller stresslet configura- 



tion that had been used in the literature ( Drescher et al. 



2010) and have further shown complex time- dependent 
oscillatory patterns. In contrast it was found that the 
flow field of E.coli is well described by a pusher stresslet, 
but its strength is very small and it is washed out by 



rotational diffusion of the swimming direction (Drescher 



et 



2011j). These results have important implications 



for the behavior of microorganisms near surfaces and 
open the way to new quantitative investigations of the 
role of hydrodynamic interactions. 

Finally, an area that is receiving increasing attention 
is the study of active matter at high density, where ac- 
tive glassy or solid states may emerge. As mentioned in 
section |II.C| the collective dynamics of confluent layers 
of epithelial cells has been likened to that of glassy and 
supercooled systems. Mesoscopic models of interacting 



cell layers and tissues are beginning to emerge (Henkes 



et al, 2011), although much more work remains to be 



done to understand the complex interplay of contractile 
stresses and substrate adhesion in controlling the build 
up of cellular stresses in collective cell migration. 
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VI. CONCLUSIONS, 
DIRECTIONS 



OUTLOOK AND FUTURE 



2001 Reichenbach et a/., 2006). We note, however, that 



In this review we have discussed two dimensional and 
three dimensional active systems, that are maintained 
out of equilibrium by a permanent energy consumption 
that takes place locally in each active unit. In these 
systems local polarity, when present, yields spontaneous 
motion of the polar entities with respect to the center of 
mass motion. Typically, models of active matter describe 
the collective behavior of systems like fish shoals, birds 
fiock or animals herds as well as vibrated granular mat- 
ter, bacterial colonies or the cellular cytoskeleton. From a 
spatial symmetry point of view these systems are not dif- 
ferent from ferroelectric and nematic liquid crystals. The 
fundamental difference stems from the fact that active 
systems consume and dissipate energy at all times. This 
feature has a number of profound consequences which we 
have addressed in this review: 

• Polar fiocks on a substrate, e.g., animal herds, ker- 
atocytes on a glass slide, display long-range order 
in two dimensions. This behavior is fundamentally 
different form that of equilibrium systems that are 
bound to obey the Mermin- Wagner theorem and 
can at best exhibit quasi-long range order; 

• nematic as well as polar ordered phases systems 
on a substrate exhibit giant density fiuct nations, 
breaking the familiar equilibrium scaling of 
number fiuctuations in subregions with N particles 
on average; 

• the uniform ordered state of bulk momentum- 
conserving systems of polar and nematic symmetry 
is generically unstable in the Stokesian regime; 

• the end result of this instability for the typical ac- 
tive suspension appears to be turbulence at low 



gemuth 



Reynolds number (Dombrowski et al. 2004 



Wol- 



2008) driven by the competition between 
forcing by active stress and relaxation by orienta- 
tional diffusion. 

• all these systems can support a new type of sound- 
like propagating waves that can have different prop- 
agation laws in opposite directions due to polarity 

We have described how macroscopic hydrodynamic 
equations for active systems can be obtained either from 
symmetry arguments, from generalized thermodynam- 
ics close to equilibrium, or from microscopic models. 
In each case we have presented the simplest descrip- 
tion that highlights the generic features uniquely associ- 
ated with activity. Our analysis is restricted to two and 
three dimensions and we have deliberately left aside one- 
dimensional and quasi one-dimensional systems, which 



the mechanisms leading to the onset of traffic jams are 
probably related to the giant density fluctuations pre- 
dicted in active systems in two dimensions [see also ( Fily 



et al 2012 Yang et al. 2010)] 



The rewarding feature is that all these unusual features 
have either been observed experimentally or conflrmed in 
careful numerical simulations. It is very diflicult to es- 
tablish the existence of long range order experimentally, 
but simulations with a large number of active particles 
show unambiguously its existence. This is an important 
result that the statistical physics community was not in- 
clined to believe. Giant density fluctuations have now 
been observed in simulations (Chate et al, 2008a), ex- 



perimentally in vibrated granular media (Narayan et al. 



2007) and in bacterial suspensions (Zhang et a/., 2010). 
These observations have conflrmed the theoretically pre- 
dicted scaling. We note, however, that this scaling has 
recently been questioned by Toner (2012). These large 



fluctuations are also present in the cellular cytoskeleton 



and their role is starting to be recognized ( Gowrishankar 



and Rao, 2012) and their potential physiological rele- 



vance is beginning to be explored (jChaudhuri et al. , 2011 



Goswami et al. 2008 Gowrishankar et al. 2012). 



The long wavelength instability of momentum conserv- 
ing systems is conflrmed by numerical simulation and in 
fact can give rise to "low Reynolds number turbulence" 
in some instances (Dombrowski et al. 20041). It is a 



good paradigm for explaining bacterial swirls but also 
the rotating microtubule spindles observed both in vivo 



and in vitro (Nedelec[ |1998| [Nedelec et al.\ |1997[ ). Fur- 
thermore, the existence of "low Reynolds number" waves 
have also been predicted in crystals moving without in- 
ertia through a dissipative medium, such as sediment- 



ing colloids (Lahiri and Ramaswamy, 


1997) or drifting 


Abrikosov lattices ( 


Balents et al. 


1998 


Ling et al. , 1998 


Simha and Ramaswamy, 1999 


), and also observed and 



understood in microfluidic drifting drop arrays (Beatus 



et al. 2007 2006). Such waves propagating in only one 
direction have been observed in cell lamellipodia, as ex- 
pected from active gel theory (Giannone et al. 2004). 



Eventually, one may be able to discuss with order of 
magnitude accuracy phenomena such as cell wound heal- 



ing (Petitjean et al. 2010), cell division (Salbreux et al. 



2009), and ceh oscillations (Salbreux et al.\ 2007). One 



have already been investigated extensively (Helbing 



can thus say that we currently have a good qualitative 
understanding of active systems. What do we need to get 
to the next stage, that is to get to real quantitative un- 
derstanding? One ingredient which we have omitted on 
purpose up to now is signaling: biochemical signaling for 
cell biology or bacterial colonies external factors like food, 
smell or sun in flocks, shoals and herds. Such factors have 
been discussed in bacterial colonies and shown to give rise 
to beautifully ordered structures (Gates et al. 2010). A 
direction that looks particularly promising, is coupling of 
active gel physics and chemical reactions in the context 
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of cells and tissues, resulting for instance, in increasing 
considerably the range of existence of Turing structures 



(Bois et a/., 2011 Giomi et al 2010, 2012), and in the 



spatiotemporal control and enhancement of chemical re- 



action rates in signaling and sorting (Chaudhuri et al. 



2011). A complete understanding of cell biology will ob- 
viously require including detailed biochemical signaling 
networks. However this cannot be achieved without an 
extensive experimental input. Symmetry considerations 
and conservation laws are no longer sufficient and this 
task requires the cooperation of several disciplines. It is 
worth the effort though, since it could be extended to 
tissue dynamics and developmental biology. 

Another natural extension in the case of herds, flocks 
and shoal is the investigation of a putative leader 

120121). Up to 



role (Couzin et al. 2011 



Leonard et 



now, completely "democratic" rules have been assumed. 
Within this approach, a leader could appear like a delta 
function, and its role involve the response function of 
the collection of individuals. In another promising direc- 



tion, Guttal and Couzin (2010) treat the parameters in 



a Vicsek-style model of flocking and migration as heri- 
table, selectable attributes, and study their dynamics on 
evolutionary timescales, flnding a remarkable range of 
evolutionarily stable strategies including coexistence of 
distinct behaviors in a herd. 

Another frontier concerns active quantum systems. 
A beautiful classical mimic has been introduced re- 



cently ( ^Couder and Fort 2006| ) and some Al-Ga-Se het- 
erostructures give rise to zero resistivity behavior de- 
scribed by equations very similar to those of Toner- Tu 
(Alicea et al.^ |2005| ) . A general investigation of the type 
we present here would be very useful. 

Eventually working with real systems is necessary and 
important, but it is often difficult. Artiflcial systems 
such as provided by layers of vibrated granular mat- 



ter ( 


Deseigne et al. 2010 Kudrolli et al. , 


2008 


Ku- 


mar et al.\ 2011 Narayan et al. 


2007), noisy walkers 


(Kumar et al. 2008), artiflcial swimmers ( 


Bart olo and 


Lauga 


2010 Dreyfus et al. 


2005 ) and colloids propelled 


by catalytic reactions (Gibbs and Zhao| 201 1| Golesta- 


nianl 20091 Golestanian et al.\ 2005| Howse et al.\ |2007| 


Palacci et al. , 2010 Paxton et al. 


, [2004|) are very use- 



ful but are often restricted to two dimensions (granular 
layers) or have been studied only at relatively low den- 
sity (artiflcial self-propelled particles). It would be nice 
to have three dimensional active systems driven by light 
or chemical reactions (Buttinoni et al. 2011). Attempts 



have been made (Prost, 2010) to raise interest among 



chemists, but so far progress towards realizing artiflcial 
active gels has been limited. Such systems would allow 
one to vary parameters over a wide range and in a con- 
trolled manner and would therefore provide extremely 
valuable artiflcial realizations of active matter. 
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